hdu 1215 七夕节

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
struct node
{
int x;
int n;
}pp[42000];
int T,N,ans,tmp;
int P[42000];
int quick_pow(int x,int n)
{
if(n==0)
return 1;
int cnt=1;
while(n)
{
if(n&1)
cnt=cnt*x;
x=x*x;
n>>=1;
}
return cnt;
}
bool isp(int x)
{
if(x==2)
return true;
if(x%2==0)
return false;
double upd=sqrt(x*1.0);
for(int i=3;i<upd;i+=2)
if(x%i==0)
return false;
return true;
}
void init()
{
int nn=0;
for(int i=2;i<=500000;i++)
if(isp(i))
P[nn++]=i;
//cout<<nn<<"\n";
}
int main()
{
init();
scanf("%d",&T);
while(T--)
{
scanf("%d",&N);
int num=0;
int m=N;
for(int i=0;i<41663;i++)
{
if(P[i]>m)
break;
if(m%P[i]==0)
{
pp[num].x=P[i];
pp[num].n=0;
while(m%P[i]==0)
{
pp[num].n++;
m/=P[i];
}
num++;
}
}
ans=1;
for(int i=0;i<num;i++)
{
tmp=0;
for(int j=0;j<=pp[i].n;j++)
tmp+=quick_pow(pp[i].x,j);
ans*=tmp;
}
ans-=N;
printf("%d\n",ans);
}
}