HDOJ 4135 Co-prime && NOJ 1922 count_prime

HDOJ 4135 Co-prime 题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4135

NOJ 1922 count_prime 题目链接：https://icpc.njust.edu.cn/Problem/Local/1922/

NOJ 1922 count_prime是这次南京理工大学第八届程序设计大赛中的C题，这道题在HDOJ上做过一道基本上一模一样的题也就是4135的Co-prime，所以当时虽然还是贡献了几次罚时，但还是很快做出来了，这两道题的题意就是让我们求出在一个区间中某一个数互质的数的个数，由于是一个范围很多人拿到就会不知所措，但是其实我们可以分成两块来解答，区间l~r可以拆成1~r - 1~l，而如果是1~r和1~l的范围就要好求的多，一个质数x在1~r中的倍数个数就是r/x，同理在1~l中的倍数个数就是l/x，那么在l~r中的个数就是r/x - l/x。所以我们首先先求出所给数的质因数然后再运用这个方法，但是还是有一个问题，那就是会有重叠，也就是说可能一个数会既是某个质数的倍数又是另一个质数的倍数，所以我们得用容斥原理来解决这个问题，容斥原理是用到集合的知识，比如AVB = A+B-A^B。这就是最基础的容斥原理，运用这个原理我们就可以求得在l~r中间的个数了。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <string>
#include <cstring>
using namespace std;
typedef long long ll;
ll prime[70];
ll m,n;
void get_prime()
{
    for(ll i=2;i*i<=n;i++)
        if(n&&n%i==0) //对n进行素因数分解
        {
            prime[m++]=i;
            while(n&&n%i==0) n/=i;
        }
    if(n>1) prime[m++]=n;
}
ll solve(ll num)
{
    ll i,j;
    ll ans=0,tem,flag;
    for(i=1;i<1<<m;i++)
    {
        tem=1,flag=0;
        for(j=0;j<m;j++)
            if(i&1<<j)
                flag++,tem*=prime[j];
        if(flag&1) ans += num/tem; //容斥原理，奇加偶减
        else ans -= num/tem;
    }
    return ans;
}
int main()
{
	int t;
    ll a,b,k=1;
    scanf("%d",&t);
    while(t--)
    {
        m=0;
        scanf("%lld%lld%lld",&a,&b,&n);
        get_prime();
        printf("Case #%d: %lld\n",k++,b-solve(b)-(a-1-solve(a-1)));
    }
    return 0;
}