UESTC 618 无平方因子数 (容斥 + 莫比乌斯反演)
来源:互联网 发布:js获取当前元素index 编辑:程序博客网 时间:2024/06/06 20:18
无平方因子数
Time Limit: 4000/2000MS (Java/Others)
Memory Limit: 65535/65535KB (Java/Others)
Submit Status
无平方因子数即对于任意一个素数p,p2都不会整除那个数,如1 , 5=5 , 15=3×5都是无平方因子数,而20=22×5不是。现在给定一个n (1≤n<1012) ,求区间[1,n]中无平方因子数的个数。
Input
第一行有个整数T,代表数据组数(T≤10)
接下来有T行,每行有个整数n (1≤n<1012)
Output
输出T行,每行输出一个整数代表区间[1,n]内的无平方因子数的个数。
Sample Input
3
1
10
30
Sample Output
7
19
Source
UESTC Training for Math
题目链接:http://acm.uestc.edu.cn/#/problem/show/618
题目分析:又是无平方因子数,比BZOJ那题简单很多,直接算就行了,参照BZOJ 2440
#include <cstdio>#include <cstring>#include <algorithm>#define ll long longusing namespace std;int const MAX = 1e6 + 5;int mob[MAX], p[MAX];bool prime[MAX];void Mobius(){ int pnum = 0; memset(prime, true, sizeof(prime)); mob[1] = 1; for(int i = 2; i < MAX; i++) { if(prime[i]) { p[pnum ++] = i; mob[i] = -1; } for(int j = 0; j < pnum && i * p[j] < MAX; j++) { prime[i * p[j]] = false; if(i % p[j] == 0) { mob[i * p[j]] = 0; break; } mob[i * p[j]] = -mob[i]; } }}ll cal(ll n){ ll cnt = 0; for(ll i = 1; i * i <= n; i++) cnt += (ll) mob[i] * (n / (i * i)); return cnt;}int main(){ Mobius(); int T; scanf("%d", &T); while(T --) { ll n; scanf("%lld", &n); printf("%lld\n", cal(n)); }}
0 0
- UESTC 618 无平方因子数 (容斥 + 莫比乌斯反演)
- UESTC 618 无平方因子数 ( 莫比乌斯)
- UESTC 1720 无平方因子数(数论,容斥)
- BZOJ 2440 完全平方数 (容斥+莫比乌斯反演+二分)
- 【bzoj2440】[中山市选2011]完全平方数 莫比乌斯反演+容斥原理
- BZOJ_P2440 [中山市选2011]完全平方数(数论+莫比乌斯反演+容斥原理)
- 【BZOJ2440】完全平方数,莫比乌斯反演+二分答案+容斥思想
- BZOJ 2440 完全平方数(莫比乌斯反演+容斥原理+二分)
- 【BZOJ 2440】[中山市选2011]完全平方数 莫比乌斯反演+容斥
- BZOJ 2440: [中山市选2011]完全平方数 莫比乌斯反演 +容斥原理
- BZOJ 2440 中山市选2011 完全平方数 二分答案+容斥原理+莫比乌斯反演
- BZOJ-2440 中山市选2011 完全平方数 二分查找 + 莫比乌斯反演 + 容斥原理
- bzoj 2440: [中山市选2011]完全平方数(二分,容斥,莫比乌斯反演)
- 无平方因子数 UESTC 1720
- [容斥+莫比乌斯]bzoj2440_完全平方数
- uestc 1720 square-free integer(无平方因子数)&nyoj 580square-free integer(容斥,数论)
- 【莫比乌斯反演】[BZOJ2440]完全平方数
- bzoj2440: 完全平方数【莫比乌斯反演】
- ScrollView嵌套ListView,.measure(0,0);空指针NullException
- char *a 与char a[] 的区别
- 【Qt OpenGL教程】11:旗帜效果(飘动的纹理)
- 零售数据分析要点
- git rev tips
- UESTC 618 无平方因子数 (容斥 + 莫比乌斯反演)
- rsync关于adler-32的描述
- iOS控件封装-(UIScrollView)图片轮播器
- iOS 的 APP 如何适应 iPhone 5s/6/6Plus 三种屏幕的尺寸?
- Exponentiation
- 更新最后一条记录
- gcc编译运行一个或多个文件
- [分享]iTOP-4412开发板QT 持 HDMI 显示
- javascript运动框架---1