好多的树

时间限制：3000 ms  |  内存限制：65535 KB
难度：5
 
描述

在那遥远的地方有一片神奇的森林，它的神奇之处是：森林中的每棵树都长在一对整数确定的坐标上。有一个小红人站在(0,0)的位置上，放眼望去，看到很多的树，眼前的树是那么的多，以至于它一直数不清。那么就来写个程序帮它数数吧！

 

 
输入
第一行一个整数n,代表测试数据组数。
接下来有n(n<=20)行数，每行数有两个整数i，j代表任何的（x，y）(0<x<=i<=100000,0<y<j<=100000)的整数坐标上都有一棵树。
小红人始终是站在（0，0) 点上看树的。
输出
输出小红人看到的树的个数。每个结果占一行。
样例输入
1
7 4
样例输出
20
View Code 
  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cmath>
  4 #define Max 100001
  5 using namespace std;
  6 int prime[10000];
  7 bool isprime[Max];
  8 int cont=0;
  9 void prm()
 10 {
 11     int i,j;
 12     for(i=0;i<Max;i++) isprime[i]=true;
 13 
 14     for(i=2;i<Max;i++)
 15     {
 16         if(isprime[i])
 17         {
 18             prime[cont++]=i;
 19             for(j=2*i;j<Max;j+=i) isprime[j]=false;
 20         }
 21     }
 22 }
 23 int main()
 24 {
 25     int t,T;
 26     int i,k,j,p,q,r,totle;
 27     long long n,m,ans;
 28     int index[10];
 29     prm();
 30     cin>>T;
 31     while(T--)
 32     {
 33         cin>>n>>m;
 34         //scanf("%lld%lld",&n,&m);
 35         ans=n*m;
 36         for(i=2;i<=n;i++)
 37         {
 38             if(isprime[i]) {ans-=m/i;continue;}
 39             k=0;
 40             for(j=0,t=i;prime[j]<=sqrt(i);j++)
 41             {
 42                 if(t%prime[j]==0)
 43                 {
 44                     index[k++]=prime[j];
 45 
 46                     while(t%prime[j]==0)  t/=prime[j];
 47                     if(t==1) break;
 48                     if(isprime[t]) {index[k++]=t;break;}
 49                 }
 50             }
 51             if(k==1) {ans-=m/index[0];continue;}
 52             if(k==2) {ans-=( m/index[0]+m/index[1]-m/(index[0]*index[1]) );continue;}
 53             if(k==3)
 54             {
 55                 totle=1;
 56                 for(p=0;p<k;p++) ans-=m/index[p],totle*=index[p];///1
 57                 for(p=0;p<k;p++) ans+=m/(totle/index[p]);///2
 58                 ans-=m/totle;////3
 59                 continue;
 60             }
 61             if(k==4)
 62             {
 63                 totle=1;
 64                 for(p=0;p<k;p++) {ans-=m/index[p];totle*=index[p];}////1
 65 
 66                 for(p=0;p<k;p++)///2ºÍ3
 67                 {
 68                     for(q=p+1;q<k;q++)
 69                     {
 70                         ans+=m/( index[p]*index[q] );///2
 71                     }
 72                 }
 73                 for(p=0;p<k;p++) ans-=m/(totle/index[p]);////3
 74                 ans+=m/totle;///4
 75                 continue;
 76             }
 77             if(k==5)
 78             {
 79                 totle=1;
 80                 for(p=0;p<k;p++) ans-=m/index[p],totle*=index[p];////1
 81 
 82                 for(p=0;p<k;p++)///2///3
 83                 {
 84                     for(q=p+1;q<k;q++)
 85                     {
 86                         ans+=m/( index[p]*index[q] );////2
 87                         ans-=m/(totle/(index[p]*index[q]));/////3
 88                     }
 89                 }
 90                 for(p=0;p<k;p++) ans+=m/(totle/index[p]);///4
 91                 ans-=m/totle;////5
 92                 continue;
 93             }
 94             if(k==6)
 95             {
 96                 totle=1;
 97                 for(p=0;p<k;p++) ans-=m/index[p],totle*=index[p];///1
 98 
 99                 for(p=0;p<k;p++)//2
100                 {
101                     for(q=p+1;q<k;q++)
102                     {
103                         ans+=m/( index[p]*index[q] );////2
104                         ans+=m/(totle/(index[p]*index[q]));////4
105                     }
106                 }
107 
108                 for(p=0;p<k;p++)///3
109                 {
110                     for(q=p+1;q<k;q++)
111                         for(r=q+1;r<k;r++)
112                             ans-=m/( index[p]*index[q]*index[r] );
113                 }
114                 for(p=0;p<k;p++) ans-=m/(totle/index[p]);///5
115                 ans+=m/totle;///6
116                 continue;
117             }
118         }
119         //printf("%lld\n",ans);
120         cout<<ans<<endl;
121     }
122 
123     return 0;
124 }
posted @ 2012-03-06 15:45  qijinbiao1  阅读(343)  评论(0)    收藏  举报