poj 1142

题目：http://poj.org/problem?id=1142

如果一个数本身是素数，那么这个数不是Smith Numbers，如果一个数不是素数，那么10 ^ 4( n 最多 8 位)以内一定有它的质因子，打出10000以内的素数，分解求和，从 n + 1往后枚举

const int  N = 10001;
int prime[N];
bool vis[N];
int num,tnum;
void is_prime()
{
    int i,j;
    num = 0;
    for(i = 2; i <= N; i++)
    {
        if(!vis[i]) prime[num ++] = i;
        for(j = 2; j * i <= N; j++)
        {
            if(vis[j * i]) continue;
            vis[j * i] = true;
        }
    }
}
bool is_prim(int n)
{
    int i;
    for(i = 2; i * i <= n; i++)
    {
        if(n % i == 0) return false;
    }
    return true;
}
int calsum(int n)
{
    int tsum = 0;
    while(n)
    {
        tsum += (n % 10);
        n /= 10;
    }
    return tsum;
}
int cal(int n)
{
    int i;
    int tsum = 0;
    for(i = 0; i < num; i++)
    {
        if(n == 1) break;
        while(n % prime[i] == 0)
        {
            int temp = prime[i];
            n /= prime[i];
            tsum += calsum(temp);
        }
    }
    if(n != 1) tsum += calsum(n);
    return tsum;
}
int main()
{
    int i,j;
    is_prime();
    int n;
    while(scanf("%d",&n) != EOF)
    {
        if(!n) break;
        i = n + 1;
        while(1)
        {
            if(is_prim(i)) {i++; continue;}
            if(calsum(i) == cal(i)) break;
            i++;
        }
        cout<<i<<endl;
    }
    return 0;
}

题目：http://codeforces.com/problemset/problem/233/B

被上面误导了一下，我看标签写的是二分，数学，然后就写二分，可是写好了发现不对呀，根本就不是单调函数怎么可以分出结果呢，然后开始想枚举x，但是发现太大了，但是s(x) 最大是81，于是就枚举 s(x)，这样就形成一个 x ^ 2 +  b * x + (-n)  = 0 2次方程，然后根据求根公式 x = (-b +(-) sqrt(b * b - 4 * a * c)) / 2 * a;然后求出 x，然后再判断

typedef long long ll;
ll cal(ll x)
{
    ll temp = 0;
    while(x)
    {
        temp += (x % 10);
        x /= 10;
    }
    return temp;
}
int main()
{
    ll n;
    int i,flag;
    while(cin>>n)
    {
        flag = 0;
        for(i = 1; i <= 100; i++)
        {
            ll m = (sqrt((i * i + 4 * n) * 1.0) - i) / 2;  
            ll ans = cal(m);
            if(cal(m) == i && m * m + i * m == n)
            {
                cout<<m<<endl;flag = 1;break;
            }
        }
        if(!flag) printf("-1\n");
    }
    return 0;
}