POJ 1204 Word Puzzles AC自动机题解

AC自动机的灵活运用,本题关键是灵活二字。

因为数据不是很大,时间要求也不高的缘故,所以本题有人使用暴力法也过了,有人使用Trie,然后枚举八个方向,也过了,不过速度很慢。

当然有人使用AC自动机没AC的,在讨论区里喊AC自动机超时的,那是因为不会灵活运用,或者是硬套模板的,AC了速度也不会快。


给出本人的算法思路:

1 把需要查找的关键字建立Trie, 然后构造AC自动机

2 查找的时候分八个方向查找,比如棋盘是board[N][M],那么就可以循环i(0->N-1),然后每次把board[i]当做一个文本,做过HDU的key word那道题目的人都知道如何搜索一个文本中有多少关键字出现,这样不就把本题转换为查找文本关键字的问题了。 不过本题是需要八个方向的八个文本都搜索一遍,八是个不大的常数。其他方向如何构造一个文本,就可以自己思考一下,不难。

3 记录结果就看个人的编程能力了,这里的技巧是逆转关键字建立Trie,那么搜索到尾就是每个关键字的开头了。


AC自动机虽然高级,但是也还是一个算法工具,套模板,一大抄的代码没什么意义的,还是彻底消化掉,然后灵活运用才是关键。

对于本题写出高效的算法,那么是超过5星级难度的。


下面给出个上面思路的实现算法,速度很快的。

出处: http://blog.csdn.net/kenden23

#include <stdio.h>
#include <string.h>
#include <queue>
using std::queue;

const int MAX_LEN = 1001;
const int ARR_SIZE = 26;
char board[MAX_LEN][MAX_LEN];
char word[MAX_LEN], gDir[MAX_LEN];
int gIndices[MAX_LEN][2];
int L, C, W;

inline int getIndex(char ch) { return ch-'A'; }

struct Node
{
	int n, index;
	Node *fail;
	Node *arr[ARR_SIZE];
};

void clearNode(Node *p)
{
	p->n = 0;
	p->index = 0;
	p->fail = NULL;
	for (int i = 0; i < ARR_SIZE; i++)
	{
		p->arr[i] = NULL;
	}
}

Node gPool[MAX_LEN*MAX_LEN];
int gPoolId;
Node *Trie;

void insertTrie(char w[], int index)
{
	Node *pCrawl = Trie;
	int len = strlen(w);
	for (int i = len-1; i >= 0; i--)//逆向插入字符串,方便找到起头节点
	{
		int id = getIndex(w[i]);
		if (!pCrawl->arr[id])
		{
			pCrawl->arr[id] = &gPool[gPoolId++];
			clearNode(pCrawl->arr[id]);
		}
		pCrawl = pCrawl->arr[id];
	}
	pCrawl->n++;
	pCrawl->index = index;
}

void buildFail()
{
	queue<Node *> qu;
	qu.push(Trie);
	while (!qu.empty())
	{
		Node *pCrawl = qu.front();
		qu.pop();
		for (int i = 0; i < ARR_SIZE; i++)
		{
			if (!pCrawl->arr[i]) continue;
			pCrawl->arr[i]->fail = Trie;
			Node *fail = pCrawl->fail;
			while (fail)
			{
				if (fail->arr[i])
				{
					pCrawl->arr[i]->fail = fail->arr[i];
					break;
				}
				fail = fail->fail;
			}
			qu.push(pCrawl->arr[i]);
		}
	}
}

void searchDirect(int r, int c, int rDir, int cDir, char dir)
{
	Node *pCrawl = Trie;
	while (0<=r && 0<=c && r < L && c < C)
	{
		int id = getIndex(board[r][c]);		
		while (!pCrawl->arr[id] && pCrawl != Trie) pCrawl = pCrawl->fail;
		if (pCrawl->arr[id])
		{
			pCrawl = pCrawl->arr[id];
			Node *tmp = pCrawl;
			while (tmp && tmp->n != -1)
			{
				if (tmp->n)//record results
				{
					gIndices[tmp->index][0] = r;
					gIndices[tmp->index][1] = c;
					gDir[tmp->index] = dir;
				}
				tmp->n = -1;
				tmp = tmp->fail;
			}
		}
		r += rDir, c += cDir;	//Watch out, if we use "continue"!
	}
}

void searchAll()
{
	char dir = 'E', counterDir = 'A';
	for (int c = 0; c < C; c++)
	{
		searchDirect(L-1, c, -1, 0, dir);
		searchDirect(0, c, 1, 0, counterDir);
	}
	dir = 'G', counterDir = 'C';
	for (int r = 0; r < L; r++)
	{
		searchDirect(r, 0, 0, 1, dir);
		searchDirect(r, C-1, 0, -1, counterDir);
	}
	dir = 'H', counterDir = 'D';
	for (int c = 0; c < C; c++)
	{
		searchDirect(0, c, 1, 1, dir);
		searchDirect(L-1, c, -1, -1, counterDir);
	}
	for (int r = 1; r < L; r++)
	{
		searchDirect(r, 0, 1, 1, dir);
		searchDirect(r, C-1, -1, -1, counterDir);
	}
	dir = 'F', counterDir = 'B';
	for (int c = 0; c < C; c++)
	{
		searchDirect(L-1, c, -1, 1, dir);
		searchDirect(0, c, 1, -1, counterDir);
	}
	for (int r = 0; r < L; r++)
	{
		searchDirect(r, 0, -1, 1, dir);
		searchDirect(r, C-1, 1, -1, counterDir);
	}
}

int main()
{
	Trie = &gPool[0];
	while (scanf("%d %d %d", &L, &C, &W) != EOF)
	{
		clearNode(Trie);
		gPoolId = 1;
		getchar();
		for (int i = 0; i < L; i++)
		{
			gets(board[i]);
		}
		for (int i = 0; i < W; i++)
		{
			gets(word);
			insertTrie(word, i);
		}
		buildFail();
		searchAll();
		for (int i = 0; i < W; i++)
		{
			printf("%d %d %c\n", gIndices[i][0], gIndices[i][1], gDir[i]);
		}
	}
	return 0;
}



  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
POJ3635是一道经典的数学题,需要使用一些数学知识和算法进行解决。 题目描述: 给定四个正整数 a、b、p 和 k,求 a^b^p mod k 的值。 解题思路: 首先,我们可以将指数 b^p 写成二进制形式:b^p = c0 * 2^0 + c1 * 2^1 + c2 * 2^2 + ... + ck * 2^k,其中 ci 为二进制数的第 i 位。 然后,我们可以通过快速幂算法来计算 a^(2^i) mod k 的值。具体来说,我们可以用一个变量 x 来存储 a^(2^i) mod k 的值,然后每次将 i 加 1,如果 ci 为 1,则将 x 乘上 a^(2^i) mod k,最后得到 a^b^p mod k 的值。 代码实现: 以下是 Java 的代码实现: import java.util.*; import java.math.*; public class Main { public static void main(String[] args) { Scanner sc = new Scanner(System.in); BigInteger a = sc.nextBigInteger(); BigInteger b = sc.nextBigInteger(); BigInteger p = sc.nextBigInteger(); BigInteger k = sc.nextBigInteger(); BigInteger ans = BigInteger.ONE; for (int i = 0; i < p.bitLength(); i++) { if (b.testBit(i)) { ans = ans.multiply(a.modPow(BigInteger.ONE.shiftLeft(i), k)).mod(k); } } System.out.println(ans); } } 其中,bitLength() 函数用于获取二进制数的位数,testBit() 函数用于判断二进制数的第 i 位是否为 1,modPow() 函数用于计算 a^(2^i) mod k 的值,multiply() 函数用于计算两个 BigInteger 对象的乘积,mod() 函数用于计算模数。 时间复杂度: 快速幂算法的时间复杂度为 O(log b^p),其中 b^p 为指数。由于 b^p 的位数不超过 32,因此时间复杂度为 O(log 32) = O(1)。 总结: POJ3635 是一道经典的数学题,需要使用快速幂算法来求解。在实现时,需要注意 BigInteger 类的使用方法,以及快速幂算法的细节。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值