hdu 4135 co-prime

裸的容斥原理题

学习用

代码：

#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);
    }
 }