华东交通大学2016年ACM“双基”程序设计竞赛 1010

Problem Description

LB是个十分喜欢钻研的人，对什么事都要搞明白。有一天他学习到了阶乘，他十分喜欢，所以他在想一个问题。如果给定一个数n，求n！能不能被2016整除。LB算了好久都没有算出来，所以他向你求助，你能不能帮他解决这个问题呢？

Input

第一行只包含一个整数T（T≤1000），表示测试用例的个数。
对于每个测试用例，第一行只包含一个整数N（0≤N≤10000），N如描述所示。

Output

对于每个测试用例，输出一行，如果N！能被2016整除输出"YES"（不带引号），如果不能输出"NO"（不带引号）。

Sample Input

1
2016

Sample Output

YES

Author

zhengjinke2123

解法：此处扩展 http://blog.csdn.net/xieshimao/article/details/6342322

当然这里并没有那么麻烦

?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120

#include<stdio.h> #include<string.h>   struct node {     long p;     long n;     node()     {         n=0;     } }; long prime[50000]; bool flag[50000]; node stc[1000]; long num1=0,num2;   void init() {     unsigned long i,j;     memset(flag,0,sizeof(flag));     for(i=2; i<50000; i++)     {         if(!flag[i])         {             prime[num1++]=i;             for(j=i*i; j<50000; j=j+i)                 flag[j]=true;         }     } }   void devide(long m) {     long i;     num2=0;     for(i=0; prime[i]*prime[i]<=m; i++)     {         if(m%prime[i]==0)         {             stc[num2].p=prime[i];             stc[num2].n++;             m=m/prime[i];             while(m%prime[i]==0)             {                 stc[num2].n++;                 m=m/prime[i];             }             num2++;         }         if(m==1)             break;     }     if(m>1)     {         stc[num2].p=m;         stc[num2].n++;         num2++;     } }   bool judge(long n,long m) {     long i,sum,temp;     for(i=num2-1; i>=0; i--)     {         sum=0;         temp=stc[i].p;         while(n>=temp)         {             sum=sum+n/temp;             temp=temp*stc[i].p;         }         if(sum<stc[i].n)             return false;     }     return true; }   int main() {     long n,m,i;     init();     int t;     while(~scanf("%d",&t))     {         while(t--)         {             scanf("%d",&n);             m=2016;             for(i=0; i<1000; i++)                 stc[i].n=0;             if(n<=1&&m>1)             {                 puts("NO");                 continue;             }             if(m==0)             {                 puts("NO");                 continue;             }             if(n>=m)             {                 puts("YES");                 continue;             }             devide(m);             if(judge(n,m))             {                 puts("YES");             }             else             {                 puts("NO");             }         }     }     return 0; }