hdu 1159

唉。。。这两天真是多灾多难啊。。。净碰到蛋疼的题。犯一些蛋疼的错误（这题是在一个判断的地方关系弄错了）。

这题的主要原理是多个集合的容斥原理。

代码如下：

#include"stdio.h"

#include"math.h"



__int64 prime[10005]={2,3},count=2;

__int64 cnt,factor[20];

__int64 n,m,ans;



void cal()

{

    int i,j,flag;

    for(i=5;i<100000;i+=2)

    {

        flag=0;

        for(j=0;j<count;j++)

            if(i%prime[j]==0)

            {flag=1;break;}

        if(flag==0)

            prime[count++]=i;

    }

}



void cal1(int n)

{

    int i,lim=(int)sqrt(n)+1;

    cnt=0;

    for(i=0;prime[i]<=lim&&n!=1;i++)

    {

        if(n%prime[i]==0)

            factor[cnt++]=prime[i];

        while(n%prime[i]==0)

            n/=prime[i];

    }

    if(n!=1)

        factor[cnt++]=n;

}



void rescure(int x,int sign,int num,int dcn,int d)

{

    int i;

    if(sign==1)

    {

        num*=factor[d];

        dcn++;

    }

    if(d+1==cnt)

    {

        if(num==1)

            return ;

        if(dcn%2==1)

            ans+=(n-x)/num;

        else

            ans-=(n-x)/num;

        return ;

    }

    for(i=0;i<=1;i++)

        rescure(x,i,num,dcn,d+1);

}



int main( )

{

    __int64 a,b,c,d,k,i,j,t,ncase=0;

    __int64 sum;

    cal();

    scanf("%I64d",&t);

    while(t--)

    {

        ncase++;

        scanf("%I64d%I64d%I64d%I64d%I64d",&a,&b,&c,&d,&k);

        printf("Case %d: ",ncase);

        if(k>b||k>d||k==0)

        {

            printf("0\n");

            continue;

        }

        m=b<d?b/k:d/k;

        n=b>d?b/k:d/k;

        sum=n*m-m*(m+1)/2+1;

        ans=0;

        for(i=2;i<=m;i++)

        {

            cal1(i);

            for(j=0;j<=1;j++)

                rescure(i,j,1,0,0);

        }

        printf("%I64d\n",sum-ans);

    }

    return 0;

}