题目描述


输入描述:

输出描述:

示例1
输入
3
输出
5

分析:

1.这个题其实考的是一个莫比乌斯反演题,但是由于我知识储备不够,没有看出来,题目给的范围可以瞎搞一下,所以下面容斥可以过。
2.转换一下就是一道经典的反演题,参考

AC代码:

容斥:

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long ll;
const int M=1e5+10;
int f[M];
int cal(int x){int num=0;while(x){num+=x%10;x/=10;}return num;
}
int main(){int n;cin>>n;ll ans=0;for(int i=2;i<=n;i++){int m=i,kk=0;for(int j=2;j*j<=n;j++){///唯一分解定理,求其素因子if(m%j==0){f[kk++]=j;while(m%j==0)m/=j;}}if(m>1) f[kk++]=m;int lim=1<<kk,sum=0;for(int j=1;j<lim;j++){///枚举哪些质因子被拿出来了int x=j,su=0,sm=1;for(int k=0;k<kk;k++){///找出哪些素因子被拿出if(x&(1<<k)){su++;sm*=f[k];}}if(su&1) sum+=n/sm-i/sm;///规定,当素因子拿出奇数个,则加上,在(i~n)里存在该质因子数的个数。else sum-=n/sm-i/sm;///当为偶数时,此时之前已经拿出,重复,需减掉。}///sum为不互质的数;///n-i,为总数,原为找寻的for(j=i+1,j<=n;j++)中(i,j)互质的数ans+=(n-i-sum)*cal(i);}ans+=n;///没有枚举1,因为全部会被标记,所以+1cout<<ans<<endl;return 0;
}

莫比乌斯反演:

#include<stdio.h>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int M=1e5+10;
int prime[M],mu[M],f[M];
int cnt;
bool book[M];
void sieve(int N)
{mu[1]=1;for(int i=2; i<=N; i++){if(!book[i])prime[cnt++]=i,mu[i]=-1;for(int j=0; prime[j]<=N/i; j++){book[i*prime[j]]=true;if(i%prime[j]==0)break;mu[i*prime[j]]=-mu[i];}}for(int i = 0; i <= N; i++){int x = i;while(x){f[i]+=x%10;x/=10;}}
}
int main()
{sieve(M - 1);int n;scanf("%d",&n);ll ans=0;for(int i=1; i<=n; i++){if(!mu[i])continue;ll res=0;for(int j=1; j<=n/i; j++)res+=(ll)f[j*i]*(n/i-j+1);//1LL,长整型1,和*1.0一样,改变数据类型用;ans+=mu[i]*res;}printf("%lld\n",ans);return 0;
}

莫比乌斯反演/容斥 +2020ICPC 江西省大学生程序设计竞赛 A Simple Math Problem相关推荐

  1. BZOJ2301:[HAOI2011]Problem b(莫比乌斯反演,容斥)

    Description 对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函数为x和y的最大公约数. Input 第一行一个整数 ...

  2. BZOJ 2301 - Problem b(莫比乌斯反演+容斥)

    题目链接 https://cn.vjudge.net/problem/HYSBZ-2301 [题意] 对于给出的 n 个询问,每次求有多少个数对(x,y)(x,y) ,满足 a≤x≤b,c≤y≤da≤ ...

  3. 2021 ICPC 江西省大学生程序设计竞赛(热身赛)【完结】

    热身赛,都是练手速的水题 目录 Abnormal unity Dripping Water Wears Through a Stone Follow me Strength team Abnormal ...

  4. 2021 ICPC 江西省大学生程序设计竞赛(正式赛)

    K.Many Littles Make a Mickle(签到题) 题目大意:有T组测试样例,每一组测试样例给一个n,一个m,从1-n层,第i层有i*i个房间,每个房间居住m个人,问一共能住多少人? ...

  5. 2021年中国大学生程序设计竞赛 女生专场 - 热身赛 Problem C. 口算训练(质因子分解)

    题目分析:判断 xxx 是 yyy 的倍数,等价于质因子分解后, yyy 中的每个质因子的出现次数都小于等于其在 xxx 中的出现次数. 那么对于每次询问 [l,r,d][l,r,d][l,r,d], ...

  6. Comet OJ - Contest #8 E神奇函数(莫比乌斯函数容斥)

    Comet OJ - Contest #8 E神奇函数(莫比乌斯函数容斥) 题目大意 定义d(x)d(x)d(x)的值为x的最小素因子,定义 f(x)={1x=1d(x)f(xd2(x))x>1 ...

  7. “亚信科技杯”南邮第七届大学生程序设计竞赛之网络预赛 (部分题解)

    "亚信科技杯"南邮第七届大学生程序设计竞赛之网络预赛 比赛链接http://acm.njupt.edu.cn/acmhome/contest.do?&method=cont ...

  8. 2018中国大学生程序设计竞赛-网络选拔赛题解

    以下所有AC题解程序来自"仙客传奇"团队. A. Buy and Resell AC的C++语言程序: #include<iostream> #include<c ...

  9. “统信杯” 第十七届黑龙江省大学生程序设计竞赛

    "统信杯" 第十七届黑龙江省大学生程序设计竞赛 前言 题目 A. Bookshelf Filling C. Tree Division D. Collision Detector ...

最新文章

  1. 2022-2028中国工作服市场现状及未来发展趋势
  2. ccf压缩编码java_压缩编码(CCF)——动态规划
  3. Sql Server2005一則怪現象
  4. python numpy库安装-Python Numpy库安装与基本操作示例
  5. python modbus类封装_Python 中引入一个文件,模块的概念
  6. 017_Cookie
  7. ArrayBlockingQueue原理分析-dequeue方法
  8. Java对象的引用类型
  9. mysql properties文件路径_读取web项目properties文件路径 解决tomcat服务器找不到properties路径问题...
  10. 检查键是否存在于JavaScript对象中?
  11. 实战CSS:小米商城静态实现
  12. 海康摄像头不记得登陆密码了怎么办?2020年联系人工客服最新版解答方案
  13. 使用计算机时按什么键打开,计算机进入BIOS的快捷键是什么|打开计算机以进入BIOS设置时按下的键...
  14. 数字逻辑基础实验二—时序逻辑电路的设计
  15. Android开发艺术探索完结篇——天道酬勤
  16. 字典大全(修改,添加,删除)所有遍历
  17. C# 在图片上写字
  18. Substance Painter插件添加
  19. 三公到底指太师太傅太保还是…
  20. vue实现点击标题,显示隐藏兄弟元素(类似电脑打印设备显示打印设备)

热门文章

  1. linux之fdisk查看分区和mkfs.ext3删除分区和mount挂载和e2label添加卷标使用总结
  2. 用非递归方式实现二叉树先序便利
  3. Android之android studio如何解决‘:app:packageDebug‘.(Duplicate files copied in APK META-INF/DEPENDENCIES)
  4. python process 函数_Python Process创建进程的2种方法详解
  5. iphone换机数据迁移_iPhone迁移数据到Android(相册与短信)
  6. 软件测试实验报告下载 实验一到实验五
  7. 第1章 数据分析概述
  8. 男女共厕,考验男生心理素质的时候到了......
  9. 和哪个专业的男生谈恋爱最惨?
  10. python怎么获取lol皮肤名称_LOL手游免费皮肤获得方法 LOL手游皮肤怎么获得