Ring 【HDU - 2296】【AC自动机+DP】

题目链接


  问的是可以包含情况下的最大的匹配,就是譬如说"vve"可以覆盖"v"、"ve"的值,所以是值的叠加,求最后的最大值的最小长度,如果长度相等就是去取最小字典序。

  思路:可还真是捣腾的难呀,弄了一下午终于dubug出来的了,没想到竟然是建立fail指针的时候出了少许错误,哎……那么,怎么去解这道题?毕竟要输出答案的字符串,所以我还去补了下下"strcpy()"、"strcmp()"这些C函数的用法,然后就是解题思路了。

  可以去推我们所要建立的长度为N的字符串的DP,那么就是首先利用字典树Trie来得到每个状态的对应的位置,然后去推之后的状态,有相等、大于的关系,分别去判断即可。

  总之,还是道AC自动机+DP的入门好题了,有很多的想法。


#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
#define MP(a, b) make_pair(a, b)
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
int N, M, tot, dp[55][1005];
char virtus[105][15];
char str[55][1005][55];
struct node
{
    int nex[26], val, fail;
    node() { memset(nex, 0, sizeof(nex)); fail = val = 0; }
    void clear() { memset(nex, 0, sizeof(nex)); fail = val = 0; }
}a[1005];
void insert(int ti, int val)
{
    int len = (int)strlen(virtus[ti]), temp = 0;
    for(int i=0; i Q;
    Q.push(0);
    int temp, p, son;
    while(!Q.empty())
    {
        temp = Q.front();   Q.pop();
        for(int i=0; i<26; i++)
        {
            son = a[temp].nex[i];
            if(son)
            {
                if(!temp) a[son].fail = 0;
                else
                {
                    p = a[temp].fail;
                    while(p && !a[p].nex[i]) p = a[p].fail;
                    a[son].fail = a[p].nex[i];
                }
                if(a[a[son].fail].val) a[son].val += a[a[son].fail].val;
                Q.push(son);
            }
            else a[temp].nex[i] = a[a[temp].fail].nex[i];
        }
    }
}
void solve()
{
    dp[0][0] = 0;
    char s[55], ans[55], temp;
    int son;
    strcpy(str[0][0], "\0");
    strcpy(ans, "\0");
    for(int i=1; i<=N; i++)
    {
        for(int j=0; j<=tot; j++)
        {
            if(dp[i-1][j] == -1) continue;
            strcpy(s, str[i-1][j]);
            for(int k=0; k<26; k++)
            {
                son = a[j].nex[k];
                temp = 'a' + k;
                s[i-1] = temp;
                s[i] = '\0';
                if(dp[i][son] < dp[i-1][j] + a[son].val)
                {
                    dp[i][son] = dp[i-1][j] + a[son].val;
                    strcpy(str[i][son], s);
                }
                else if(dp[i][son] == dp[i-1][j] + a[son].val && strcmp(str[i][son], s) > 0) strcpy(str[i][son], s);
            }
        }
    }
    int maxx = 0;
    for(int i=1; i<=N; i++)
    {
        for(int j=0; j<=tot; j++)
        {
            if(dp[i][j] > maxx)
            {
                maxx = dp[i][j];
                strcpy(ans, str[i][j]);
            }
            else if(dp[i][j] == maxx && strlen(ans) >= strlen(str[i][j]) && strcmp(ans, str[i][j]) > 0) strcpy(ans, str[i][j]);
        }
    }
    printf("%s\n", ans);
}
inline void init()
{
    tot = 0;
    a[0].clear();
    memset(dp, -1, sizeof(dp));
}
int main()
{
    int T;  scanf("%d", &T);
    while(T--)
    {
        init();
        scanf("%d%d", &N, &M);
        for(int i=1; i<=M; i++) scanf("%s", virtus[i]);
        for(int i=1; i<=M; i++)
        {
            int e1; scanf("%d", &e1);
            insert(i, e1);
        }
        build_fail();
        solve();
    }
    return 0;
}
/*
1
4 2
love ve
6 5
答案:love
 
1
49 3
loveufoever
ever
al
12 7 3
答案:aloveufoeveraloveufoeveraloveufoeveraloveufoever
*/

 

你可能感兴趣的:(DP动态规划,数据结构,AC自动机)