POJ1365——Prime Land(素数筛选加暴力)

题目链接

     这是一道素数的水题，唯一的难点就是理解题意，大致是说给你一列数，a1,a2,a3,a4,a5,a6,.......然后按(a1^a2)*(a3^a4)*(a5^a6).......得到一个和值sum，然后把sum-1，再按(a1^a2)*(a3^a4)*(a5^a6).......的格式对sum-1分解质因式，比如样例中的5 1 2 1，就得到sum=5^1+2^1=10,然后sum-1得9，然后对9分解质因式得到3*3，也就是3^2，这也就是输出的结果为什么是3 2了。还有一点就是输出按质因数从大到小的顺序输出。我的做法是先输入数据并得到和值sum，再对从2到sum-1的范围内筛出质数，再依次从大到小去除sum-1，并不断更新sum，最后直到除到无法再除为止，由于是按从大到小来除的，最后输出也无需排序，直接把保存结果的数组依次输出即可。(ps:这段代码用C++交，G++交会超时)。

  

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

using namespace std;

int vis[100000];
int flag[100000];

int main()
{
    long long sum;
    int a, b;
    while(true)
    {
        sum=1;
        while(true)
        {
            scanf("%d",&a);
            char c=getchar();
            if(c=='\n')break;
            scanf("%d",&b);
            c=getchar();
            for(int i=1;i<=b;i++)sum*=a;
            if(c=='\n')break;
        }
        if(a==0)break;
        sum--;
        //cout<<sum<<endl;
        long long  m=sqrt(sum+0.5);
        memset(vis,0,sizeof(vis));
        for(long long i=2;i<=m;i++)if(!vis[i])
            for(long long j=i*i;j<=sum;j+=i)vis[j]=1;
        /*for(int i=1;i<=sum;i++)
            if(!vis[i])cout<<i<<endl;*/
        long long j1=0,j2=1;
        bool sign=false;
        memset(flag,0,sizeof(flag));
        for(long long i=sum;i>=2;i--)
        {
            sign=false;
            while(!vis[i])
            {
                if(sum%i==0)flag[j1]=i;
                else break;
                sum/=i;
                flag[j2]++;
                sign=true;
            }
            if(sign)
            {
                j1+=2;
                j2+=2;
            }
        }
        for(long long i=0;i<j2-2;i++)
            printf("%d ",flag[i]);
        printf("%d\n",flag[j2-2]);
            /*cout<<flag[i]<<' ';
        cout<<flag[j2-2]<<endl;*/
    }
    return 0;
}