【动态规划】XMU 1588 01序列计数

题目链接：

　　http://acm.xmu.edu.cn/JudgeOnline/problem.php?id=1588

题目大意：

　　给n1个0和n2个1，连续的0不超过k1个，连续的1不超过k2个。问序列有几种（对1000000001(10⁹+1)取模）

　　n1,n2<=100 k1,k2<=10

题目思路：

　　【动态规划】

　　动态规划比较明显。

　　f[i][j][k][0]表示 i个0 j个1 末尾连续k个0的方案数

　　f[i][j][k][1]表示 i个0 j个1 末尾连续k个1的方案数

　　根据最后末尾取0还是取1可以推出状态转移方程。

　　初始化就是全为0或1的情况。

　　以下是f[][][][0]的状态转移方程。

　　当k>1时，f[i][j][k][0]=f[i-1][j][k-1][0]

　　当k=1时，f[i][j][k][0]=Σ(f[i-1][j][p][1])，p=0..k2;

　　

 //
 //by coolxxx
 //
 #include<iostream>
 #include<algorithm>
 #include<string>
 #include<iomanip>
 #include<memory.h>
 #include<time.h>
 #include<stdio.h>
 #include<stdlib.h>
 #include<string.h>
 #include<stdbool.h>
 #include<math.h>
 #define min(a,b) ((a)<(b)?(a):(b))
 #define max(a,b) ((a)>(b)?(a):(b))
 #define abs(a) ((a)>0?(a):(-(a)))
 #define lowbit(a) (a&(-a))
 #define sqr(a) ((a)*(a))
 #define swap(a,b) ((a)^=(b),(b)^=(a),(a)^=(b))
 #define eps 1e-8
 #define J 10
 #define MAX 0x7f7f7f7f
 #define PI 3.1415926535897
 #define mod 1000000001
 #define N 104
 #define M 14
 using namespace std;
 int n,m,lll,ans,cas;
 int nn,mm;
 int f[N][N][M][];
 int main()
 {
     #ifndef ONLINE_JUDGE
 //    freopen("1.txt","r",stdin);
 //    freopen("2.txt","w",stdout);
     #endif
     int i,j,k;
 //    while(~scanf("%s",s1))
     while(~scanf("%d",&n))
     {
         scanf("%d%d%d",&m,&nn,&mm);
         for(i=;i<=nn;i++)f[i][][i][]=;
         for(i=;i<=mm;i++)f[][i][i][]=;
         for(i=;i<=n;i++)
         {
             for(j=;j<=m;j++)
             {
                 f[i][j][][]=;
                 for(k=;k<=mm;k++)
                     f[i][j][][]=(f[i][j][][]+f[i-][j][k][])%mod;
                 for(k=;k<=nn;k++)
                     f[i][j][k][]=max(f[i-][j][k-][],f[i][j][k][]);

                 f[i][j][][]=;
                 for(k=;k<=nn;k++)
                     f[i][j][][]=(f[i][j][][]+f[i][j-][k][])%mod;
                 for(k=;k<=mm;k++)
                     f[i][j][k][]=max(f[i][j-][k-][],f[i][j][k][]);
             }
         }
         ans=;
         for(i=;i<=nn;i++)ans=(ans+f[n][m][i][])%mod;
         for(i=;i<=mm;i++)ans=(ans+f[n][m][i][])%mod;
         printf("%d\n",ans);
     }
     return ;
 }

 /*
 //

 //
 */