POJ1509 Glass Beads(最小表示法 后缀自动机)
2018-06-29 06:14:38来源:博客园 阅读 ()
Time Limit: 3000MS | Memory Limit: 10000K | |
Total Submissions: 4901 | Accepted: 2765 |
Description
The necklace should be made of glass beads of different sizes connected to each other but without any thread running through the beads, so that means the beads can be disconnected at any point. The actress chose the succession of beads she wants to have and the IBM promised to make the necklace. But then he realized a problem. The joint between two neighbouring beads is not very robust so it is possible that the necklace will get torn by its own weight. The situation becomes even worse when the necklace is disjoined. Moreover, the point of disconnection is very important. If there are small beads at the beginning, the possibility of tearing is much higher than if there were large beads. IBM wants to test the robustness of a necklace so he needs a program that will be able to determine the worst possible point of disjoining the beads.
The description of the necklace is a string A = a1a2 ... am specifying sizes of the particular beads, where the last character am is considered to precede character a1 in circular fashion.
The disjoint point i is said to be worse than the disjoint point j if and only if the string aiai+1 ... ana1 ... ai-1 is lexicografically smaller than the string ajaj+1 ... ana1 ... aj-1. String a1a2 ... an is lexicografically smaller than the string b1b2 ... bn if and only if there exists an integer i, i <= n, so that aj=bj, for each j, 1 <= j < i and ai < bi
Input
Output
Sample Input
4 helloworld amandamanda dontcallmebfu aaabaaa
Sample Output
10 11 6 5
Source
题目大意:对于给定的字符串,输出其最小表示法的第一个字符在第几位
最小表示法:最小表示法又叫做最小循环表示。
你可以直观的理解为对于一个字符串,选一个位置把它劈开,把前一半接到后一半,形成一个新的字符串,在这些新的字符串中字典序最小的即为字符串的最小表示。
最小表示法有专门的算法(三指针法?https://www.cnblogs.com/XGHeaven/p/4009210.html)
但是它可以轻松的被SAM解决
我们先把SAM建出来,然后从根节点开始,每次走最小的转移边,走$|S|$次,所得的串即为最小表示
那么第一个字母可用通过$len-|S|$找到
#include<cstdio> #include<cstring> using namespace std; const int MAXN = 40001; int N; char s[MAXN]; int fa[MAXN], len[MAXN], ch[MAXN][26], tot = 1, last = 1, root = 1; void insert(int x) { int now = ++tot, pre = last; last = now; len[now] = len[pre] + 1; for(; pre && !ch[pre][x]; pre = fa[pre]) ch[pre][x] = now; if(!pre) fa[now] = root; else { int q = ch[pre][x]; if(len[q] == len[pre] + 1) fa[now] = q; else { int nows = ++tot; len[nows] = len[pre] + 1; memcpy(ch[nows], ch[q], sizeof(ch[q])); fa[nows] = fa[q]; fa[q] = fa[now] = nows; for(; pre && ch[pre][x] == q; pre = fa[pre]) ch[pre][x] = nows; } } } int main() { #ifdef WIN32 freopen("a.in", "r", stdin); #endif int QwQ; scanf("%d", &QwQ); while(QwQ--) { memset(fa, 0, sizeof(fa)); memset(ch, 0, sizeof(ch)); memset(len , 0, sizeof(len)); tot = last = root = 1; scanf("%s", s + 1); N = strlen(s + 1); for(int i = 1; i <= N; i++) s[i + N] = s[i]; N <<= 1; for(int i = 1; i <= N; i++) insert(s[i] - 'a'); int now = root, tot = 0; for(; tot <= N / 2; tot++) { for(int i = 0; i <= 25; i++) if(ch[now][i]) {now = ch[now][i]; putchar(i + 'a'); break;} } printf("%d\n", len[now] - N / 2); } return 0; }
标签:
版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有
- [AtCoder][ARC082]Sandglass 题解 2018-08-02
- POJ1286 Necklace of Beads(Polya定理) 2018-07-12
- 【洛谷P3498】 [POI2010]KOR-Beads 2018-06-17
IDC资讯: 主机资讯 注册资讯 托管资讯 vps资讯 网站建设
网站运营: 建站经验 策划盈利 搜索优化 网站推广 免费资源
网络编程: Asp.Net编程 Asp编程 Php编程 Xml编程 Access Mssql Mysql 其它
服务器技术: Web服务器 Ftp服务器 Mail服务器 Dns服务器 安全防护
软件技巧: 其它软件 Word Excel Powerpoint Ghost Vista QQ空间 QQ FlashGet 迅雷
网页制作: FrontPages Dreamweaver Javascript css photoshop fireworks Flash