做了个非搜索题,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;
}