HDU 5976 Detachment

#include<stdio.h>
#include <vector>
#include <string.h>
#include<iostream>
using namespace std;
const long long mod=1000000007;
long long k,a[100000],sum[100000];
/*
和正确答案对比了几百组数据，思路也是对的，不知道为啥wa了，恶心死了
*/
void init()
{
    k=0;
    a[0]=0;
    sum[0]=1;
    for(int i=2;i<=100000;i++)
    {
        k++;
        a[k]=a[k-1]+i;
        sum[k]=sum[k-1]*i%mod;
        if(a[k]>=1000000007)break;
    }
}
long long   ni(long long  x,long long y)
{
    long long sum=1;
    while(y)
    {
          if(y&1)sum=(sum*x)%mod;
        x=(x*x)%mod;
        y/=2;
    }
    return sum;
}
int main()
{

 init();
 long long b[10]={0};
 b[1]=1;
 b[2]=2;
 b[3]=3;
 int t;
 long long n;
 scanf("%d",&t);
 while(t--)
 {
     scanf("%I64d",&n);
     if(n<=3)printf("%I64d\n",b[n]);
     else
     {
         int l=1,r=k;
         while(l<=r)
         {
             int mid=(r+l)/2;
             if(a[mid]<n)
                l=mid+1;
             else r=mid-1;
         }
         while(a[l]>n)
         {
             l--;
         }
         if(a[l]==n)printf("%I64d\n",sum[l]);
         else
         {int p=n-a[l];
         if(p<=l)
         printf("%I64d\n",((sum[l]*(l+2)%mod)*ni(l-p+2,mod-2))%mod);
         else  printf("%I64d\n",((sum[l]*(l+3)*ni(l-p+2+1,mod-2))%mod));
         }
     }
 }
    return 0;
    }