DP 路径

来源:互联网 发布:p2p网络管理器 编辑:程序博客网 时间:2024/05/21 11:32

题面去内网找

每一个坏点对答案的负贡献就是走到他的总方案数×他到(n,m)的总方案数。
计算方法:f[i]=C(x+y,x)-sigma f[j]*C(x1,y1)
x1,y1就是i,j的x,y之差。而j的要求就是,xj<=xi,yj<=yi。
对坏点排个序,N^2搞即可。

#include<cstdio>#include<cstring>#include<cstdlib>#include<algorithm>#include<iostream>#define mod 1000000007#define ll long long#define N 100005using namespace std;int n,m,k;struct node{int x,y;}s[3005];ll a[N*3],b[N*3],f[3005],ans;ll cheng(ll x,ll m){    ll ans=1;    while(m)    {         if(m&1)ans=ans*x%mod;         x=x*x%mod;         m>>=1;    }    return ans;}void init(){      b[0]=a[0]=1;      for(ll i=1;i<=n+m;i++)a[i]=a[i-1]*i%mod;      b[n+m]=cheng(a[n+m],mod-2);      for(ll i=n+m-1;i>=1;i--)b[i]=b[i+1]*(i+1)%mod;}ll C(int n,int m){return (a[n]*b[m])%mod*b[n-m]%mod;}bool cmp(node x,node y){return x.x!=y.x?x.x<y.x:x.y<y.y;}int main(){    scanf("%d%d%d",&n,&m,&k);    init();    for(int i=1;i<=k;i++)        scanf("%d%d",&s[i].x,&s[i].y);    s[k+1]=(node){n,m};    sort(s+1,s+k+2,cmp);    for(int i=1;i<=k+1;i++)    {        f[i]=C(s[i].x+s[i].y,s[i].x);        for(int j=i-1;j>=1;j--)        {             if(s[j].y>s[i].y)continue;             int x=s[i].x-s[j].x,y=s[i].y-s[j].y;             f[i]-=f[j]*C(x+y,x)%mod;             if(f[i]<0)f[i]+=mod;        }    }     cout<<f[k+1]%mod;}