DP 路径

36 篇文章 0 订阅
36 篇文章 0 订阅

题面去内网找

每一个坏点对答案的负贡献就是走到他的总方案数×他到(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 100005
using 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;
}
  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值