• 博客园logo
  • 会员
  • 众包
  • 新闻
  • 博问
  • 闪存
  • 赞助商
  • HarmonyOS
  • Chat2DB
    • 搜索
      所有博客
    • 搜索
      当前博客
  • 写随笔 我的博客 短消息 简洁模式
    用户头像
    我的博客 我的园子 账号设置 会员中心 简洁模式 ... 退出登录
    注册 登录
XD-TEST
博客园    首页    新随笔    联系   管理    订阅  订阅

bzoj1925 [Sdoi2010]地精部落

  dp,f[i]表示放置了1到i的方案数,g[i]表示放置了1到i且最左侧为山谷的方案数,可以得到dp方程f[i]=Σg[j]*g[i-1-j]*c[i-1][j](0<=j<i),求组合数直接开数组会爆,可以用滚动数组或者开个vector计算,不会爆。

     代码

 1 #include<cstdio>
 2 #include<vector>
 3 using namespace std;
 4 const int N  = 4210;
 5 int n,P,i,j;
 6 vector<int> c[N];
 7 long long f[N],g[N];
 8 int main()
 9 {
10     scanf("%d%d",&n,&P);
11     for (i=0;i<=n;i++)
12     {
13         c[i].push_back(1);
14         for (j=1;j<=i;j++)
15         {
16             int tmp=0;
17             tmp=c[i-1][j-1];
18             if (j<=i-1) tmp=(tmp+c[i-1][j])%P;
19             c[i].push_back(tmp);
20         }
21     }
22     g[1]=g[0]=1;
23     for (i=2;i<=n;i++)
24     for (j=0;j<i;j++)
25     {
26         (f[i]+=(g[j]*g[i-1-j])%P*c[i-1][j])%=P;
27         if (j&1)
28         (g[i]+=(g[j]*g[i-1-j])%P*c[i-1][j])%=P;
29     }
30     printf("%lld",f[n]);
31 }

 

posted @ 2016-05-17 20:44  fzmh  阅读(422)  评论(0)    收藏  举报
刷新页面返回顶部
博客园  ©  2004-2025
浙公网安备 33010602011771号 浙ICP备2021040463号-3