Password Suspects LA4126

78 篇文章 2 订阅
55 篇文章 0 订阅


AC自动机+状态压缩递推,通过这题主要学习了一下输出答案的方法,DP的话输出状态比较好办,因为最优决策一般只是一条路径,所以记录一下前驱可以用线性的复杂度输出路径,而递推确是多条路径的合集,输出路径比较麻烦,需要根据最终状态向前递推,并把经过的路径节点全部标记出来,由于递推状态设计本身是不会有重复的路径的,所以最后再从头dfs一遍访问过的节点并且记录,当递归深度合适时,输出记录即可。



#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <queue>
#include <algorithm>
#include <vector>
#include <cstring>
#include <stack>
#include <cctype>
#include <utility>   
#include <map>
#include <string>  
#include <climits> 
#include <set>
#include <string>    
#include <sstream>
#include <utility>   
#include <ctime>

using std::priority_queue;
using std::vector;
using std::swap;
using std::stack;
using std::sort;
using std::max;
using std::min;
using std::pair;
using std::map;
using std::string;
using std::cin;
using std::cout;
using std::set;
using std::queue;
using std::string;
using std::istringstream;
using std::make_pair;
using std::getline;
using std::greater;
using std::endl;
using std::multimap;

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int, int> PAIR;
typedef multimap<int, int> MMAP;

const int MAXN(110);
const int SIGMA_SIZE(26);
const int MAXM(110);
const int MAXE(4000010);
const int MAXH(18);
const int INFI((INT_MAX-1) >> 2);
const int MOD(10007);
const ULL LIM(1000000000000000ull);

int point[12];

int que[MAXN];
int front, back;

struct AC
{
    int ch[MAXN][SIGMA_SIZE];
    int f[MAXN];
    int val[MAXN];
    int size;
    
    void init()
    {
        memset(ch[0], 0, sizeof(ch[0]));
        f[0] = val[0] = 0;
        size = 1;
    }
    inline int idx(char temp)
    {
        return temp-'a';        
    }
    
    void insert(char *S, int tv)
    {
        int u = 0, id;
        for(; *S; ++S)
        {
            id = idx(*S);
            if(!ch[u][id])
            {
                memset(ch[size], 0, sizeof(ch[size]));
                val[size] = 0;
                ch[u][id] = size++;
            }
            u = ch[u][id];
        }
		val[u] |= (1 << tv);
    }  
    void construct()
    {
        front = back = 0;
        int cur, u;
        for(int i = 0; i < SIGMA_SIZE; ++i)
        {
            u = ch[0][i];
            if(u)
            {
                que[back++] = u;
                f[u] = 0;
            }
        }
        while(front < back)
        {
            cur = que[front++];
            for(int i = 0; i < SIGMA_SIZE; ++i)
            {
                u = ch[cur][i];
                if(u)
                {
                    que[back++] = u;
                    f[u] = ch[f[cur]][i];
                    val[u] |= val[f[u]];
                }
                else
                    ch[cur][i] = ch[f[cur]][i];
            }
        }
    }
};

AC ac;

LL table[26][MAXN][1 << 10];
bool vis[26][MAXN][1 << 10];
char rec[27];
int n;

void dfs(int dep, int cur, int state)
{
	if(dep == n)
	{
		rec[dep] = '\0';
		printf("%s\n", rec);
		return;
	}
	for(int i = 0; i < SIGMA_SIZE; ++i)
		if(vis[dep+1][ac.ch[cur][i]][state|ac.val[ac.ch[cur][i]]])
		{
			rec[dep] = 'a'+i;
			dfs(dep+1, ac.ch[cur][i], state|ac.val[ac.ch[cur][i]]);
		}
}

void solve(int m)
{
	int lim = (1 << m)-1, temp;
	LL tv;
	for(int i = 0; i <= n; ++i)
		for(int j = 0; j < ac.size; ++j)
			for(int k = 0; k <= lim; ++k)
				table[i][j][k] = 0;
	table[0][0][0] = 1;
	for(int i = 0; i < n; ++i)
	{
		for(int j = 0; j < ac.size; ++j)
			for(int k = 0; k <= lim; ++k)
				if(table[i][j][k])
				{
					tv = table[i][j][k];
					for(int k2 = 0; k2 < SIGMA_SIZE; ++k2)
						table[i+1][ac.ch[j][k2]][k|ac.val[ac.ch[j][k2]]] += tv;
				}
	}
	LL ans = 0;
	for(int i = 0; i < ac.size; ++i)
		ans += table[n][i][lim];
	printf(" %lld suspects\n", ans);
	if(ans <= 42)
	{
		for(int i = 0; i <= n; ++i)
			for(int j = 0; j < ac.size; ++j)
				for(int k = 0; k <= lim; ++k)
					vis[i][j][k] = false;
		for(int i = 0; i < ac.size; ++i)
			if(table[n][i][lim])
				vis[n][i][lim] = true;
		for(int i = n-1; i >= 0; --i)
			for(int j = 0; j < ac.size; ++j)
				for(int k = 0; k <= lim; ++k)
					if(table[i][j][k])
					{
						for(int k2 = 0; k2 < SIGMA_SIZE; ++k2)
						{
							temp = ac.ch[j][k2];
							if(vis[i+1][temp][k|ac.val[temp]])
								vis[i][j][k] = true;
						}
					}
		dfs(0, 0, 0);
	}
}

char str[15];

int main()
{
	int m, n_case(0);
	while(scanf("%d%d", &n, &m), n+m)
	{
		ac.init();
		for(int i = 0; i < m; ++i)
		{
			scanf("%s", str);
			ac.insert(str, i);
		}
		ac.construct();
		printf("Case %d:", ++n_case);
		solve(m);
	}
	return 0;	
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值