Gym
来源:互联网 发布:java log4j2.xml 编辑:程序博客网 时间:2024/06/18 12:25
题目:给出A,B两个字符串,求A,B中相同的回文串的对数
思路:对A,B分别建一个回文树,然后分别从奇数根节点和偶数做一次dfs
代码:
#pragma comment(linker, "/STACK:1024000000,1024000000")#include<iostream>#include<algorithm>#include<ctime>#include<cstdio>#include<cmath>#include<cstring>#include<string>#include<vector>#include<map>#include<set>#include<queue>#include<stack>#include<list>#include<numeric>using namespace std;#define LL long long#define ULL unsigned long long#define INF 0x3f3f3f3f3f3f3f3f#define mm(a,b) memset(a,b,sizeof(a))#define PP puts("*********************");template<class T> T f_abs(T a){ return a > 0 ? a : -a; }template<class T> T gcd(T a, T b){ return b ? gcd(b, a%b) : a; }template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}// 0x3f3f3f3f3f3f3f3fconst int MAXN = 200050, SIZE = 26;struct Palindromic_Tree {int next[MAXN][SIZE];//next指针int fail[MAXN];//fail指针int cnt[MAXN];//表示节点i表示的本质不同的串的个数(建树时求出的不是完全的,最后Count()函数跑一遍以后才是正确的)int num[MAXN];//表示以节点i表示的最长回文串的最右端点为回文串结尾的回文串个数int len[MAXN];//len[i]表示节点i表示的回文串的长度int S[MAXN];//存放添加的字符 S[MAXN << 1]int last;//指向上一个字符所在的节点,方便下一次addint n;//字符数组指针 //lbd, rbdint p;//节点指针/*int lbd,rbd;int last[2];//0表示在字符串的首部加上一个字符,1表示在尾部加上一个字符*/int NewNode(int l) {//新建节点for(int i = 0; i < SIZE; ++i)next[p][i] = 0;cnt[p] = num[p] = 0;len[p] = l;return p++;}void Init(/*int x*/) {//初始化p = n = 0;NewNode(0); NewNode(-1);last = 0;S[n] = -1;//开头放一个字符集中没有的字符,减少特判/*p=0;NewNode(0); NewNode(-1);last[0] = last[1] = 0;memset(S, -1, sizeof(S));lbd = x; rbd = x - 1;//x是字符串的最大长度*/fail[0] = 1;}int GetFail(int x/*, int d*/) {//和KMP一样,失配后找一个尽量最长的while(S[n - len[x] - 1] != S[n]) x = fail[x];/*if(d == 0) while(S[lbd + len[x] + 1] != S[lbd]) x = fail[x];else while(S[rbd - len[x] - 1] != S[rbd]) x = fail[x];*/return x;}void Add(int c/*, int d*/) {//d=0表示在首部添加字符,d=1表示在尾部添加字符S[++n] = c;/*if(d == 0) S[--lbd] = c; else S[++rbd] = c;int cur = GetFail(last[d], d);*/int cur = GetFail(last);//通过上一个回文串找这个回文串的匹配位置if(!next[cur][c]) {//如果这个回文串没有出现过,说明出现了一个新的本质不同的回文串int now = NewNode(len[cur] + 2);//新建节点fail[now] = next[GetFail(fail[cur]/*, d*/)][c];//和AC自动机一样建立fail指针,以便失配后跳转next[cur][c] = now;num[now] = num[fail[now]] + 1;}last = next[cur][c];cnt[last]++; /*last[d]=next[cur][c];if(len[last[d]] == rbd - lbd + 1) last[d ^ 1] = last[d];return num[last[d]];*/}void Count() {for(int i = p - 1; i >= 0; --i)cnt[fail[i]] += cnt[i];//父亲累加儿子的cnt,因为如果fail[v]=u,则u一定是v的子回文串!}}treeA,treeB;LL dfs(int posA,int posB){ LL sum=0; for(int i=0;i<SIZE;i++) if(treeA.next[posA][i]!=0&&treeB.next[posB][i]!=0){ int x=treeA.next[posA][i]; int y=treeB.next[posB][i]; sum+=(LL)treeA.cnt[x]*treeB.cnt[y]+dfs(x,y); } return sum;}char str[MAXN];int main(){ int T,cas=0; scanf("%d",&T); while(T--){ treeA.Init(); treeB.Init(); scanf("%s",str); for(int i=0;str[i]!='\0';i++) treeA.Add(str[i]-'a'); treeA.Count(); scanf("%s",str); for(int i=0;str[i]!='\0';i++) treeB.Add(str[i]-'a'); treeB.Count(); LL ans=0; ans+=dfs(0,0); ans+=dfs(1,1); printf("Case #%d: %lld\n",++cas,ans); } return 0;}
阅读全文
0 0
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- Gym
- maven工程引入jstl,且出现uri异常的解决
- 在命令行中使用wsimport命令在MyEclipse中生成webservice文件,无法生成,出现directory not found这种情况
- java获取路径问题
- Android 开源项目分类汇总
- 将博客搬至简书,两边都会同时更新
- Gym
- StringBuffer
- 利用SSH工具上传文件到CentOS 6.5
- Oracle通过Navicat建表查询时表不存在的问题解决
- Java中基础数据类型的面试相关+String类
- 探索 ConcurrentHashMap 高并发性的实现机制
- javaweb中servlet生命周期介绍
- GBDT与梯度的理解
- 在maven中启动tomcat报错java.lang.ClassNotFoundException: org.springframework.web.context.ContextLoad