問題解決策を数えるパスシミュレーションゲーム

問題の意味:k個の障害は(M、N)に(0,0)を探し、あるはこれらの障壁を介してプログラムの数を通過しません。

直接DPのタイムアウト。バリアストリップを介してパスの数 - バリアストリップを通るパスの数、すなわち、C(N + M、N) - パスの総数の障害のパス番号を通過しない(N、M)に(0,0)、見つけることができます。

DP(i)を設定する他の障害のパス番号を介さずに、i番目の障害物への(0,0)から表します。

i番目のバリア障壁のJ番目のいずれかのパスを上がることができるように、バリアを介してバーを転送するときに列挙jは、パスの数を計算する(障害物が経過した)、j番目の障壁を通して障壁を表します。

スルーi番目の障害物ので(0,0)から得られた他の障害、すなわち、DP(J)、および他の障害を通過できJ番目の障害物に繰り返し、(0,0)からの時間を避けるためにパスの数。
[値パスの総数を減算し、それはDP(i)はこのような答えを得ました。時間複雑\(^ O 2(K) \)

コード:

#include <stdio.h>
#define ll long long
int md=1000000007;
ll ny[200010],jn[200010],jc[200010];
int x[3010],y[3010],k=0;
int C(int n,int m)
{
    ll jg=(jc[n]*jn[m])%md;
    jg=(jg*jn[n-m])%md;
    return int(jg);
}
int dp[3010];
bool js[3010];
void ycl()
{
    ny[1]=jc[0]=jc[1]=jn[0]=jn[1]=1;
    for(int i=2;i<=200000;i++)
    {
        ny[i]=((md-md/i)*ny[md%i])%md;
        jc[i]=(i*jc[i-1])%md;
        jn[i]=(jn[i-1]*ny[i])%md;
    }
}
void dfs(int u)
{
    js[u]=true;
    dp[u]=0;
    for(int i=0;i<k;i++)
    {
        if(i!=u&&x[i]<=x[u]&&y[i]<=y[u])
        {
            if(!js[i])
                dfs(i);
            int t=((ll)dp[i]*C(x[u]-x[i]+y[u]-y[i],x[u]-x[i]))%md;
            dp[u]=(dp[u]+t)%md;
        }
    }
    dp[u]=(C(x[u]+y[u],x[u])-dp[u]+md)%md;
}
int main()
{
    freopen("path.in","r",stdin);
    freopen("path.out","w",stdout);
    ycl();
    int n,m,s;
    scanf("%d%d%d",&n,&m,&s);
    for(int i=0;i<s;i++)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        if((a==0&&b==0)||(a==n&&b==m))
        {
            printf("0");
            return 0;
        }
        bool zd=false;
        for(int j=0;j<k;j++)
        {
            if(x[j]==a&&y[j]==b)
            {
                zd=true;
                break;
            }
        }
        if(!zd)
        {
            x[k]=a;
            y[k]=b;
            k+=1;
        }
    }
    x[k]=n;
    y[k]=m;
    dfs(k);
    printf("%d",dp[k]);
    fclose(stdin);
    fclose(stdout);
    return 0;
}

おすすめ

転載: www.cnblogs.com/lnzwz/p/11246789.html