Another Easy Problem fzu1753

Another Easy Problem

Time Limit: 1000MS

Memory Limit: 32768KB

64bit IO Format: %I64d & %I64u

[Submit]   [Go Back]   [Status]

Description

xtt最近学习了高斯消元法解方程组,现在他的问题来了，如果是以下的方程，那么应该如何解呢?

C(n1,m1)==0 (mod M)

C(n2,m2)==0 (mod M)

C(n3,m3)==0 (mod M)

................

C(nk,mk)==0 (mod M)

xtt希望你告诉他满足条件的最大的M

其中C(i,j)表示组合数,例如C(5,2)=10,C(4,2)=6...

Input

输入数据包括多组，每组数据的第一行是一个正整数T(1<=T<=150)表示接下来描述的T个方程

接下来T行，每行包括2个正整数ni,mi (1<=mi<=ni<=100000)

Output

输出一行答案，表示满足方程组的最大M。

Sample Input

3
100 1
50 1
60 1

Sample Output

10

 #include<stdio.h>
 #include<string.h>
 #include<stdlib.h>
 #include<time.h>
 #include<iostream>
 #include<algorithm>
 using namespace std;
 int p[]= {},t,pp[];
 void init()
 {
     long long i,j;
     t=;
     for(i=; i<; i++)
     {
         if(!p[i])
         {
             p[t++]=i;
             j=i*i;
             while(j<)
             {
                 p[j]=;
                 j+=i;
             }
         }
     }
 }
 void fun(int n,int k)
 {
     int i,j,m,sum;
     for(i=; i<t; i++)
     {
         m=n;
         sum=;
         while(m)
         {
             m/=p[i];
             sum+=m;
         }

         m=k;
         while(m)
         {
             m/=p[i];
             sum-=m;
         }

         m=n-k;
         while(m)
         {
             m/=p[i];
             sum-=m;
         }
         if(pp[i]!=-)
         pp[i]=pp[i]>sum?sum:pp[i];
         else pp[i]=sum;
     }
 }
 int pow3(int a,int b)//求a^b
 {
     int r = ,base = a;
     while(b!=  )
     {
         if(b&)
             r*= base;
         base*= base;
         b>>= ;
     }
     return r;
 }

 int main()
 {
     init();
     int tt,m,n,i;
     while(~scanf("%d",&tt))
     {
         memset(pp,-,sizeof(pp));
         while(tt--)
         {
             scanf("%d%d",&m,&n);
             fun(m,n);
         }
         long long sum=;
         for(i=; i<t; i++)
         {
             if(pp[i])
                 sum*=pow3(p[i],pp[i]);
         }
         printf("%I64d\n",sum);
     }
 }