hdu 5576 dp

题目大意：给你一个长度为 n 的 字符串表示一个乘法，一次操作随机选两个字符进行交换，进行m次操作，让你求出所有可能操作

的答案和。  (1 <= n, m <= 50)

思路：巨难。。 对于固定位置的三个字符来说，把中间的看成*, 两边的为a，b， 这个产生的贡献为，a * b * C，C只与a, b和 * 的位置有关，

所以我们只要求出a*b的贡献就好啦。

我们定义dp[ a ][ b ][ c ][ k ] 表示三个字符分别是a，b，c， b在中间，还需要交换 k次的答案。

最后枚举三个位置求答案。

#include<bits/stdc++.h>
#define LL long long
#define fi first
#define se second
#define mk make_pair
#define pii pair<int, int>
#define y1 skldjfskldjg
#define y2 skldfjsklejg

using namespace std;

const int N =  + ;
const int M = 1e5 + ;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9 +;

char s[];
int n, m, a[], cnt[];
LL f[][][][], ten[];

void add(LL &a, LL b) {
    a += b; if(a >= mod) a -= mod;
}

LL dp(int a, int b, int c, int m) {
//    printf("%d %d %d %d $$\n", a, b, c, m);
    if(!m) return b ==  ? a * c : ;
    if(a > c) swap(a, c);
    LL &ans = f[a][b][c][m];
    if(~ans) return ans;
    ans = ;

    add(ans, dp(c, b, a, m - ));   // a <-> c
    add(ans, dp(b, a, c, m - ));   // a <-> b
    add(ans, dp(a, c, b, m - ));   // b <-> c

    LL num = (n - ) * (n - ) / ;

    if(num) add(ans, num * dp(a, b, c, m - ) % mod);  // a, b, c 之外的两个数相互交换

    int t[];
    memcpy(t, cnt, sizeof(t));
    t[a]--; t[b]--; t[c]--;
    for(int i = ; i <= ; i++) { // a, b, c 之外的一个数与 a 或 b 或 c 交换
        if(!t[i]) continue;
        add(ans, t[i] * dp(i, b, c, m - ) % mod); //挑一种数与a交换
        add(ans, t[i] * dp(a, i, c, m - ) % mod); //挑一种数与b交换
        add(ans, t[i] * dp(a, b, i, m - ) % mod); //挑一种数与c交换
    }

    return ans;
}

int main() {

    ten[] = ;
    for(int i = ; i <= ; i++) ten[i] = ten[i - ] *  % mod;

    int T; scanf("%d", &T);

    for(int cas = ; cas <= T; cas++) {
        memset(f, -, sizeof(f));
        memset(cnt, , sizeof(cnt));
        scanf("%d%s", &m, s);
        n = strlen(s);
        printf("Case #%d: ", cas);

        for(int i = ; i < n; i++) {
            a[i] = (s[i] == '*' ?  : s[i] - '');
            cnt[a[i]]++;
        }

        LL ans = ;

        for(int i = ; i < n; i++) {
            for(int j = i + ; j < n; j++) {
                for(int k = j + ; k < n; k++) {
                    add(ans, ten[j - i - ] * ten[n - k - ] % mod * dp(a[i], a[j], a[k], m) % mod);
                }
            }
        }

        printf("%lld\n", ans);
    }
    return ;
}

/*
2
1
1*2
2
1*2
*/