HDU4937题解

/*题目大意：给定一个数n，若这个数在base进制下全由3，4，5，6组成的话，
则称base为n的幸运进制，给定n，求有多少个幸运进制。无穷多个的话输出-1，
单个位置上超过9用相应的字符表示。
解题思路：首先n%base=3|4|5|6，假设我们当前考虑3，即n%base=3，那么就有://这一句注意！！！ 
n=base*k+3
n-3=base*k
将枚举n-3的因子作为base，判断剩下的k是否也满足条件。注意这里枚举的因子要大于3，因为n%base=3. 
将可行的base放入set中去重即可。(set是什么？miao)
注意这里的无穷多解，只有在n=3|4|5|6的时候，因为对应34的话是一个字符，而不是现实34.*/
#include <cstdio>
#include <cstring>
#include <ctime>
#include <cstdlib>
#include <set> 
#include <algorithm>

using namespace std;//去学set去重啊o_O 
typedef long long ll;//好吧，以后用typedef 
const int maxn = 1e6;//以后就用这种姿势吧 

ll N;
int np, vis[maxn+5], prime[maxn+5];
set<ll> ans;//这是什么 

int c, cnt[maxn+5];
ll nowbase, fact[maxn+5];

void prime_table (int n)//打素数表，而且还是那种我无法理解的打素数表，我还是用我原来那个吧 
 {
    np = 0;
    memset(vis, 0, sizeof(vis));

    for (int i = 2; i <= n; i++)
	 {
        if (vis[i])
            continue;
		prime[np++] = i;//记录2-n未被访问的数 
        for (int j = 2 * i; j <= n; j += i)
            vis[j] = 1;
     }
}

void div_factor (ll n)//对n分解质因数，fact存数，cnt存个数 
{
    c = 0;
    for (int i = 0; i < np && n >= prime[i]; i++) //未被访问的那些数 
	{

        if (n % prime[i] == 0)
		 {

            cnt[c] = 0;
            fact[c] = prime[i];//为这个未被访问的数 

            while (n % prime[i] == 0)
			 {
                n /= prime[i];//分解质因数么 
                cnt[c]++;// 有质因数fact[c]的个数 
             }
            c++;//再依次分解去 
         }
    }

    if (n != 1)//分解完之后素数np范围外的素数。 
	 {
        fact[c] = n;
        cnt[c++] = 1;//这就是分解质因啊 
     }
}

bool judge (ll n, ll base) //判断剩下的k是否也满足条件。 
{

    if (base <= nowbase)//进制不可以小于这一位的数 ,注意这里枚举的因子要大于3，因为n%base=3.
        return false;

    while (n) 
	{
        if (n%base < 3 || n%base > 6)//按此进制分解，出范围则为不可行 
            return false;
        n /= base;
    }
    ans.insert(base);//这种进制可行，加入答案中 
    return true;
}

void dfs (int d, ll s, ll u) //0，1，n为原始数 
{

    judge(u/s, s);//s是进制，判断这种进制可不可以。 

    if (d == c)//c为不同的质因数个数 
        return;//不同的质因数构成base 

    for (int i = 0; i <= cnt[d]; i++)//一个一个搜索质因数 
    {
        dfs (d+1, s, u);//多了一个进制数 ，从一个不乘到乘完，用dfs列举这个数质因数的组合情况倒是挺巧的 
        s *= fact[d];//这是再把进制乘回去咩 
    }
}

void solve (ll n) 
{

    if (n <= 0)//n不符合。return 
        return;

    div_factor(n);//分解质因数 ，n-3=base*k，将枚举n-3的因子作为base， 
    dfs(0, 1, n);
}

int main () {
    prime_table(maxn);

    int cas;
    scanf("%d", &cas);

    for (int kcas = 1; kcas <= cas; kcas++)
	 {
        //scanf("%lld", &N);
        scanf("%I64d", &N);
        printf("Case #%d: ", kcas);//都喜欢在这个位置输出case啊 
        if (N >= 3 && N <= 6)
            printf("-1\n");//3-6的时候任意进制都可以，直接输出-1，表示无穷 
        else 
		{
            ans.clear();//初始化 
            for (int i = 3; i <= 6; i++)
			{
                nowbase = i;
                solve(N-i);//当前位为i,n-3=base*k,分解n-3找base 
            }
            int ret = ans.size();
            printf("%d\n", ret);
        }
    }
    return 0;
}