2021.07.14【NOIP提高B组】模拟 积木

思路：

首先很容易想到用DP做法
可是超时怎么办
我们可以用组合数，计算一下答案

$code$

#include<iostream>
#include<cstdio>
#include<cmath>

using namespace std;

long long mod=1e9+7;
long long n, a[20010];
long long f[200100];

long long qp(long long x, long long k)
{
	long long ans=1;
	while(k!=0)
	{
		if(k&1)
			ans=ans*x%mod;
		x=x*x%mod;
		k>>=1;
	}
	return ans%mod;
}

long long inv(long long x, long long k)
{
	return qp(x, k)%mod;
}

long long C(long long n, long long m)
{
	if(m>n)
		return 0;
	return f[n]*inv(f[m], mod-2)%mod*inv(f[n-m], mod-2)%mod;
}

long long solve(long long len, long long x, long long y)
{
	long long sum=0, dif=abs(x-y);
	for(long long r=0; r<=len-dif; r++)
	{
		if((len-r+dif)%2==0)
			sum=(sum+C(len, r)%mod*(C(len-r, (len-r+dif)/2)-C(len-r, (len-r+x+y+2)/2)+mod)%mod)%mod;
	}
	return sum%mod;
}

int main()
{
	freopen("brick.in", "r", stdin);
	freopen("brick.out", "w", stdout);
	f[0]=1;
	for(long long i=1; i<=100010; i++) f[i]=f[i-1]*i%mod;
	scanf("%lld", &n);
	for(register long long i=1; i<=n; ++i)
		scanf("%lld", &a[i]);
	if(a[1]>0||a[n]>0)
	{
		printf("0");
		return 0;
	}
	a[1]=0, a[n]=0;
	register long long ans=1;
	for(long long i=2, j=1; i<=n; i++)
	{
		if(a[i]>=0)
		{
			long long k=solve(i-j, a[i], a[j])%mod;
			ans=(ans*k)%mod;
			j=i;
		}
	}
	printf("%lld", ans);
	return 0;
}