裸的容斥原理题
学习用
代码:
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <map>
#include <set>
#include <queue>
using namespace std;
int stacks[1005];
int len;
long long a,b;
long long numa=0;
long long numb=0;
int ks=1;
int dfs(int lens,int flag,int num)
{
if(lens==len)
{
if(flag)
{
numa-=(a/num);
numb-=(b/num);
}
else
{
numa+=(a/num);
numb+=(b/num);
}
return 1;
}
dfs(lens+1,flag,num);
dfs(lens+1,!flag,num*stacks[lens]);
}
int main()
{
int n;
scanf("%d",&n);
int c;
while(n--)
{
len=0;
scanf("%I64d%I64d%d",&a,&b,&c);
int k=c;
a--;
for(int i=2;i<=sqrt(c);i++)
{
if(k%i==0)
{
stacks[len++]=i;
while(k%i==0)k/=i;
}
}
if(k!=1)
stacks[len++]=k;
numa=0;
numb=0;
dfs(0,0,1);
// printf("%I64d %I64d\n",numa,numb);
printf("Case #%d: %I64d\n",ks++,numb-numa);
}
}