省选模拟赛第一场 T2头目选举（组合计数+拉格朗日插值）

 

 

 

题解

官方题解：

因为一个人只要有一个技能不被碾压，它就可以成为n-j的那一部分人

所以我们就只需要计算是哪些人在第 i 个技能胜过蒟蒻舵主

其中k-j的那一部分人必须胜过舵主，剩下的就需要用已经在前面的某个技能胜过舵主的那一部分人来凑

选好了人，就直接算他们可能的技能伤害的组合（预处理）

有一个易错点，就是拉格朗日插值只能计算多项式，非多项式的函数就只能拟合了（有误差）

而我们一开始的点值不能带入

......

应该代入

.......

因为上面的那个不是多项式函数，下面的才是多项式函数（应该是吧）（被坑了好久。。。）

 

代码：٩(๑>◡<๑)۶人生第2道拉格朗日插值٩(๑>◡<๑)۶

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
inline int gi()
{
	char c;int num=0,flg=1;
	while((c=getchar())<'0'||c>'9')if(c=='-')flg=-1;
	while(c>='0'&&c<='9'){num=num*10+c-48;c=getchar();}
	return num*flg;
}
#define N 105
const int mod=1000000007;
int f[N][N],C[N][N],g[N],H[N],R[N];
int y[N];
int ksm(int x,int y)
{
	int ret=1;
	while(y){
		if(y&1)ret=1ll*ret*x%mod;
		y>>=1;x=1ll*x*x%mod;
	}
	return ret;
}
int main()
{
	freopen("leader.in","r",stdin);
	freopen("leader.out","w",stdout);
	int n,m,p,i,j,k,sum;
	C[0][0]=1;
	for(i=1;i<=100;i++){
		C[i][0]=C[i][i]=1;
		for(j=1;j<i;j++)
			C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
	}
	n=gi();m=gi();p=gi();
	for(i=1;i<=m;i++)H[i]=gi();
	for(i=1;i<=m;i++)R[i]=gi();
	for(i=1;i<=m;i++){
		for(j=1;j<=102;j++){
			sum=0;
			for(k=1;k<=j;k++)
				sum=(1ll*sum+1ll*ksm(j-k,R[i]-1)*ksm(k,n-R[i]))%mod;
			y[j]=sum;
		}
		for(j=1;j<=102;j++){
			sum=y[j];
			for(k=1;k<=102;k++){
				if(k==j)continue;
				sum=1ll*sum*(H[i]-k)%mod*ksm((j-k+mod)%mod,mod-2)%mod;
			}
			g[i]=(1ll*g[i]+1ll*sum)%mod;
		}
		g[i]=(g[i]+mod)%mod;
	}
	memset(f,0,sizeof(f));
	f[0][n-1]=1;
	for(i=1;i<=m;i++)
		for(j=0;j<n;j++)
			for(k=j;k<n;k++)
				if(R[i]-1-(k-j)>=0)
					f[i][j]=(1ll*f[i][j]+1ll*f[i-1][k]*C[k][k-j]%mod*C[n-k-1][R[i]-1-(k-j)]%mod*g[i])%mod;
	printf("%d",f[m][p]);
}