传送门
题意:给你A,B,N三个数,问你在A,B之间(include),有多少个数与N互质。
A、B已经到了1e15了,暴力是绝对不可能的,这辈子都不可能暴力的。
可以举例说明
m=12,n=30.
第一步:求出n的质因子:2,3,5;
第二步:(1,m)中是n的因子的倍数当然就不互质了(2,4,6,8,10,12)->n/2 6个 (3,6,9,12)->n/3 4个, (5,10)->n/5 2个。
答案显然不是全加起来。里面明显出现了重复的,我们现在要处理的就是如何去掉那些重复的了!
第三步:这里就需要用到容斥原理了,公式就是:n/2+n/3+n/5-n/(2*3)-n/(2*5)-n/(3*5)+n/(2*3*5)。发现除数是奇数个的时候是加,偶数个的时候是减。
(https://blog.csdn.net/spark__s/article/details/54983271)
这道题还可以学到m个数字可以组成的所有情况:
int full;
full=(1<<m)-1;
for(int i=1;i<=full;i++)
{
for(int j=1;j<=m;j++)
if(i>>(j-1)&1)
//可以取到第j个内容
}
最后发现老是WA,de了好久的bug,是因为返回类型是long long的solve返回值写成了int导致有奇奇怪怪的乱码出现。
#include<bits/stdc++.h>
#define show(a) cout<<#a<<'='<<a<<endl
using namespace std;
long long divn[1000];
long long full;
long long num=0;
void decomp(int n)
{
num=0;
for(int i=2;i*i<=n;i++)
{
if(n%i==0)
{
divn[++num]=i;
//show(i);
while(n%i==0)
n/=i;
}
}
if(n!=1)
divn[++num]=n;
}
long long solve(long long x)
{
long long res=0;
for(long long i=1;i<=full;i++)
{
long long cnt=0;
long long pro=1;
for(int j=1;j<=num;j++)
{
if((i>>(j-1))&1)
{
pro*=divn[j];
cnt++;
}
}
if(cnt&1)
res+=x/pro;
else
res-=x/pro;
}
return x-res;
}
int main()
{
int casen;
scanf("%d",&casen);
int sum=casen;
while(casen--)
{
long long a,b;
int n;
scanf("%I64d%I64d%d",&a,&b,&n);
decomp(n);
full=(1<<num)-1;
printf("Case #%d: %I64d\n",sum-casen,solve(b)-solve(a-1));
}
return 0;
}