欧拉函数
给定N,M,求1<=X<=N 且gcd(X,N)>=M条件下X的个数
#include<bits/stdc++.h>
using namespace std;
int euler(int n)
{
int ans=n;
for(int i=2;i*i<=n;i++)
{
if(n%i==0)
{
ans-=ans/i;
while(n%i==0) n/=i;
}
}
if(n>1) ans-=ans/n;
return ans;
}
int main()
{
int t,n,m;
scanf("%d",&t);
while(t--)
{
int ans=0;
scanf("%d%d",&n,&m);
for(int i=1;i*i<=n;i++)
{
if(n%i!=0) continue;
if(i>=m) ans+=euler(n/i);//i*(n/i)==n
if(n/i>=m&&i*i!=n) ans+=euler(i);
}
printf("%d\n",ans);
}
return 0;
}
/*
1<=x<=n,gcd(x,n)>=m
x=s*a<=n(a<=n/s),n=s*b(n%s==0); gcd(a,b)==1(a=prim); m<=s<=n;
*/
线性筛
给定N,1<=x,y<=N, 满足gcd(x,y)==prim 条件下可选的 ( x,y ) 对数。
#include<bits/stdc++.h>
using namespace std;
const int MAX=1e7+5;
bool mark[MAX];
int phi[MAX],prim[MAX],tot;
long long sum[MAX];
void phi_prim(int n)
{
memset(mark,0,sizeof(mark));
phi[1]=1;tot=0;
for(int i=2;i<=n;i++)
{
if(!mark[i]){prim[++tot]=i;phi[i]=i-1;}
for(int j=1;j<=tot;j++)
{
int x=prim[j];
if(i*x>n) break;
mark[i*x]=1;
if(i%x==0){phi[i*x]=phi[i]*x;break;}
else phi[i*x]=phi[i]*phi[x];
}
}
}
int main()
{
int n;
long long ans=0;
scanf("%d",&n);
phi_prim(n);
for(int i=1;i<=n;i++) sum[i]=sum[i-1]+phi[i];
for(int i=1;i<=tot;i++) ans+=sum[n/prim[i]]*2-1;
printf("%lld\n",ans);
return 0;
}
/*
1<=x,y<=n; gcd(x,y)=prim
x=s*a,y=s*b; s=prim, gcd(a,b)==1
ans+=phi[a]*2-1;
*/
杜教筛
const int MAX=1e7+5;
bool vis[MAX];
int prim[MAX],mu[MAX],fac[MAX],tot,pcnt;
map<long long,long long>dp;
void moblus()
{
mu[1]=1;tot=0;
for(int i=2;i<MAX;i++)
{
if(!vis[i]) {prim[tot++]=i;mu[i]=-1;}
for(int j=0;j<tot&&i*prim[j]<MAX;j++)
{
vis[i*prim[j]]=1;
if(i%prim[j]) mu[i*prim[j]]=-mu[i];
else {mu[i*prim[j]]=0;break;}
}
}
for(int i=2;i<MAX;++i) mu[i]+=mu[i-1];
}
long long M(long long x)
{
if(x<MAX) return mu[x];
if(dp[x]) return dp[x];
long long sum=1;//phi->x*(x+1)/2;
for(long long l=2,r;l<=x;l=r+1)
{
r=x/(x/l);
sum-=M(x/l)*(r-l+1);
}
return dp[x]=sum;
}
2018徐州ICPC网络赛D Easy Math 莫比乌斯反演/杜教筛
莫比乌斯反演
莫比乌斯反演的本质是容斥
反演分为两种形式:
①约数型
对于等式
可以得到
②倍数型
对于等式
可以得到
#include<bits/stdc++.h>
using namespace std;
const int MAX=1e5+5;
bool vis[MAX];
int prim[MAX],mu[MAX],tot;
void moblus()
{
memset(mu,0,sizeof(mu));
memset(prim,0,sizeof(prim));
memset(vis,0,sizeof(vis));
mu[1]=1;tot=0;
for(int i=2;i<MAX;i++)
{
if(!vis[i]) {prim[tot++]=i;mu[i]=-1;}
for(int j=0;j<tot&&i*prim[j]<MAX;j++)
{
vis[i*prim[j]]=1;
if(i%prim[j]) mu[i*prim[j]]=-mu[i];
else {mu[i*prim[j]]=0;break;}
}
}
}
int main()
{
int t,a,b,c,d,K;
moblus();
scanf("%d",&t);
for(int tc=1;tc<=t;tc++)
{
long long ans1=0,ans2=0;
scanf("%d%d%d%d%d",&a,&b,&c,&d,&K);
if(K==0) {printf("Case %d: 0\n",tc);continue;}
b/=K,d/=K;
int n=min(b,d);
for(int i=1;i<=n;i++)
ans1+=(long long)mu[i]*(b/i)*(d/i);
for(int i=1;i<=n;i++)
ans2+=(long long)mu[i]*(n/i)*(n/i);
printf("Case %d: %I64d\n",tc,ans1-ans2/2);
}
return 0;
}
/*
1<=x<=b ,1<=y<=d ,gcd(x,y)==k
1<=x/k<=b/k,1<=y/k<=d/k,gcd(x/k,y/k)==1
a=x/k ,b=y/k ,gcd(a,b)==1
a与b互质
*/
翻译:vici@cust