前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布

HDU4832

作者头像
triplebee
发布2018-01-12 15:24:13
4000
发布2018-01-12 15:24:13
举报
文章被收录于专栏:计算机视觉与深度学习基础

由于水平和竖直相互独立,所以可以分开计数,最后再用组合数算一下,万年老坑long long

代码语言:javascript
复制
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int n,m,k,x,y;
int dpx[1005][1005],dpy[1005][1005],c[1005][1005];
long long sum1[1005],sum2[1005];//这里不用longlong在最后相乘过程中会爆
bool judge(int t,int n)
{
    if(t>=1 && t<=n)
    return true;
    return false;
}
void init()
{
    memset(c,0,sizeof(c));
    for(int i=0;i<1005;i++)
    {
        c[i][0]=1;
        c[i][i]=1;
    }
    for(int i=1;i<1005;i++)
    {
        for(int j=1;j<i;j++)
        {
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%9999991;
        }
    }
}
int main()
{
    init();
    int t;
    cin>>t;
    for(int te=1;te<=t;te++)
    {
        cin>>n>>m>>k>>x>>y;
        memset(dpx,0,sizeof(dpx));
        memset(dpy,0,sizeof(dpy));
        memset(sum1,0,sizeof(sum1));
        memset(sum2,0,sizeof(sum2));
        dpx[0][x]=1;
        dpy[0][y]=1;
        for(int i=1;i<=k;i++)
        {
            for(int j=1;j<=n;j++)
            {
                for(int s=-2;s<=2;s++)
                {
                    if(s==0) continue;
                    int t=j+s;
                    if(!judge(t,n))continue;
                    dpx[i][t]=(dpx[i][t]+dpx[i-1][j])%9999991;
                }
            }
        }
        for(int i=1;i<=k;i++)
        {
            for(int j=1;j<=m;j++)
            {
                 for(int s=-2;s<=2;s++)
                 {
                     if(s==0) continue;
                     int t=j+s;
                     if(!judge(t,m))continue;
                     dpy[i][t]=(dpy[i][t]+dpy[i-1][j])%9999991;
                 }
            }
        }
        for(int i=0;i<=k;i++)
        {
            for(int j=1;j<=n;j++)
            {
                sum1[i]=(dpx[i][j]+sum1[i])%9999991;
                //cout<<dpx[i][j]<<" ";
            }
            //cout<<sum1[i]<<endl;
            //cout<<endl;
        }
        for(int i=0;i<=k;i++)
        {
            for(int j=1;j<=m;j++)
            {
                sum2[i]=(dpy[i][j]+sum2[i])%9999991;
            }
        }
        long long ans=0;
        for(int i=0;i<=k;i++)
        {
            long long t=c[k][i]*((sum1[i]*sum2[k-i])%9999991)%9999991;
            ans=(ans+t)%9999991;
        }
        cout<<"Case #"<<te<<":"<<endl;
        cout<<ans<<endl;
    }
    return 0;
}
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2014-08-04 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档