首页 > 代码库 > 华东交通大学2016年ACM“双基”程序设计竞赛 1010

华东交通大学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
当然这里并没有那么麻烦
#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;
}

  

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