两人刚结束电话交谈，史密斯突然灵感来临，对威兰斯基的电话号码“4937775”产生了兴趣，总觉得这是个特别的数。可它的特殊之处究竟在哪儿呢？史密斯开始思索考证起来，他先把 4937775 分解质因数：4937775=3×5×5×65837，然后再把 4937775 所有质因数各位上的数字相加得：3＋5＋5＋6＋5＋8＋3＋7=42，接着他又把 4937775 各位上的数字相加得：4＋9＋3＋7＋7＋7

＋5=42，秘密终于找到了，原来这两个和相等。这真有意思，难道是巧合么？有没有其他的数也有此特点呢？结果发现，所有质数都是具有如此特点。

后来的数学家们把这样的数叫做“史密斯数”，而且还决定质数（简单不复杂）不属于斯密斯数。除质数之外还有许多数具有这样独特的性质，其中最小的数是 4。大家不妨检查一下，4=2×2，2＋2=4。类似有，22=2×11,2＋2=2＋1＋1；27=3×3×3，2＋7=3＋3＋3。

你的任务是寻找最接近而且大于给定的数的斯密斯数。

smith.in

4937774

smith.out

4937775

#include<map>
#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#define FORa(i,s,e) for(int i=s;i<=e;i++)
#define FORs(i,s,e) for(int i=s;i>=e;i--)
#define gc pa==pb&&(pb=(pa=buf)+fread(buf,1,10000,stdin),pa==pb)?EOF:*pa++
#define File(name) freopen(name".in","r",stdin),freopen(name".out","w",stdout)using namespace std;
static char buf[10000],*pa=buf,*pb=buf;
inline int read();const int N=4000000,BN=278412;
int n,cnt,prime[N+1],v[N+1];
int a[BN+4]={4,22,27,58,85,94,121,166,202,}//后面的省略，不然上传不了啊！
map<int,bool> mp;
void Oulashai()
{FORa(i,2,n){if(!v[i]) v[i]=i,prime[++cnt]=i,mp[i]=1;FORa(j,1,cnt){if(prime[j]>v[i]||prime[j]*i>n) break;v[prime[j]*i]=prime[j];}}
}
int Back(int x)
{int cnt=0;while(x) cnt=cnt+x%10,x/=10;return cnt;
}
bool Check(int x)
{int fx=x,cnt1=0,cnt2=Back(x);FORa(i,2,sqrt(fx)) {int ct=0,backi=Back(i);if(fx%i==0){while(fx%i==0){fx/=i,ct++;    if(cnt1+ct*backi>cnt2) return 0;} cnt1=cnt1+ct*backi;        }}if(fx>1) cnt1+=Back(fx);return cnt1==cnt2;
}
int main()
{File("smith");n=read(),Oulashai();
/*    FORa(i,1000000,n){if(!mp[i])if(Check(i))printf("%d,",i);}*/FORa(i,0,BN+2) mp[a[i]]=1;FORa(i,n+1,N+20){if(mp[i]) {printf("%d",i);exit(0);    }}return 0;
}
inline int read()
{register int x(0);register int f(1);register char c(gc);while(c<'0'||c>'9') f=c=='-'?-1:1,c=gc;while(c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=gc;return x*f;
}

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#define FORa(i,s,e) for(int i=s;i<=e;i++)
#define FORs(i,s,e) for(int i=s;i>=e;i--)
#define gc pa==pb&&(pb=(pa=buf)+fread(buf,1,10000,stdin),pa==pb)?EOF:*pa++
#define File(name) freopen(name".in","r",stdin),freopen(name".out","w",stdout)using namespace std;
static char buf[10000],*pa=buf,*pb=buf;
inline int read();const int N1=100000000,N2=5762455;
int n,cnt,prime[N2+1];
bool bz[N1+1];
void Oulashai(int fn)
{FORa(i,2,fn){if(!bz[i]) prime[++cnt]=i;FORa(j,1,cnt){if(i*prime[j]>fn) break;         bz[i*prime[j]]=true;if(i%prime[j]==0) break;}}
}
int Back(int x)
{int ct=0;while(x) ct=ct+x%10,x/=10;return ct;
}
bool Check(int x)
{int fx=x,cnt1=0,cnt2=Back(x);FORa(i,2,sqrt(fx)) {int ct=0,backi=Back(i);if(fx%i==0){while(fx%i==0){fx/=i,ct++;    if(cnt1+ct*backi>cnt2) return 0;} cnt1=cnt1+ct*backi;        }}if(fx>1) cnt1+=Back(fx);return cnt1==cnt2;
}
int main()
{File("smith");n=read(),Oulashai(n+2000);int p=n;while(1){++p;if(bz[p]&&Check(p)) {printf("%d",p);exit(0);}}return 0;
}
inline int read()
{register int x(0);register int f(1);register char c(gc);while(c<'0'||c>'9') f=c=='-'?-1:1,c=gc;while(c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=gc;return x*f;
}