HDU 3758 Factorial Simplification 数论

这题就考一个知识点，就是求阶乘中某个素数因子的个数

题意不再说了

具体思路就是把分子的所有素因子以及其个数存起来。

分母也是。

然后看对应的素因子的个数分子是否都不小于分母就可以判断是不是整数

如果是整数进行下一步

每次二分寻找能整除的最大阶乘的数。直到找不到大于1的数

用二分是因为阶乘的素因子个数一定是非递减的。

二分的上下界多大呢。

下界是2，题目给定了

上界就不好说了。

但是我们注意到10007是素数，而题目的数据范围是10000

也就是说我们上界设为10006就行。因为10000以内的数乘来乘去也不会能模10007去

#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <cstdio>
#include <cmath>
#include <queue>
#include <map>
#include <set>
#define MAXN 100005
#define MAXM 400005
#define INF 1000000011
#define lch(x) x<<1
#define rch(x) x<<1|1
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define eps 1e-7
using namespace std;
bool tag[10001];
int p[10001];
int cnt;
void get_prime()
{
    cnt = 0;
    tag[1] = 1;
    for (int i = 2; i < 10000; i++)
    {
        if (!tag[i])
        p[cnt++] = i;
        for (int j = 0; j < cnt && p[j] * i < 10000; j++)
        {
            tag[i*p[j]] = 1;
            if (i % p[j] == 0)
            break;
        }
    }
}
int n, m;
int a[11005], b[11005];
int numa[10005], numb[10005];
int ans[1000005];
int num[1001111];
void get(int *x, int n, int *num)
{
    for(int i = 1; i <= n; i++)
    {
        for(int j = 0; j < cnt; j++)
        {
            int tn = x[i];
            if(tn < p[j]) break;
            while(tn){num[j] += tn / p[j]; tn /= p[j];}
        }
    }
}
bool ok()
{
    for(int i = 0; i < cnt; i++)
        if(numb[i] > numa[i]) return 0;
    return 1;
}
int main()
{
    int T;
    get_prime();
    scanf("%d", &T);
    while(T--)
    {
        scanf("%d%d", &n, &m);
        for(int i = 1; i <= n; i++) scanf("%d", &a[i]);
        for(int i = 1; i <= m; i++) scanf("%d", &b[i]);
        memset(numa, 0, sizeof(numa));
        memset(numb, 0, sizeof(numb));
        get(a, n, numa);
        get(b, m, numb);
        if(!ok())
        {
            printf("-1\n");
            continue;
        }
        for(int i = 0; i < cnt; i++) numa[i] -= numb[i];
        int tp = 0;
        memset(num, 0, sizeof(num));
        while(1)
        {
            int low = 2, high = 10006;
            int can = 0;
            while(low <= high)
            {
                int mid = (low + high) >> 1;
                int fg = 1;
                for(int j = 0; j < cnt; j++)
                {
                    int tn = mid;
                    if(tn < p[j]) break;
                    int sum = 0;
                    while(tn){sum += tn / p[j]; tn /= p[j];}
                    if(sum > numa[j])
                    {
                        fg = 0;
                        break;
                    }
                }
                if(fg) {low = mid + 1; can = max(can, mid);}
                else high = mid - 1;
            }
            int x = can;
            if(x == 0) break;
            ans[tp] = x;
            int ot = INF;
            for(int j = 0; j < cnt; j++)
            {
                int tn = x;
                if(tn < p[j]) break;
                int sum = 0;
                while(tn){sum += tn / p[j]; tn /= p[j];}
                if(sum > 0) ot = min(ot, numa[j] / sum);
            }
            num[tp] += ot;
            tp++;
            for(int j = 0; j < cnt; j++)
            {
                int tn = x;
                if(tn < p[j]) break;
                int sum = 0;
                while(tn){sum += tn / p[j]; tn /= p[j];}
                numa[j] -= ot * sum;
            }
        }
        printf("%d\n", tp);
        for(int i =0; i < tp; i++) printf("%d %d\n", ans[i], num[i]);
    }
    return 0;
}