HDU：不可模数

不可摸数 

Time Limit : 2000/1000ms (Java/Other)   Memory Limit : 65536/32768K (Java/Other)

Total Submission(s) : 25   Accepted Submission(s) : 6

Problem Descriptions(n)是正整数n的真因子之和，即小于n且整除n的因子和.例如s(12)=1+2+3+4+6=16.如果任何

数m，s(m)都不等于n,则称n为不可摸数.

Input包含多组数据，首先输入T,表示有T组数据.每组数据1行给出n(2<=n<=1000)是整数。

Output如果n是不可摸数，输出yes，否则输出no 

Sample Input3

2

5

8 

Sample Outputyes

yes

no

解题思路：这道题要打表做，先将所有数的真因子之和制表。首先我们要确定的是我们要对那些数制表，这道题目打表的数打到五十

多万就可以。现在我来推导一个五十多万是怎么推出来的，一个很大的数，它的真因子和可能会很小，因为输入的n是2~1000之间的数，

假设存在这么一个数X，1~1000中的每一个都是它的因子（当然这只是假设，极限情况），那么S（x)=2+3+4...+n...+1000=500499

首先1是所有数的真因子：所有元素中先加上1

然后2是所有偶数的真因子：所有逻辑下标为2的倍数的位置的元素中再加上2

然后是3:在所有逻辑下标为3的位置的元素中加3

。。。。。。。

代码：

#include<string.h>

#include<stdio.h>

#define M 500003

int a[M];

int b[M];

void biao()

{

    int i,j;

    memset(a,0,sizeof(a));

    memset(b,0,sizeof(b));

    for(i=1;i<M/2;i++)

        for(j=2*i;j<M;j=j+i)

        {

            a[j]=a[j]+i;

        }

    for(i=2;i<M;i++)//每一个数中都存放了他的真因子和

        if(a[i]<=1000)//如果某个数的真因子和时小于或等于1000，说明存在真因子为a[i]的数字（不论几个）那么a[i]都不是不可模数

            b[a[i]]=1;//将a[i]作为b数组的下标，并附值为1，代表a[i]不是不可摸数

}

int main()

{

    int T,N;

    biao();//打表

    scanf("%d",&T);

    while(T--)

    {

        scanf("%d",&N);

        if(b[N]==1)

            printf("no\n");

        else

            printf("yes\n");

    }

    return 0;

}