hdu 5976 Detachment

http://acm.hdu.edu.cn/showproblem.php?pid=5976

拆数字，并求所拆的数字乘积的最大值，显然是拆成2*3*4*5*。。。*n这样的形式，但会存在

多出来的数字，如果多出来的数字加上2能够等于n+1，则答案就是2*3*。。。*n/2*(n+1),

如果小于n+1,则可以将多出来的数加到前面的3或者4或者5或者后面的数上,是的所加的那个数加上多出来的数字后正好等于n+1,则答案就是2*3*4*5*....*n/（i）*n+1，可以预处理打表前缀和和前缀积，并且除法取模需要用到乘法逆元。

#include<iostream>
#include<cstdio>
using namespace std;
long long int MOD=7+1e9;
int
4000
qzh[55555];
long long int qzj[55555];
int ny[55555];
void ycl()
{
qzj[1]=1;
qzh[1]=0;
ny[1]=1;
for(int i=2;i<=48888;i++)
{
ny[i]=(MOD-MOD/i)*ny[MOD%i]%MOD;
qzj[i]=(i*qzj[i-1])%MOD;
qzh[i]=i+qzh[i-1];
}
}
int main()
{

int t;long long int ans;	int a;
int lef,rig,mid;
while(scanf("%d",&t)!=EOF)
{
ycl();
while(t--)
{

scanf("%d",&a);
if(a<5)printf("%d\n",a);
else {
lef=2;
rig=48888;
mid=(lef+rig)>>1;
while(lef+1<rig)
{
if(qzh[mid]>a)
{
rig=mid;
mid=(lef+rig)>>1;
}
else{
lef=mid;
mid=(lef+rig)>>1;
}
}
int dy=a-qzh[lef];
if(2+dy>=lef+1)ans=qzj[lef]*ny[2]%MOD*(2+dy)%MOD;
else ans=qzj[lef]*ny[lef+1-dy]%MOD*(lef+1)%MOD;
ans=(ans+MOD)%MOD;
printf("%lld\n",ans);
}
}
}
return 0;
}