D - GCD HDU - 1695 -模板-莫比烏斯容斥
阿新 • • 發佈:2019-01-13
- 思路: 都 除以 k 後轉化為 1-b/k 1-d/k中找互質的對數,但是需要去重一下 (x,y) (y,x) 這種情況。
- 這種情況出現 x ,y 肯定 都在 min (b/k, d/k) ,所以 奇數 最後 減去 一半 即可。
-
#include<bits/stdc++.h> using namespace std; #define ll long long #define maxn 1234567 bool vis[maxn+10]; ll t,n,m,prime[maxn+10]; ll mu[maxn+10],ans,a,b,c,d,k,sum; void getphi() { int cnt=0; mu[1]=1; for(int i=2; i<maxn; i++) { if(!vis[i]) { prime[++cnt]=i; mu[i]=-1; } for(int j=1; j<=cnt&&i*prime[j]<maxn; j++) { vis[i*prime[j]]=1; if(i%prime[j]==0) { mu[i*prime[j]]=0; break; } else mu[i*prime[j]]=-mu[i]; } } } int main() { getphi(); scanf("%lld",&t); for(int q=1; q<=t; q++) { ans=sum=0; scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&k); if(k==0) { printf("Case %d: 0\n",q); continue; } b/=k,d/=k; n=min(b,d); for(int i=1; i<=n; i++) { ans+=(ll)mu[i]*(b/i)*(d/i); sum+=(ll)mu[i]*(n/i)*(n/i); } printf("Case %d: %lld\n",q,ans-sum/2); } return 0; }