hdu 1695 欧拉筛+容斥原理

题目传送门

//欧拉筛+容斥原理
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn=100000+5;
int phi[maxn],tot=0,prime[maxn];
LL sum[maxn],ans;
bool check[maxn],done[maxn];
int fac[maxn],cnt,x,y,sum1;
void dfs(int i,int nu,int f,int mu)
{
    if(nu==f) {
        sum1+=x/mu;
        return;
    }
    if(i==cnt) return;
    dfs(i+1,nu+1,f,mu*fac[i]);
    dfs(i+1,nu,f,mu);
}
int rongchi()                       //容斥原理
{
    int s=0;
    for(int i=1;i<=cnt;i++)
    {
        sum1=0;
        dfs(0,0,i,1);
        if(i&1) s+=sum1;
        else s-=sum1;
    }
    return x-s;
}
int main()
{
    //欧拉筛+前缀和
    memset(check,false,sizeof(check));
    phi[1]=1;
    for(int i=2;i<maxn-2;i++)
    {
        if(!check[i]) {
            phi[i]=i-1;
            prime[tot++]=i;
        }
        for(int j=0;j<tot;j++)
        {
            if(i*prime[j]>maxn-2) break;
            check[i*prime[j]]=true;
            if(i%prime[j]==0) {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else {
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
            }
        }
    }
    sum[1]=phi[1];
    for(int i=2;i<=maxn-2;i++)
        sum[i]=sum[i-1]+phi[i];
    //cout<<sum[3]<<endl;
    int T,kase=0;
    scanf("%d",&T);
    while(T--)
    {
        int a,b,c,d,k;
        ans=0;
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        if(k==0) {
            printf("Case %d: 0\n",++kase);
            continue;
        }
        x=b,y=d;
        if(x>y) swap(x,y);
        x/=k,y/=k;
        ans+=sum[x];
        int t=x+1;
        //容斥原理     在[1,x]中寻找后面的区间数的互质数的个数
        for(int i=t;i<=y;i++)
        {
            //整数分解
            int t=i;
            cnt=0;
            for(int j=0;j<tot&&prime[j]*prime[j]<=i;j++)
            {
                if(t%prime[j]==0) {
                    fac[cnt++]=prime[j];
                    while(t%prime[j]==0) {
                        t/=prime[j];
                    }
                }
                if(t==1) break;
            }
            if(t>1) {
                fac[cnt++]=t;
            }
            //枚举因子    奇加偶减
            ans+=(LL)rongchi();
        }
        printf("Case %d: %lld\n",++kase,ans);
    }
    return 0;
}