1215 七夕节

七夕节

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

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

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

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

Sample Input
3
2
10
20

Sample Output
1
8
22

这道题用到了素数打表的方法，首先筛选出素数然后判断假如是素数输出1，否则将他所有因子找到
最后测试了一下自己的两个素数快速打表的模板
第一个：

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e6+5;
bool visit[101000000];
int prime[100000000];

void init_prim()
{
    memset(prime, 0, sizeof(prime));
    memset(visit, true, sizeof(visit));
    int num = 0;
    for (int i = 2; i <= maxn; ++i)
    {
        if (visit[i] == true)
        {
            num++;
            prime[num] = i;
        }
        for (int j = 1; ((j <= num) && (i * prime[j] <= maxn));  ++j)
        {
            visit[i * prime[j]] = false;
            if (i % prime[j] == 0) break; //µã¾¦Ö®±Ê
        }
    }
}
int main()
{
    int n,kase=0;
    init_prim();
    while(scanf("%d",&n)&&n) {
        printf("Case %d: %d\n",++kase,visit[n]);
    }
    return 0;
}

第二个：

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e6+5;
vector<int>prime(maxn,1),p(maxn,0);

void isprime()
{
    int cnt=0;
    for(int i=2;i<maxn;i++) {
        if(prime[i]) {
            p[cnt++]=i;
            for(int j=i*2;j<maxn;j+=i) {
                prime[j]=0;
            }
        }
    }
    return ;
}

int main()
{
    isprime();
    for(int i=0;i<1000;i++)
        cout<<i<<" "<<prime[i]<<endl;
    /*for(int i=0;i<1000;i++)
        cout<<i<<" "<<p[i]<<endl;
        for(int i=0;i<10000;i++)
        cout<<" "<<p[i];*/
    return 0;
}

这两个模板在大数方面各有各的好处
不过应该差不了多少，刚开始的时候一直T，后来加速了输入流，就过了
附上两个代码
第一个：

#include<bits/stdc++.h>
using namespace std;
using LL =int64_t;
const int maxn=5e5+5;
bool visit[1010000];
int prime[1000000];

void init_prim()
{
    memset(prime, 0, sizeof(prime));
    memset(visit, true, sizeof(visit));
    int num = 0;
    for (int i = 2; i <= maxn; ++i)
    {
        if (visit[i] == true)
        {
            num++;
            prime[num] = i;
        }
        for (int j = 1; ((j <= num) && (i * prime[j] <= maxn));  ++j)
        {
            visit[i * prime[j]] = false;
            if (i % prime[j] == 0) break; //µã¾¦Ö®±Ê
        }
    }
}
int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    LL n;
    cin>>n;
    init_prim();
    for(int i=0;i<n;i++) {
        LL temp;
        cin>>temp;
        if(visit[temp]) cout<<1<<endl;
        else {
            LL sum=1;
            for(int i=2;i*i<=temp;i++) {
                if(temp%i==0) {sum+=i;
                if(temp%(temp/i)==0&&temp/i!=i) sum+=temp/i;}
            }
            cout<<sum<<endl;
        }
    }
    return 0;
}

第二个：

#include<bits/stdc++.h>
using namespace std;
using LL =int64_t;
const int maxn=5e5+6;
vector<int>prime(maxn,1),p(maxn,0);
void isprime()
{
    int cnt=0;
    for(int i=2;i<maxn;i++) {
        if(prime[i]) {
            p[cnt++]=i;
            for(int j=i*2;j<maxn;j+=i) {
                prime[j]=0;
            }
        }
    }
    return ;
}

int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    LL n;
    cin>>n;
    isprime();
    for(int i=0;i<n;i++) {
        LL temp;
        cin>>temp;
        if(prime[temp]) cout<<1<<endl;
        else {
            LL sum=1;
            for(int i=2;i*i<=temp;i++) {
                if(temp%i==0) {sum+=i;
                if(temp%(temp/i)==0&&temp/i!=i) sum+=temp/i;}
            }
            cout<<sum<<endl;
        }
    }
    return 0;
}

第一个用了967 MS，第二个用了795 MS，果然第二个比较好用而且内存比较小