hdu 1215 七夕节

参考博客：http://blog.csdn.net/niushuai666/article/details/7346408
收获：线性素数筛法、素数分解法
公式博客：http://www.cnblogs.com/baidongtan/archive/2012/08/30/2663015.html
还有关于素数筛法的三种比较：http://blog.csdn.net/sjf0115/article/details/8693756

#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<cmath>
#include<algorithm>
using namespace std;
#define MAX 500000
bool vis[MAX+5];
int prime[MAX+5];
void print_prime(){//线性求素数法
    memset(vis,false,sizeof(vis));
    int i,index=0;//注意初始化！！
    int half=sqrt(MAX*1.0)+1;
    for(i=2;i<half;i++){
        if(!vis[i]){
            prime[index++]=i;
        }
        for(int j=0;j<index&&prime[j]*i<=MAX;j++){
            vis[prime[j]*i]=true;
            if(i%prime[j]==0){
                break;
            }
        }
    }
}
int work(int n){//素数分解
    int i=0,total=1,temp=n;
    for(;prime[i]*prime[i]<=n;i++){
        int sum=1,num=1;
        while(n%prime[i]==0){
            num*=prime[i];
            n/=prime[i];
            sum+=num;
        }
        total*=sum;
    }
    if(n!=1){//可能剩下的落单的最大素数
        total*=n+1;
    }
    return total-temp;
}
int main()
{
    int t;
    cin>>t;
    print_prime();
    //cout<<prime[0]<<endl;
    while(t--){
        int n;
        cin>>n;
        cout<<work(n)<<endl;
    }
    return 0;
}