POJ-1142 & HDOJ-1333 Smith Numbers 解题报告

       分解质因数题。题意：所谓史密斯数就是说这个数分解成质因数后每一个质因数各个数位上的数字之和相加的结果与这个数各个数位上的数字之和相等，并且这个数不能是质数。现在给你一个不超过8位的数字n，求出大于n的最小史密斯数。

       我的解题思路：分解质因数后判断就行了，注意素数不是史密斯数，另外因为给的是不超过8位的数字，所以当给的数字是99999999的时候答案是会超过八位的，因此计算数位的时候至少要计算9位。在POJ上有99999999这样的极限数据。

       我的解题代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <cmath>
#include <algorithm>

using namespace std;

const int N = 10000;

bool isprime[N];
int primes[N], pn;
int prmft[N], power[N], pfn;    //分解的质因数，质因数的个数，不同质因数的个数
int num;

void InitRead();

void DataProcess();

void FastSieve(int maxn);

void Factor(int now);

bool IsSmith(int now);

int GetDigit(int now);  //获得各个数位之和

int main()
{
    //freopen("/home/jzqt/Code/in.txt", "r", stdin);
    //freopen("/home/jzqt/Code/out2.txt", "w", stdout);
    InitRead();
    while (~scanf("%d", &num))
    {
        if (num == 0) break;
        DataProcess();
    }
    return 0;
}

void InitRead()
{
    memset(isprime, true, sizeof(isprime));
    isprime[0] = isprime[1] = false;
    pn = 0;
    FastSieve(N - 1);
    return;
}

void DataProcess()
{
    for (int i=num+1; ; ++i)
    {
        Factor(i);
        if (IsSmith(i))
        {
            printf("%d\n", i);
            break;
        }
    }
    return;
}

void FastSieve(int maxn)
{
    for (int i=2; i<=maxn; ++i)
    {
        if (isprime[i]) primes[pn++] = i;
        for (int j=0; j<pn; ++j)
        {
            if (i * primes[j] > maxn) break;
            isprime[i * primes[j]] = false;
            if (i % primes[j] == 0) break;
        }
    }
    return;
}

void Factor(int now)
{
    pfn = 0;
    memset(power, 0, sizeof(power));
    for (int i=0; i<pn; ++i)
    {
        if (now % primes[i] == 0)
        {
            while (now % primes[i] == 0)
            {
                now /= primes[i];
                power[pfn]++;
            }
            prmft[pfn++] = primes[i];
        }
    }
    if (now != 1)
    {
        power[pfn]++;
        prmft[pfn++] = now;
    }
    return;
}

bool IsSmith(int now)
{
    if (pfn == 1 && power[0] == 1) return false;
    int a, b = 0;
    a = GetDigit(now);
    for (int i=0; i<pfn; ++i)
    {
        b += GetDigit(prmft[i]) * power[i];
    }
    return a == b;
}

int GetDigit(int now)
{
    int ans = 0, temp = 1;
    for (int i=0; i<9; ++i)
    {
        ans += now / temp % 10;
        temp *= 10;
    }
    return ans;
}