题意,给出一个n(1≤n≤10^9),将n分成若干不相同的部分a1,a2,a3,a4,a5......,s=a1*a2*a3*a4...*am,让求s的最大值。
和相同,这部分数为连续的数的话乘积最大。拆分数以2,3,4,5....,pos,....这样的格式。
当不能恰好这样分完n时:加到pos时和比n大就停止,将pos倒着均分到每个ai上。以此使得连续的数变大。
此时有两种情况:
1,当n=23时,2,3,4,5,6此时sum[6]=20,下一个7不能选择(找出这个pos用二分快),需要向前面均分3 (23-20),得到2 3 5 6 7。
结果就是7的阶乘除以4,(除法取模要用乘它的逆元取模)
2,当n=26时,2,3,4,5,6,sum[6]=20,下一个7不能选择,需要向前面均分6(26-20) ,得到3,4,5,6,8。
结果·为8的阶乘除7,再除2。
模拟这两种情况。
代码:
#include<bits/stdc++.h>
typedef long long ll;
const int mod=1e9+7;
const int maxn=1e5+5;
using namespace std;
ll a[maxn];//前缀和
ll sum[maxn];//前缀积
void extgcd(ll a,ll b,ll&d,ll&x,ll&y)
{
if(!b) {d=a;x=1;y=0;}
else {extgcd(b,a%b,d,y,x);y-=x*(a/b);}
}
ll insverse(ll a,ll n)
{
ll d,x,y;
extgcd(a,n,d,x,y);
return d==1?(x+n)%n:-1;
}
void init()
{
a[1]=0;
sum[1]=1;
for(int i=2;i<=maxn;i++)
{
a[i]=a[i-1]+(ll)i;
sum[i]=(sum[i-1]*i)%mod;//i!
}
}
int main()
{
int t;ll n;
scanf("%d",&t);
init();
while(t--)
{
scanf("%lld",&n);
if(n==1)
{
printf("1\n");
continue;
}
ll ans=0;
ll pos= upper_bound(a,a+maxn,n)-a;//2,3,4,5,6; 23->pos=7
ll dis=n-a[pos-1];//n-a[6]=23-20
///pos+1的阶乘结果除以pos1;
if(dis==pos-1) ans=sum[pos+1]%mod*insverse(pos,mod)%mod*insverse((ll)2,mod)%mod;
else ans=sum[pos]%mod*insverse(pos-dis,mod)%mod;
printf("%lld\n",ans);
}
}