HDU 5459 Jesus Is Here (递推，组合数学)

有点麻烦的递推，递推的原则：向小的问题方向分解，注意边界。

字符串的递推式为

定义f为Si中的总方案数

首先可以得到

f_i=f_i-1+f_i-2+组合(s_i-2,s_i-1)

然后考虑Si-2和Si-1之间的组合

为了得到小的问题，进行拆分

为了以后表示的方便和逻辑上的清晰，把Si~Si之间的组合总长度定义出来

因为这里的si-2和si-2的中间还有一段Si-3

所以其组合总长度就可以表示为

Ci表示Si中cff出现的次数，Li表示Si的长度

定义一个函数ccl

最后还剩下一个部分Si-2和Si-3

定义

至此，总方案已经可以表示出来

然后再从顶往下的逐步细化

考虑g(i)的计算

这里面也有个递推，首先可以得到

为了方便再定义一个函数

显然这个容易计算得多

剩下的组合也可以表示出来了

最后一部分f12(i)

#include<bits/stdc++.h>
using namespace std;

#define PB push_back
#define MP make_pair
#define fi first
#define se second

typedef long long ll;


const int N = 201314+5;
ll f[N], c[N], l[N], g[N];
const int Mod =  530600414;

inline ll f21(int i) { return f[i]-f[i-2]-f[i-1]; }
inline ll ccl(int i,int j,int k) { return c[i]*c[j]%Mod*l[k]; }
inline ll f12(int i) { return g[i-2]+f21(i-1)+ccl(i-3,i-2,i-2); }//用到i-3

//#define LOCAL
int main()
{
#ifdef LOCAL
    freopen("in.txt","r",stdin);
#endif
    l[1] = 1; l[2] = 2;
    l[3] = 3; l[4] = 5;
    c[3] = 1; c[4] = 1;
    g[3] = 3; g[4] = 5;
    for(int i = 5; i <= 201314; i++){
        l[i] = (l[i-1]+l[i-2])%Mod;
        c[i] = (c[i-1]+c[i-2])%Mod;
        f[i] = (f[i-1]+f[i-2]+g[i-2]+ccl(i-2,i-2,i-3)+f12(i-1)+Mod)%Mod;//用到i-4,所以从5开始
        g[i] = (g[i-1]+g[i-2]+ccl(i-1,i-1,i-2)+ccl(i-2,i-2,i-1)+f12(i)+f21(i)+ccl(i-2,i-1,i))%Mod;//用到f[i]和l[i],所以f[i]l[i]在前面
    }
    int T, kas = 0; scanf("%d",&T);
    while(T--){
        int n; scanf("%d",&n);
        printf("Case #%d: %I64d
",++kas,f[n]);
    }
    return 0;
}