hdu-5976 Detachment 解题报告（乘法逆元、贪心？）

题目描述

题目：https://vjudge.net/problem/HDU-5976
好像有翻译来着，我就不翻译了。

思路分析

这道题初看起来就让我们联想到高中数学的基本不等式：和一定，分解为相同的数时，他们的积最大，即：aa+bb>=2ab，当且仅当a==b时，ab最大。
但是这道题有个要求：互不相等，所以，我们取个较近的情况：他们是阶乘的情况最大，这时这道题首先的贪心策略（其实我不大确定是不是贪心。。求大佬轻喷。）。
之后，并不是拿到每个数都要分解一次。因为这道题1<=T<=1e6，这么多次询问，每次都去分解很费时间，所以需要通过打表，将分解的结果和用到的结果保存下来，因为阶乘有递推关系，所以计算一个新的数的时候，之前算过的数还能利用。另外，超前一下，这道题取余的时候用到了这个情况(a/b)%m,若直接算，误差会比较大，所以需要用到逆元的性质：
(a/b)%m=(ak)%m，k为b（b在这里是某一个mul的元素）模mod的逆。
为了防超时这个也要打表。所以，打表的代码就是：

void make_list()
{
	for(int i=2;i<maxn;i++)
	{
		sum[i]=sum[i-1]+i;
		mul[i]=(mul[i-1]%mod*i%mod)%mod;
		inv[i]=(mod-mod/i)*inv[mod%i]%mod;
	}
}

接下来的难点是怎么找到这个阶乘。
因为分成1对乘积没有贡献，所以从2开始分起来，也就是要分成2+3+4+5+……+n的形式。这个查找建议用stl库中的lower_bound()函数，它采用的是二分查找（返回第一个大于等于它的数的下标），复杂度为logn，速度较快。
情况1：
若存在n，使得：1+2+3+……+n==长度，那么，刚刚好，我们只要输出对应的mul[n]就可以了。
情况2：
若很尴尬的：sum[n]<长度<sum[n+1]（注：此时lower_bound()返回的是n+1），我们就还需要按情况处理：先算出t=长度-sum[n]的值。原因如下：
对于一组阶乘，若要再将一个数分配给这个阶乘中的数，则将这个数分配给小的数，最后的乘积会更大。证明如下：

因此，在保证没有重复的前提下，应该优先分配给最小的数
情况2.1 若2+t>n(2加上t后超过了阶乘内最大的因数，不会重复）

ll result=((mul[id]*inv[2])%mod*(2+t))%mod;

容我偷个懒。。。但注意，这里要用上面逆元的性质，避免精度损失。
情况2.2 若不符合上面的条件：
这时候可以将这个数分给后面的几个数（避免重复）//这里原理其实还有点含糊，之后会再琢磨琢磨

完整代码

#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
const ll maxn=50000;
ll sum[maxn]={0};
ll mul[maxn]={1,1};
ll inv[maxn]={0,1};
void init()
{
	for(int i=2;i<maxn;i++)
	{
		sum[i]=sum[i-1]+i;
		mul[i]=(mul[i-1]%mod*i%mod)%mod;
		inv[i]=(mod-mod/i)*inv[mod%i]%mod;
	}
}

int main()
{
	ll kase;
	scanf("%lld",&kase);
	init();
	while(kase--)
	{
		ll x;
		scanf("%lld",&x);
		if(x==1)
		{
			cout<<1<<endl;
		}
		else
		{
			ll id=lower_bound(sum+2,sum+maxn,x)-sum;
			if(x-sum[id]==0)
			{
				printf("%lld\n",mul[id]%mod);
				continue;
			}
			else
			{
				id--;
				ll t=x-sum[id];
				if(2+t>id)
				{
					ll result=((mul[id]*inv[2])%mod*(2+t))%mod;
					printf("%lld\n",result);
				}
				else
				{
					ll result=((mul[id]*inv[id+1-t])%mod*(id+1))%mod;
					printf("%lld\n",result);
				}
			}
		}
	}
	return 0;
}