F - Co-prime HDU - 4135 （容斥原理）

题目

题意，求 【l ，r】中有多少和n互质的数。

首先暴力会超时（数据范围是10^15)

另外，互质也就是除了1没有别的公共因子，也就是说n的因子或者因子的倍数和n都不互质

这样就转化为了求n的因子的问题了，当然要先找到素因数，继而找所有素因数的倍数

n分解之后是多个素数相乘（忽略1）的形式，这里我用的埃式筛，当然，选择任何有效的算法都可行。

素因数的倍数和n也是不互质，所以要把所有素因数的倍数都排除掉，这里要用到容斥原理了。

例如12 素因数为2 3 那么2 的倍数有6个，3的倍数有4个，但是6的倍数被重复计算了两次，所以要减去一个6的倍数

也就是 sum=n/2+n/3-n/6.

区间的话， sum[r]-sum[l-1]就可以了。

代码

#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;
#include<vector>
#define ll long long
ll a[10000]; //素因数
int s;
int haha(ll n)  //筛选出素因数
{
    s=0;
    for(ll i=2; i*i<=n; i++)
    {
        if(n%i==0)
        {
            while(n%i==0)
                n/=i;
            a[s++]=i;
        }
    }
    if(n>1)
        a[s++]=n;
}
ll nop(ll m)  //容斥
{
    ll que[1000];  //素因数及素因数的倍数
    int top=0;
    que[top++]=-1;
    for(int i=0; i<s; i++)
    {
        int t=top;
        for(int j=0; j<t; j++)
        {
            que[top++]=que[j]*a[i]*(-1);
//            printf("%lld \n",que[top-1]);
        }
    }
    ll sum=0;
    for(int i=1; i<top; i++)
    {
        sum+=m/que[i];
    }
    return sum;
}
int main()
{
    ll n,m1,m2,t,op=0;
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld %lld %lld",&m1,&m2,&n);
        m1-=1;
        haha(n);
        ll sum=0,ans=0;
        sum=m1-nop(m1);
        ans=m2-nop(m2);
        op++;
        printf("Case #%lld: %lld\n",op,ans-sum);
    }
}