hdu 1695 欧拉筛+容斥原理

题目传送门

//欧拉筛+容斥原理
#include
#include
#include
#include
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;imaxn-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<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;j1) {
                fac[cnt++]=t;
            }
            //枚举因子    奇加偶减
            ans+=(LL)rongchi();
        }
        printf("Case %d: %lld\n",++kase,ans);
    }
    return 0;
}