HDU-5663-Hillan and the girl(Mobius反演)
来源:互联网 发布:linux shell 定时任务 编辑:程序博客网 时间:2024/06/07 23:45
题目链接:http://bestcoder.hdu.edu.cn/contests/contest_chineseproblem.php?cid=687&pid=1003
输出描述
对于每组测试数据,仅一行,一个整数,即问题的答案。
输入样例
21 233333310 10
输出样例
033
Hint
第一组数据中,f(i,j)f\left(i,j\right)f(i,j)显然始终为0,因为iii始终为1,gcd(i,j)\gcd\left(i,j\right)gcd(i,j)始终为完全平方数(始终为1)。
说白了这题跟BZOJ2820一毛一样,只是2820枚举质数,这题枚举平方数。。。
把2820的拿过来。。。前面枚举质数换成平方数就好。。。
对于F(k)可以在求mu的时候顺便求出来,
分块处理加速计算,因为对于(n/i)和(m/i)在一定区间内的值是一定的,根据这点可以每次跳过这段的计算。
这里的分块有点类似求和sum=n/1+n/2+n/3+...+n/(n-1)+n/n.直接去算O(n),分块O(sqrt(n)).
//#include <bits/stdc++.h>#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#define INF 0x3f3f3f3f#define LL long long#define bug cout<<"bug\n"using namespace std;const int MAXN = 1e7+7;const int MAXM = 1e9+7;int mu[MAXN],prime[MAXN],num_prime;long long F[MAXN];bool check[MAXN];void get_Mobius(){ memset(check,0,sizeof(check)); memset(F,0,sizeof(F)); mu[1]=1; num_prime=0; for(int i=2; i<MAXN; ++i) { if(!check[i]) { prime[num_prime++]=i; mu[i]=-1; } for(int j=0; j<num_prime; ++j) { if(i*prime[j]>MAXN)break; check[i*prime[j]]=1; if(i%prime[j]==0) { mu[i*prime[j]]=0; break; } else mu[i*prime[j]]-=mu[i]; } } for(int i=1; i*i<MAXN; ++i) for(int j=i*i,k=1; j<MAXN;++k, j+=i*i) F[j]+=mu[k]; for(int i=1; i<MAXN; ++i) F[i]+=F[i-1];}int main(){ int T; get_Mobius(); scanf("%d",&T); while(T--) { int n,m; long long ans=0; scanf("%d%d",&n,&m); if(n>m)swap(n,m); for(int i=1,j; i<=n; i=j+1) { j=min(n/(n/i),m/(m/i)); ans+=1LL*(F[j]-F[i-1])*(n/i)*(m/i); } /* for(int i=1; i<=n; ++i) ans+=F[i]*(n/i)*(m/i);*/ printf("%I64d\n",(LL)n*m-ans); } return 0;}
对于F(k)可以在求mu的时候顺便求出来,
分块处理加速计算,因为对于(n/i)和(m/i)在一定区间内的值是一定的,根据这点可以每次跳过这段的计算。
这里的分块有点类似求和sum=n/1+n/2+n/3+...+n/(n-1)+n/n.直接去算O(n),分块O(sqrt(n)).
0 0
- HDU-5663-Hillan and the girl(Mobius反演)
- HDU5663 Hillan and the girl[莫比乌斯反演]
- BestCoder Round #79 (div.2) 1001 1004 HDU 5663 Hillan and the girl
- 解题报告:HDU_5663 Hillan and the girl 莫比乌斯反演
- HDU-1695-GCD(Mobius反演/容斥)
- hdu 6134 Battlestation Operational mobius反演
- SPOJ PGCD (mobius反演 + 分块)
- NJUST-1923-triple(Mobius反演)
- BZOJ 2005 能量采集 (Mobius反演)
- BZOJ2005:能量采集(Mobius反演)
- hdu5663Hillan and the girl
- 【hdu 1695 】GCD 【容斥+欧拉 or mobius反演】
- hdu1695gcd(Mobius反演)
- Mobius反演总结
- BZOJ2154,BZOJ2693 (Mobius反演)
- BZOJ-2818-Gcd(欧拉函数/Mobius反演)
- BZOJ-2820-YY的GCD-(Mobius反演)
- HYSZB PROBLEM B(mobius反演+分块+容斥原理)
- 匈牙利算法
- poj1459——Power Network(网络流,超级源点、汇点)
- LTE物理传输资源(3)-时频资源
- 二分图的最大匹配、完美匹配和匈牙利算法
- EventBus3.0的入门和使用
- HDU-5663-Hillan and the girl(Mobius反演)
- Android提醒微技巧,你真的了解Dialog、Toast和Snackbar吗?
- C++堆、栈、自由存储区、全局/静态存储区和常量存储区收藏(转)
- 排序算法java三--堆排序
- JS之运动框架
- 剑指offer 面试题34:寻找丑数(Leetcode 263.ugly number)解题报告
- C++ Primer 学习笔记_82_模板与泛型编程 -类模板成员[续二]
- hdu 3810 Magina 搜索+队列模拟分组背包
- spring mvc整合shiro登录 权限验证