BZOJ 4767 两双手 题解

题目概括

\qquad 就是在矩阵中跳马，给定两种跳马的方法$(x+Ax,y+Ay)$和$(x+Bx,y+By)$，矩阵大小在$(-500,-500)\to (+500,+500)$，然后给个终点，还有$n$个不能跳的点，求方案数嘛。

整理

初始思路

\qquad 看完题面，首先想到的是$dp$，毕竟某道只能向右、向下走题和这题有着不可描述的关联，题面看着就像。
\qquad 考虑到走法只有两种，那么就把它们当成向右或向下就好了，那从原点到某点的次数是唯一确定的（感觉说不清啊，用图解可还行）
像这张，这是最简单的走法：

然后看这个，假设两种箭头，其实从原点到某点的次数也是唯一确定的：

这下可以理解把走法转化成简单的向下、向右了吧。
\qquad 把网格旋转扭曲一下，格子可能是四边形的，但不妨碍它已经是个矩阵了，把每个特殊点和终点按照两种走法转化一下……然而这样很麻烦。
\qquad 再想一下，直接开始$dp$就好了，碰到特殊点避过就好了，用$Ax,Ay$作为纵轴，一个一个扫过去，然后只要判断到$(x-Ex)/Bx$和$(y-Ey)/By$都是整数且相等的时候就切到下一行，当然每一行的开头都要标记，$dp$到终点就完事了。

这是在没构代码的时候觉得完事了，构完代码之后想到的更多了，才发现，要是走法是右两格和右一格呢，那些斜率相同的走法是没办法用这种$dp$的，下面展示一下错误代码，我快疯了，莫名RE

错误代码

#include<bits/stdc++.h>
using namespace std;
//#define read(a){read(a),a+=500;}
inline int read(int &a){int f=1;char c=getchar();if(c==EOF) return a=EOF;
for(;c<'0' || c>'9';c=getchar()){if(c==EOF) return a=EOF;if(c=='-') f=-1;}
for(;c>='0' && c<='9';a=a*10+c-'0',c=getchar());return a*=f;}
const int mo=1e9+7;
int ex,ey,n;
int ax,ay,bx,by;
int xy[1010][1010];
int f[2010][2010];

int main()
{
	read(ex),read(ey),read(n);
	ex+=500,ey+=500;//因为有负坐标，全部加500 
	read(ax),read(ay),read(bx),read(by);
//	ax+=500,ay+=500,bx+=500,by+=500;
	for(int i=1;i<=n;i++)
	{
		int sx,sy;
		read(sx),read(sy);
		xy[sx+500][sy+500]=1;
	}
	int i=1,j=1,x=500,y=500,tx=x,ty=y;//原点是(500,500)，(tx,ty)是行的开头坐标 
	f[i][j]=1;
	while(x!=ex && y!=ey)
	{
//		cout<<i<<' '<<j<<' '<<x<<' '<<y<<endl;
		if(bx==0 && x==ex) i++,j=1,tx+=bx,ty+=by,x=tx,y=ty;//可以直接通过B走法到达终点就换行 
		else if(by==0 && y==ey) i++,j=1,tx+=bx,ty+=by,x=tx,y=ty;
		else if((x-ex)/bx==(y-ey)/by) i++,j=1,tx+=bx,ty+=by,x=tx,y=ty;
		else j++,x+=ax,y+=ay;//没到边界列，那就接着走 
		if(x<0 || x>1000 || y<0 || y>1000) f[i][j]=(1LL*f[i-1][j]+f[i][j-1])%mo;//这里是很朴素的dp 
		else if(xy[x][y]==0) f[i][j]=(1LL*f[i-1][j]+f[i][j-1])%mo;
		else f[i][j]=0;
	}
	printf("%d\n",f[i][j]);
	return 0;
}

正解的思路

\qquad 有了上面的解释，正确的解法就可以十分简洁了。
\qquad 考虑到没有障碍的情况，从原点到$(a,b)$的方案数是$C_{n+m}^n$，那有障碍的情况要用到容斥，先把特殊点排好序，然后一个一个搜，设个$f[i]$表示从原点到第$i$个特殊点不经过其他特殊点的方案数，$g[i][j]$表示从$i$到$j$的所有路径条数，那么$f[i]=g[0][i]-{\sum }_{j=1}^{i-1}g[j][i]\ast f[j]$。

代码就这个样子吧：

#include <bits/stdc++.h>
#define ll long long
using namespace std;
int n,m,_x,_y;
int ax,ay,bx,by;
int tot;
ll num[1001000],sum[1001000];//num[i]表示i!,sum[i]表示i!的逆元 
ll mod=1e9+7;
struct eto
{
	int x,y;
}e[510];
bool cmp(eto a,eto b){return a.x==b.x?a.y<b.y:a.x<b.x;}
ll work(ll x,ll y)//快速幂求逆元
{
	ll cnt=1;
	for(;y>0;)
	{
		if(y%2)cnt=cnt*x%mod;
		x=x*x%mod;
		y>>=1;
	}
	return cnt;
}
ll C(ll x,ll y){return num[x]*sum[y]%mod*sum[x-y]%mod;}
void add(int &x,int &y)//求从（0,0）走到该点分别要使用两种走法的次数 
{
	ll x1,x2,y1,y2;
	x1=x*by-y*bx;x2=ax*by-ay*bx;
	y1=x*ay-y*ax;y2=bx*ay-by*ax;
	if(x2==0||y2==0||x1/x2*x2!=x1||y1/y2*y2!=y1)x=-1,y=-1;
	else x=int(x1/x2),y=int(y1/y2);
	return;
}
ll f[1001000];
int main()
{
	freopen("hand.in","r",stdin);
	freopen("hand.out","w",stdout);
	scanf("%d%d%d",&_x,&_y,&n);
	scanf("%d%d%d%d",&ax,&ay,&bx,&by);
	add(_x,_y);
	int x,y;
	for(int i=1;i<=n;i++)
	{
		scanf("%d%d",&x,&y);
		add(x,y);
		if(x!=-1&&x<=_x&&y<=_y)
		{
			e[++tot].x=x;
			e[tot].y=y;
		}
	}
	e[++tot].x=_x;
	e[tot].y=_y;
	sort(e+1,e+tot+1,cmp);
	num[0]=1;sum[0]=1;
	for(int i=1;i<=1000000;i++)num[i]=num[i-1]*i%mod;
	sum[1000000]=work(num[1000000],mod-2);
	for(int i=999999;i>=1;i--)sum[i]=sum[i+1]*(i+1)%mod;
	for(int i=1;i<=tot;i++)
	{
		f[i]=C(e[i].x+e[i].y,e[i].x);//从（0,0）直接走到第i个特殊点不计前i-1个障碍的总方案数 
		if(f[i])
			for(int j=1;j<i;j++)
				f[i]=(f[i]-f[j]*C(e[i].x-e[j].x+e[i].y-e[j].y,e[i].x-e[j].x)%mod+mod)%mod;
	}
	printf("%lld",f[tot]);
	return 0; 
}