跳跃【题解】

---

难题直接摆烂

不考虑 $f_i$ 的限制

假设不考虑 $f_i$ 的限制，考虑用最简单的 $DP$ 算法：

设 $F[r][x][y]$ 表示跳 $r$ 步后到达 $(x,y)$ 的方案数。

F [ r ] [ x ] [ y ] = s u m { F [ r − 1 ] [ x ′ ] [ y ′ ] } ( 0 ≤ x − x ′ ≤ M x , 0 ≤ y − y ′ ≤ M y ) F[r][x][y]=sum\{F[r-1][x'][y']\}(0≤x-x'≤Mx,0≤y-y'≤My) F[r][x][y]=sum{F[r−1][x′][y′]}(0≤x−x′≤Mx,0≤y−y′≤My)

但这样会超时，我们可以把 $2$ 个状态 $x$ 和 $y$ 拆出来，形成 $Fx$ 和 $Fy$ 的 $2$ 个动态规划数组。

F x [ r ] [ x ] = s u m { F x [ r − 1 ] [ x ′ ] } ( 0 ≤ x − x ′ ≤ M x ) Fx[r][x]=sum\{Fx[r-1][x']\}(0≤x-x'≤Mx) Fx[r][x]=sum{Fx[r−1][x′]}(0≤x−x′≤Mx) F y [ r ] [ y ] = s u m { F y [ r − 1 ] [ y ′ ] } ( 0 ≤ y − y ′ ≤ M y ) Fy[r][y]=sum\{Fy[r-1][y']\}(0≤y-y'≤My) Fy[r][y]=sum{Fy[r−1][y′]}(0≤y−y′≤My)

---

0 的限制

我们在计算的时候，要保证 $(x,y)$ 和 $(x^{\prime },y^{\prime })$ 是不同的 $2$ 个点，这时想到，

可以直接令 $f_n$${}_{+}$${}_1=0$

---

考虑 $f_i$ 的限制

我们知道：正确答案 $=$ 不考虑限制的答案 $-$ 不合法的答案。

考虑跳跃的过程，可以得到等式：

$$x_1+x_2+x_3+...+x_R=Tx$$$$y_1+y_2+y_3+...+y_R=Ty$$

其中 $x_1,x_2,...,x_R,y_1,y_2,...,y_R$ 为每次跳跃的距离。

只要有一对 $(x_i,y_i)(1\le i\le R)$ 符合 $x_i=y_i$ 且为 $f$ 中的一个数，那么整个跳跃的过程都不合法。

---

容斥原理

设 $dp[i][j]$ 表示至少有 $i$ 对 $(x,y)$ 不合法，其不合法的 $f$ 值得和为 $j$ 时的方案数。

d p [ i ] [ j ] = s u m { d p [ i − 1 ] [ j − f k ] } ( f k ≤ j , 1 ≤ k ≤ n ) dp[i][j]=sum\{dp[i-1][j-f_k]\}(f_k≤j,1≤k≤n) dp[i][j]=sum{dp[i−1][j−fk​]}(fk​≤j,1≤k≤n)

设 $Fis[i]$ 表示至少有 $i$ 对 $(x,y)$ 不合法的方案数。

F i s [ i ] = C ( n , i ) ∗ s u m { d p [ i ] [ j ] ∗ F x [ R − i ] [ T x − j ] ∗ F y [ R − i ] [ T y − j ] } ( 0 ≤ j ≤ m i n { T x , T y } , 1 ≤ i ≤ R ) Fis[i]=C(n,i)*sum\{dp[i][j]*Fx[R-i][Tx-j]*Fy[R-i][Ty-j]\}(0≤j≤min\{Tx,Ty\},1≤i≤R) Fis[i]=C(n,i)∗sum{dp[i][j]∗Fx[R−i][Tx−j]∗Fy[R−i][Ty−j]}(0≤j≤min{Tx,Ty},1≤i≤R)

根据容斥原理，不合法总方案数为：

$$Fis[1]-Fis[2]+Fis[3]-Fis[4]+...+(-1{)}^R\ast Fis[R]$$

---

#include<bits/stdc++.h>
using namespace std;
int Tx,Ty,Mx,My,R,n,ans;
int bad[60];
int fx[1610][810],fy[1610][810];
int s[810];
int dp[1610][810],fis[1610];
int tri[1610][1610];
void DP()
{
	fx[0][0]=fy[0][0]=1;
	for(int i=1;i<=R;i++)
	{
		memset(s,0,sizeof s);
		for(int j=0;j<=Tx;j++)
		{
			if(j)s[j]=s[j-1];
			s[j]=(s[j]+fx[i-1][j])%10007;
			fx[i][j]=s[j]%10007;
			if(j>Mx)
			{
				fx[i][j]=(fx[i][j]-s[j-Mx-1]+10007)%10007;
			}
		}
		memset(s,0,sizeof s);
		for(int j=0;j<=Ty;j++)
		{
			if(j)s[j]=s[j-1];
			s[j]=(s[j]+fy[i-1][j])%10007;
			fy[i][j]=s[j]%10007;
			if(j>My)
			{
				fy[i][j]=(fy[i][j]-s[j-My-1]+10007)%10007;
			}
		}
	}
	ans=(fx[R][Tx]*fy[R][Ty])%10007;
	dp[0][0]=1;
	for(int j=1;j<=R;j++)
	{
		for(int k=0;k<=min(Tx,Ty);k+=10)
		{
			for(int i=1;i<=n;i++)
			{
				if(k>=bad[i])dp[j][k]=(dp[j][k]+dp[j-1][k-bad[i]])%10007;
			}
		}
	}
	for(int i=1;i<=R;i++)
	{
		for(int j=0;j<=min(Tx,Ty);j+=10)
		{
			fis[i]=(fis[i]+dp[i][j]*fx[R-i][Tx-j]%10007*fy[R-i][Ty-j]%10007+10007)%10007;
		}
		fis[i]=(fis[i]*tri[R+1][i+1])%10007;
		if(i&1)
		{
			ans=(ans-fis[i]+10007)%10007;
		}
		else
		{
			ans=(ans+fis[i])%10007;
		}
	}
}
int main()
{
	tri[1][1]=1;
	for(int i=2;i<=1601;i++)
	{
		for(int j=1;j<=i;j++)
		{
			tri[i][j]=(tri[i-1][j]+tri[i-1][j-1])%10007;
		}
	}
	scanf("%d%d%d%d%d%d",&Tx,&Ty,&Mx,&My,&R,&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&bad[i]);
	}
	n++;
	DP();
	printf("%d",ans);
	return 0;
}