HDU 1695（GCD）

#include<iostream>
using namespace std;
#define max  100005
long  long uler[max];
int num[max];
int p[max][30];
void init()
{
    uler[1]=1;
    for(int i=2;i<max;i++)
    {
        
    if(!uler[i])
    {
        for(int j=i;j<max;j+=i)
        {
            if(!uler[j]) uler[j]=j;
            uler[j]=uler[j]*(i-1)/i;
            p[j][num[j]++]=i;
        }
        
    }
    uler[i]+=uler[i-1];
   }
}
int dfs(int k,int b,int now)
{
    int ans=0;
    for(int i=k;i<num[now];i++)
    ans+=b/p[now][i]-dfs(i+1,b/p[now][i],now);
    return ans;
}
int main()
{
    int cas,cas1=1;
    init();
    cin>>cas;
    while(cas--)
    {
        int a,b,c,d,k;
       scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        if(k==0) {cout<<"Case "<<cas1++<<": "<<0<<endl;continue;}
        if(b>d) swap(b,d);
        int m=b/k,n=d/k;
           long long ans=0;
        ans+=uler[m];       
        for(int i=m+1;i<=n;i++)
        ans+=m-dfs(0,m,i);
        cout<<"Case "<<cas1++<<": "<<ans<<endl;  
    }
    return 0;
}