POJ1142 暴力质因数分解+素数测试

#include <iostream>
#include "stdio.h"
#include "stdlib.h"
#include "string.h"
#include "math.h"
using namespace std;
#define N 50005
typedef   long long int ll;
int primelist[1000]={0};
int isprime(int n);
void buildprimelist()
{
    int pos=0;
    for(int i=2;i<=10000;i++)
        if(isprime(i))  primelist[pos++]=i;
}
int ans1(int x)
{
    int ans=0;
    while(x)
    {
        ans+=x%10;
        x/=10;
    }
    return ans;
}
int ans2(int x)
{
    int ans=0;
    int temp[100],pos=0;
    for(int i=0;x!=1;)
    {
        if(isprime(x)) {
            temp[pos++]=x;
            break;
        }
        if(x%primelist[i]==0) {
            temp[pos++]=primelist[i];
            x=x/primelist[i];
        }
        else i++;
    }
    for(int i=0;i<pos;i++)
    {
        while(temp[i])
        {
            ans+=temp[i]%10;
            temp[i]/=10;
        }
    }
    return ans;
}
ll qpow(int  a,int  b,int  r)//快速幂
{
    ll ans=1,buff=a;
    while(b)
    {
        if(b&1)ans=(ans*buff)%r;
        buff=(buff*buff)%r;
        b>>=1;
    }
    return ans;
}
int miller_rabbin(int n,int a) //miller_rabbin算法
{
    int r=0,s=n-1;
    if(n%a==0) return 0;
    while(!(s&1)){ //将其转化为二进制的形式即 n=s*2^r
        s>>=1;
        r++;
    }
    ll k=qpow(a,s,n);
    if(k==1) return 1;
    for(int i=0;i<r;i++,k=k*k%n)
        if(k==n-1) return 1;
    return 0;
}
int isprime(int n) // 判断是否为素数
{
    int temp[]={2,3,5,7};
    if(n==1) return 0;
    for(int i=0;i<4;i++)
    {
        if(n==temp[i]) return 1;
        else if(!miller_rabbin(n,temp[i]))
            return 0;
    }
    return 1;
}

 int main()
 {
    int m;
    buildprimelist();
    while(scanf("%d",&m)!=EOF&&m)
    {
        for(int i=m+1;;i++)
        {
            if(isprime(i)) continue;
            if(ans1(i)==ans2(i)){
                printf("%d\n",i);
                break;
            }
        }
    }
        return 0;
 }