hdu1215七夕节

[align=left]Problem Description[/align]
七夕节那天,月老来到数字王国,他在城门上贴了一张告示,并且和数字王国的人们说:"你们想知道你们的另一半是谁吗?那就按照告示上的方法去找吧!" 人们纷纷来到告示前,都想知道谁才是自己的另一半.告示如下:

数字N的因子就是所有比N小又能被N整除的所有正整数,如12的因子有1,2,3,4,6. 你想知道你的另一半吗?

[align=left]Input[/align]
输入数据的第一行是一个数字T(1<=T<=500000),它表明测试数据的组数.然后是T组测试数据,每组测试数据只有一个数字N(1<=N<=500000).

[align=left]Output[/align]
对于每组测试数据,请输出一个代表输入数据N的另一半的编号.

[align=left]Sample Input[/align]

3 2 10 20

[align=left]Sample Output[/align]

1 8 22

[align=left]Author[/align]
Ignatius.L

[align=left]Source[/align]
杭电ACM省赛集训队选拔赛之热身赛

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>

using namespace std;

int an[510000],num;
bool bn[510000];
void prime()
{
num = 0;
memset(bn,1,sizeof(bn));
bn[0] = bn[1] = 0;
int i,j;
for(i = 2;i<=500000;i++)
{
if(bn[i]) an[num++] = i;
for(j = 0;an[j]*i<=500000;j++)
{
if(i%an[j] == 0)break;
bn[i*an[j]] = 0;
}
}
}

int main()
{
int z;
prime();
cin>>z;
while(z--)
{
int a;
scanf("%d",&a);
int aa = a;
int i,j,k;
int t,sum,ans = 1;
for(i = 0;a!=1;i++)
{
if(a%an[i] == 0)
{
t = 1,sum = 1;
while(a%an[i] == 0)
{
a/=an[i];
t *= an[i];
sum+=t;
}
ans *= sum;
}
}
ans -= aa;
printf("%d\n",ans);
}
return 0;
}