HDU 5976 Detachment

链接

• HDU 5976 Detachment - https://acm.hdu.edu.cn/showproblem.php?pid=5976

问题

将正整数 $n$ 分拆成若干不等的正整数 $a_i$，使得 $s={\prod }_{i=1}^r{a}_i$最大，求最大值对的$10^9+7$模

分拆方案

分拆方案举例

• 例1：$n=14$，分拆方案：
  • 2、3、4、5 共4项，$s=2\ast 3\ast 4\ast 5=120$
  • 3、5、6 共3项，$s=3\ast 5\ast 6=90$
  • $\cdots \cdots$
• 例2：$n=16$，分拆方案：
  • 2、3、5、6 共4项，$s=2\ast 3\ast 5\ast 6=180$
  • 2、3、4、7 共4项，$s=2\ast 3\ast 4\ast 7=168$
  • $\cdots \cdots$
• 例3：$n=18$，分拆方案：
  • 3、4、5、6 共4项，$s=3\ast 4\ast 5\ast 6=360$
  • 4、5、9 共3项，$s=4\ast 5\ast 9=180$
  • $\cdots \cdots$
• 例4：$n=19$，分拆方案：
  • 3、4、5、7 共4项，$s=2\ast 3\ast 4\ast 5=120$
  • 3、5、6 共3项，$s=3\ast 5\ast 6=90$
• 例5：$n=20$，分拆方案：
  • 2、3、4、5、6 共5项，$s=2\ast 3\ast 4\ast 5\ast 6=720$
  • $\cdots \cdots$

分拆方案总结

• 分拆方案不含1（贪心思想：1对积无贡献）
• 对于给定的n，可以得到不含1的最大的分拆数k及其（唯一）方案（贪心思想：每一项尽量小）
  $$\sum _{i=2}^{k+1}i\le n<\sum _{i=2}^{k+2}i或n=\sum _{i=2}^{k+1}i+r，0\le r\le k+1，即：$$
  $$n=\sum _{i=2}^{k+1}i+r=\{\begin{array}{ll}\stackrel{k项}{\overbrace{2+3+\cdots +k+(k+1)}},& \text{当 r=0}\\ \stackrel{k-r项}{\overbrace{2+\cdots +(k+1-r)}}+\stackrel{r项}{\overbrace{(k+3-r)+\cdots +(k+2)}},& \text{当 1≤r<k}\\ \stackrel{k项}{\overbrace{3+4+\cdots +(k+1)+(k+2)}},& \text{当 r=k}\\ \stackrel{k-1项}{\overbrace{3+4+\cdots +(k+1)}}+(k+3),& \text{当 r=k+1}\end{array}$$
  特别的，当 $n=4$ 时，则 $k=1$，即分拆方案为 { 4 } \{4\} {4}
• 有不含1的最大的分拆数k时，积 $s$ 最大。证明如下：
  • 设 $n=a_1+a_2+\cdots +a_k且a_1<a_2<\cdots <a_k且a_1=2或3或4$，是不含1的最大分拆数方案
  • 设 $n=b_1+b_2+\cdots +b_j且1<b_1<b_2<\cdots <b_j且j<k$，是不含1的分拆方案，则
    $$\prod _{i=1}^j{b}_i<\prod _{i=1}^j{a}_i，即最大分拆数方案对应的积s最大$$
  • 规则1：若 $3<a_1<a_2$，则 $2\cdot (a_1-1)\cdot (a_2-1)>a_1\cdot a_2$.
    $$\begin{gathered} \because 2\cdot (a_1-1)\cdot (a_2-1)-a_1\cdot a_2 \\ =(a_1-2)\cdot a_2-2a_1+2 \\ \ge (a_1-2)\cdot a_2-2a_2+2 \\ =(a_1-4)\cdot a_2+2\ge 2 \\ \therefore 2\cdot (a_1-1)\cdot (a_2-1)>a_1\cdot a_2，即将2项拆成3项 \end{gathered}$$
  • 规则2：若 $4\le a$，则 $a\le 2(a-2)$
    $$\begin{gathered} \because 2(a-2)-a=a-4\ge 0 \\ \therefore a\le 2(a-2)，即将1项拆成2项 \end{gathered}$$
  • 规则3：若 $1\le a_1<a_2$，则 $a_1\ast a_2\le (a_1+1)(a_2-1)$
    $$\begin{gathered} \because (a_1+1)(a_2-1)-a_1\cdot a_2=a_2-a_1-1\ge 0 \\ \therefore (a_1+1)(a_2-1)\ge a_1\cdot a_2 \end{gathered}$$
  • 通过规则1-3，将 { b 1 , ⋯   , b j } \{b_1,\cdots,b_j\} {b1​,⋯,bj​}方案转化为 { a 1 , ⋯   , a k } \{a_1,\cdots,a_k\} {a1​,⋯,ak​}方案，转化过程中积 $s$ 不变小

算法分析

• 先计算k，利用 $2，3，4，5，\cdots$ 的前缀和和二分搜索计算 $k$
• 计算$2\cdots k\%m$
• 分情况计算

代码

• 扩展欧几里得求逆元

// hdu 5976 Detachment
#include<bits/stdc++.h>
using namespace std;
#define MOD 1000000007
#define MAX 45000
int K[MAX], M[MAX];
// ax + by = (a, b)
int ex_gcd(int a, int b, int &x, int &y){
	int d;
	// a = (a, 0) = (a, b)
	if(b == 0){
		x = 1, y = 0;
		return a;
	}
	// b(y+a/bx) + a%bx = (b, a%b)
	d = ex_gcd(b, a%b, x, y); 
	int tmp = x;
	x = y;
	y = tmp - a/b*y;
	return d;
}
int main(){
	int t, n, k, r, ans, x, y;
	K[1] = 2, M[1] = 2;
	for(int i=2; i < MAX; ++i) K[i] = K[i-1]+i+1;
	for(int i=2; i < MAX; ++i) M[i] = 1LL*M[i-1]*(i+1)%MOD;
	scanf("%d", &t);
	while(t--){
		scanf("%d", &n);
		if(n < 5){
			printf("%d\n", n);
			continue;
		}
		k = lower_bound(K+1, K+MAX, n)-K;
		if(K[k] > n) k--;
		r = n - K[k];
		if(r == 0) ans = M[k];
		else if(r == k){
			ans = M[k+1];
			ex_gcd(2, MOD, x, y);
			ans = 1LL * ans * x % MOD;
		}
		else if(r == k+1){
			ans = M[k];
			ex_gcd(2, MOD, x, y);
			ans = 1LL * ans * x % MOD;
			ans = 1LL * ans * (k+3) % MOD;
		}
		else{
			ans = M[k+1];
			ex_gcd(k+2-r, MOD, x, y);
			ans = 1LL * ans * x % MOD;
		}
		printf("%d\n", (ans+MOD)%MOD);
	}
    return 0;
}

• 线性求逆元

// hdu 5976 Detachment
#include<bits/stdc++.h>
using namespace std;
#define MOD 1000000007
#define MAX 45000
int K[MAX], M[MAX], inv[MAX];
int main(){
	int t, n, k, r, ans, x, y;
	K[1] = 2, M[1] = 2, inv[1]=1;
	for(int i=2; i < MAX; ++i) K[i] = K[i-1]+i+1;
	for(int i=2; i < MAX; ++i) M[i] = 1LL*M[i-1]*(i+1)%MOD;
	for(int i=2; i < MAX; ++i) inv[i] = 1LL*(MOD-MOD/i)*inv[MOD%i]%MOD;
	scanf("%d", &t);
	while(t--){
		scanf("%d", &n);
		if(n < 5){
			printf("%d\n", n);
			continue;
		}
		k = lower_bound(K+1, K+MAX, n)-K;
		if(K[k] > n) k--;
		r = n - K[k];
		if(r == 0) ans = M[k];
		else if(r == k){
			ans = M[k+1];
			ans = 1LL * ans * inv[2] % MOD;
		}
		else if(r == k+1){
			ans = M[k];
			ans = 1LL * ans * inv[2] % MOD;
			ans = 1LL * ans * (k+3) % MOD;
		}
		else{
			ans = M[k+1];
			ans = 1LL * ans * inv[k+2-r] % MOD;
		}
		printf("%d\n", (ans+MOD)%MOD);
	}
    return 0;
}