数学之数论（同余与逆元）于2020/8/6

扩展欧几里得算法求一元同余方程
即ax于b同余于m，改写成ax+m*y=b，再解方程.

求解逆元的方法

这一段吃过不少亏
队友用的一直是费马小定理求解逆元（要求模数为素数）
当然也可以用扩展欧几里得算法求逆元

hdu5976

贪心+逆元
当然也可以说是构造，构造方法还是很巧妙的，需要一点灵光一闪

详细的推导过程点这里

似乎还要用二分来找k，时间卡得比较紧。

一整个下午花在这条题目上，讲一下最坑的点：存储2…n的乘积的数组（取余）必须用long long，这我也不知道为什么，总之是对拍试了无数次才发现的。
附上AC代码：

#include<bits/stdc++.h>
using namespace std;
int mod=1e9+7; 
void extend_gcd(int a,int b,int &x,int &y)
{
	if(b==0)
	{
		x=1,y=0;
		return;
	}
	extend_gcd(b,a%b,x,y);
	int tmp=x;
	x=y;
	y=tmp-(a/b)*y;
}
int mod_inverse(int a,int m)
{
	int x,y;
	extend_gcd(a,m,x,y);
	return(m+x%m)%m;
}
int inv[100050];
long long a[100050];
int sum[100050];
long long hanshu(int x)
{
	int l=2,r=45000,mid=(l+r)/2;
	while(l+1<r)
	{
		
		if(sum[mid]>x)
		{		
			r=mid;		
			mid=(l+r)/2;	
		}
		else 
		{		
			l=mid;
			mid=(l+r)/2;
		}
	}
	int k=x-sum[l];
//	cout<<"sum[l]:"<<sum[l]<<endl;
//	cout<<"l,mid,k:"<<l<<' '<<mid<<' '<<k<<endl;
	if(2+k>l)
	{
		long long ans=a[l]*inv[2]%mod*(k+2)%mod;
//		cout<<"a[l]:"<<a[l]<<endl;
//		cout<<"inv[2]:"<<inv[2]<<endl;
//		cout<<"ans:"<<ans<<endl;
	//	return (a[l]*inv[2]%mod*(k+2)%mod+mod)%mod;
	return (ans+mod)%mod;
	
	}
	else 
	{
		long long ans=a[l]*inv[l+1-k]%mod*(l+1)%mod;
//		cout<<"inv[l+1-k]:"<<inv[l+1-k]<<endl;
//		cout<<"ans:"<<ans<<endl;
		return (ans+mod)%mod;
		//return (a[l]*inv[l+1-k]%mod*(l+1)%mod+mod)%mod;
	}
}
int main()
{
	a[1]=1;
	inv[1]=1;
	sum[1]=0;
	for(int i=2;i<=50000;i++)
	{
		inv[i]=mod_inverse(i,mod);
	//	inv[i]=(mod-mod/i)*inv[mod%i]%mod;
		a[i]=(a[i-1]*i)%mod;
		sum[i]=sum[i-1]+i;
	}
	int t;
	scanf("%d",&t);
	while(t--)
	{
		int x;
		scanf("%d",&x);
		if(x<5)printf("%d\n",x);
		else
		printf("%d\n",hanshu(x));
	}
}