HDU 5976 Detachment 【逆元+打表】

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5976

题意：一维的大小可以想象成一条线的长度，二维的大小可以当成矩形的面积，三维的大小可以当成长方体的体积，以此类推，四维就是4个长度相乘，五维就是五个长度相乘，问最大大小是多少；

思路：

比较容易猜测的情况就是，一个线段尽可能的分成多个，分成的小线段两两之间差值尽可能小，而且最短的长度尽可能长；比如线段 9 的最佳的切分方式是 2 3 4 ，相应大小是 2*3*4；

如果按照上面的规则切分线段会发现可能存在的情况有3种；

第一种是刚好可以分成 2 3 4 .... x，相邻之间相差1；

第二种情况是 2 3 4 ... x  总长度比题目给出的长度多1，这种情况的处理方式是把2去掉，得到的总长度比题目给出的长度少1，把这个1给最大的值x，那么大小就是 3*4*...x-1* x+1；

第三种情况就是总长度比题目给的长度多 k 而且这个k一定是在 x以内的值，那么在计算大小的时候丢掉长度为k的线段；

因为在取模的过程中需要除去第二种或者第三种情况的一些数据，在累乘的取模结果中要想除掉部分数据，需要用到逆元；最后需要注意的是一个时间的问题，对于累乘的模和逆元的模（因为过程中存在大量的乘和取模的操作），一定要事先用数组存起来，不然很容易超时；

//#include <ext/pb_ds/assoc_container.hpp>
//#include <ext/pb_ds/tree_policy.hpp>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<queue>
#include<map>
#include<stack>
#include<sstream>
#include<vector>
#include<string>
#include<set>
 
using namespace std;
//using namespace __gnu_pbds;
 
#define IOS ios::sync_with_stdio(false); cin.tie(0);
#define REP(i,n) for(int i=0;i<n;++i)
 
int read(){
 
    int r=0,f=1;char p=getchar();
    while(p>'9'||p<'0'){if(p=='-')f=-1;p=getchar();}
    while(p>='0'&&p<='9'){r=r*10+p-48;p=getchar();}return r*f;
}
 
//typedef tree<pair<long long,int>,null_type,less< pair<long long,int> >,rb_tree_tag,tree_order_statistics_node_update> rbtree;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<long long,long long> pll;
const int Maxn = 1e7+5;
const long long LINF = 1e18;
const int INF = 0x3f3f3f3f;
const int Mod = 1e9+7;
const double PI = acos(-1.0);
 
int pow_mod (int a,int n) {  // 快速幂 
    if(n == 0) return 1;
    int x = pow_mod (a, n>>1);
    ll ans = (ll)x*x%Mod;
    if(n&1) ans = ans*a%Mod;
    return (int)ans;
}
 
ll A[Maxn],MOD[Maxn],n_mod[200005],sum[Maxn];
 
void init () {
	MOD[1] = 1;
	for (int i = 2; i <= 2e5; ++i) {
		sum[i] = sum[i-1]+i;   // 累加和 
		MOD[i] = (MOD[i-1]*i)%Mod; // 累乘模 
		n_mod[i] = pow_mod(i,Mod-2); // 逆元 
	} 
}
 
int main (void)
{
    int T,n;
    init ();
    T = read();
    while (T--) {
    	n = read();
        if(n <= 4) printf("%d\n",n);
        else if(n < Maxn && A[n]) printf("%lld\n",A[n]);
        else {
            int p = lower_bound (sum,sum+100000+1,n)-sum;  // 通过累加和找到对应的MOD值 
            ll ret = MOD[p],tmp;
          	if(sum[p]-n == 1) {  
          		tmp = n_mod[p*2];  // 注意下面的long long 一定要加
          		ret = ((ll)ret*tmp%Mod*(p+1)%Mod)%Mod; // 把 p和2都除掉，然后在乘p+1的模 ，p值对应的是累加过程中的最后一个数x 
			} else if (sum[p]-n > 1) {
          		tmp = n_mod[sum[p]-n];
				ret = ((ll)ret*tmp)%Mod;
			}
            printf("%lld\n",ret);
            if(n < Maxn) A[n] = ret;
        }
    }
    return 0;
}