BZOJ[4767] 两双手

向量是可以被基地唯一分解的，然后就和昨天的T3一样了

在计算x_的时候居然没注意最大可以变成500000，然后最开始设的最大值是100000，然后就一直wa

Code

#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
# define mod 1000000007
using namespace std;
typedef long long LL;
int n,m,K;
int N,M;
int ax,bx,ay,by;
LL jc[500010],njc[500010];
LL ksm(LL a,LL b){
    LL ans=1;
    while(b){
        if(b&1) ans*=a,ans%=mod;
        b=b>>1; a*=a; a%=mod;
    }
    return ans;
}
void beg(){
    int lim=500000; jc[0]=1;
    for(int i=1;i<=lim;i++) jc[i]=jc[i-1]*i%mod;
    for(int i=0;i<=lim;i++) njc[i]=ksm(jc[i],mod-2);
}
LL C(int a,int b){
    LL ret=jc[a]*njc[b]%mod*njc[a-b]%mod;
    return ret;
}
struct node{
    int x,y,x_,y_;
}g[510];
bool cmp(const node a,const node b){
    if(a.x_==b.x_) return a.y_<b.y_;
    return a.x_<b.x_;
}
void init(){
    scanf("%d%d%d",&n,&m,&K);
    scanf("%d%d%d%d",&ax,&ay,&bx,&by);
    for(int i=1;i<=K;i++) scanf("%d%d",&g[i].x,&g[i].y);
    
}
void cal(int x,int y,int &A,int &B){
    int s=x*by-y*bx,t=ax*by-ay*bx;
    if(s%t) A=1000000; else A=s/t;
    if(A<0 || A>N) A=1000000;
    s=x*ay-y*ax,t=bx*ay-by*ax;
    if(s%t) B=1000000; else B=s/t;
    if(B<0 || B>M) B=1000000;
}
LL f[510];
void work(){
    cal(n,m,N,M);
    if(N==1000000 || M==1000000){cout<<0<<endl;return;}
    for(int i=1;i<=K;i++){
        cal(g[i].x,g[i].y,g[i].x_,g[i].y_);
        if(g[i].x_==1000000 || g[i].y_==1000000) g[i].x_=1000000,g[i].y_=1000000;
    }
    sort(g+1,g+K+1,cmp); int ji=K;
    for(int i=1;i<=K;i++){
        if(g[i].x_>N || g[i].y_>M){ji=i-1; break;}
    }
    K=ji;
    LL ans=C(N+M,N);
    for(int i=1;i<=K;i++){
        f[i]=C(g[i].x_+g[i].y_,g[i].x_);
        for(int j=i-1;j>=1;j--){
            if(g[j].y_>g[i].y_) continue;
            f[i]=(f[i]-f[j]*C(g[i].x_-g[j].x_+g[i].y_-g[j].y_,g[i].x_-g[j].x_)%mod)%mod;
        }
    }
    for(int i=1;i<=K;i++){
        ans=(ans-f[i]*C(N-g[i].x_+M-g[i].y_,N-g[i].x_)%mod)%mod;
    }
    if(ans<0) ans+=mod;
    cout<<ans<<endl;
}
int main(){
//    freopen("a.in","r",stdin);
    init();
    beg();
    work();
}