【数论】【逆元】【贪心】HDU5976 Detachment

最新推荐文章于 2026-06-08 09:35:50 发布

原创最新推荐文章于 2026-06-08 09:35:50 发布 · 468 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

数论同时被 2 个专栏收录

3 篇文章

订阅专栏

贪心

1 篇文章

订阅专栏

本文介绍了一种算法问题，即如何将一个整数x分解为若干整数之和时，这些整数的乘积最大。文章给出了具体的解题思路与实现代码，包括使用预处理前缀积和扩展欧几里得算法进行模运算。

题目链接

题意:

给出一个数x，存在 $\sum_{i=1}^n a_i= x$ ， $a_i != a_j$ 使得 $\prod_{i=1}^n$ 最大，求 $\prod_{i=1}^n$ 的最大值。

解题思路：

容易发现要使得积最大，要把x分得尽可能小。所以把x分成2，3，4，5,…,i的形式。从2开始是因为分成1并不会使积更大。如果x能直接分成这个形式，就可以直接得出答案。如果多出了k,则把后k个数后移一个单位，其中后移的时候要考虑k大于已存在的数的个数的情况。
预处理前缀积来节省时间。考虑到前缀积相除不好取模，所以用后一个积乘以前一个积的逆元。

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<string>
#include<algorithm>
#define MAXN 103000
typedef long long ll;
const static ll MOD = 1000000007;
using namespace std;
ll sum[MAXN];
ll s[MAXN];
ll re[MAXN];
ll extend_euclid(ll a,ll b,ll &x,ll &y){
    if(!b){
        x = 1;
        y = 0;
        return a;
    }
    ll d = extend_euclid(b,a%b,x,y);
    ll t = x;
    x = y;
    y = t - a / b * y;
    return d;
}

int main(){
    ll T,x,tmp,ans;
    ll t = 2;
    s[0] = sum[0] = 2;
    extend_euclid(s[0],MOD,re[0],tmp);
    re[0] = (re[0] + MOD) % MOD;
    for(int i = 1;i<MAXN;i++){
        sum[i] = ++t + sum[i-1];
        s[i] = t * s[i-1]%MOD;
        extend_euclid(s[i],MOD,re[i],tmp);
        re[i] = (re[i] + MOD) % MOD;
    }
    scanf("%I64d",&T);
    while(T--){
        scanf("%I64d",&x);
        if(x < 5){
            printf("%I64d\n",x);
            continue;
        }
        int ind = lower_bound(sum,sum+MAXN,x)-sum;
        if(sum[ind] == x){
            printf("%I64d\n",s[ind]);    
            continue;
        }
        ind--;
        ll k = x - sum[ind];
        if(k != ind + 2){
            ans = s[ind + 1] * (re[ind-k+1]) % MOD;
            if(ind - k >= 0)
                ans = ans * s[ind - k] % MOD;

        }else{
            ans = s[ind] * re[0] % MOD;
            ans = ans * (ind + 4) % MOD;
        }
        printf("%I64d\n",ans);
    }
    return 0;
}