bzoj4767: 两双手

来源：互联网发布：linux bond 模式编辑：程序博客网时间：2024/04/28 14:56

注意一个条件，ax*by!=bx*ay，所以每个点所需的每只手的次数是固定的，于是就可以转换为网格图（注意这里|坐标|<=500*500）中每次向上或向右走一步。因为障碍点较少，可以计算无法到达的方案数=先走到某一个障碍，再随便走。复杂度O（n*n）。

注意有一些无解的可能需要特判。//发现bzoj上使用cerr会返回RE 2333

#include<iostream>#include<cstdio>#include<algorithm>#define P 1000000007#define ll long longusing namespace std;int Ex,Ey,n,ax,ay,bx,by,b[505],inv[1000005],p[1000005],cnt;struct T{int x,y;}a[505],t;bool cmp(T a,T b){return a.x!=b.x?a.x<b.x:a.y<b.y;}int C(int x,int y){if (x==0) return 1;return (ll)p[y]*inv[y-x]%P*inv[x]%P;}int gcd(int x,int y){return y?gcd(y,x%y):x;}T Get(int x,int y){T ans;if ((x*by-y*bx)%(ax*by-ay*bx)) return (T){666666666,666666666};ans.x=(x*by-y*bx)/(ax*by-ay*bx);if ((-x*ay+y*ax)%(ax*by-ay*bx)) return (T){666666666,666666666};ans.y=(-x*ay+y*ax)/(ax*by-ay*bx);if (ans.x<0||ans.y<0) return (T){666666666,666666666};return ans;}int main(){freopen("1.in","r",stdin);freopen("1.out","w",stdout);scanf("%d%d%d",&Ex,&Ey,&n);scanf("%d%d%d%d",&ax,&ay,&bx,&by);p[0]=1;p[1]=1;inv[0]=1;inv[1]=1;for (int i=2;i<=1000000;i++)p[i]=(ll)p[i-1]*i%P,inv[i]=(ll)(-P/i)*inv[P%i]%P;for (int i=2;i<=1000000;i++)inv[i]=(ll)inv[i]*inv[i-1]%P;//A*ax+B*bx=Ex//A*ay+B*by=Eya[cnt=1]=Get(Ex,Ey);if (a[1].x==666666666||a[1].x<0||a[1].y<0){puts("0");return 0;}for (int i=1,x,y;i<=n;i++){scanf("%d%d",&x,&y);t=Get(x,y);if (t.x<=a[1].x&&t.y<=a[1].y) a[++cnt]=t;}sort(a+1,a+cnt+1,cmp);for (int i=1;i<=cnt;i++){//cerr<<a[i].x<<' '<<a[i].y<<endl;b[i]=C(a[i].x,a[i].x+a[i].y);//cout<<i<<' '<<b[i]<<endl;//cerr<<b[i]<<endl;for (int j=1;j<i;j++)if (a[j].x<=a[i].x&&a[j].y<=a[i].y)b[i]=(b[i]-(ll)b[j]*C(a[i].x-a[j].x,a[i].x+a[i].y-a[j].x-a[j].y)%P)%P;}printf("%d\n",(b[cnt]+P)%P);//system("pause");}

0 0