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  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;}