UVa 10042 && POJ 1142 - Smith Numbers

做了个非搜索题，smith数，又是由于忽略条件导致一直卡题，题目中说smith数是非素数，所以要单独判断，smith数比较密集，暴力枚举即可

此处只给出poj的代码，两者输入略有差别，uva的就不再给出

#include <iostream>
#include <cmath>
using namespace std;
bool isprime[10005];
int prime[1500];
int factor[1500];
void init()                   //打素数表
{
    for(int i=2; i<=10000; i++)
    {
        for(int j=i*2; j<=10000; j+=i)
            isprime[j]=1;
    }
    int p=0;
    for(int i=2; i<=10000; i++)
    {
        if(!isprime[i])
            prime[p++]=i;
    }
}
int main()
{
    init();
    int n,s1,s2,ans;
    while(cin>>n&&n)
    {
        ans=0;
        for(int i=n+1;; i++)
        {
            s1=0;
            s2=0;
            int t=i,sp=0;
            int e=i;
            while(e>0)       //s1是i的所有数字之和
            {
                s1+=e%10;
                e/=10;
            }
            for(int j=0; j<1229&&prime[j]<=sqrt(i)&&t>1; j++)
            {
                while(!(t%prime[j]))
                {
                    factor[sp++]=prime[j];
                    t/=prime[j];
                }
            }
            if(sp==0)continue;       //此处判断i是否为素数
            else
            {
                if(t>1)factor[sp++]=t;
                for(int j=0; j<sp; j++)      //s2是所有素因子数字之和
                {
                    int f=factor[j];
                    while(f>0)
                    {
                        s2+=f%10;
                        f/=10;
                    }
                }
                if(s1==s2)
                {
                    ans=i;
                    break;
                }
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}