Hdu1695 (GCD) 终于写完了~

View Code  1 #include <stdio.h>
 2 #include <string.h>
 3 #define N 100010
 4 __int64 ol[N];
 5 int zys[N][20],g[N];
 6 void oula()
 7 {
 8 int i,j;
 9 ol[1]=1;
10 for (i=2;i<N;i++)
11 {
12 if (!ol[i])
13 {
14 for (j=i;j<N;j+=i)
15 {
16 if (!ol[j])
17 ol[j]=j;
18 ol[j]=ol[j]*(i-1)/i;
19 zys[j][g[j]++]=i;
20 }
21 }
22 ol[i]+=ol[i-1];
23 }
24 }
25 __int64 js(int x,int max,int n)
26 {
27 int i;
28 __int64 s=0;
29 for (i=x;i<g[n];i++)
30 s+=max/zys[n][i]-js(i+1,max/zys[n][i],n);
31 return s;
32 }
33 int main()
34 {
35 int t,f=0,i,a,b,k,max,min;
36 __int64 s;
37 memset(ol,0,sizeof(ol));
38 memset(g,0,sizeof(g));
39 memset(zys,0,sizeof(zys));
40 oula();
41 scanf("%d",&t);
42 while (t--)
43 {
44 scanf("%d%d%d%d%d",&i,&a,&i,&b,&k);
45 if (k==0)
46 {
47 printf("Case %d: 0\n",++f);
48 continue;
49 }
50 max=(a>b)?a:b;
51 min=a+b-max;
52 max/=k;
53 min/=k;
54 s=ol[min];
55 for (i=min+1;i<=max;i++)
56 s+=min-js(0,min,i);
57 printf("Case %d: %I64d\n",++f,s);
58 }
59 return 0;
60 }
61