[BZOJ2160]拉拉队排练(回文树)

原创 2016年05月30日 21:44:23

题意:求一个字符串的长度前k的奇数长度回文串的长度的乘积模一个数。

和APIO那个回文串的题差不多,可以先manacher搞出所有本质不同的回文串,再用后缀数组查查出现了多少次即可。这里偷懒写的回文树,注意跑完了之后那个cnt需要沿着后缀链接向上传递一下。

注意那个K是10^12不是10^9。。。开的int丢了5分,还好出题人比较良心只有一个点的K超了int。

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define erp(i,a,b) for(int i=a;i>=b;--i)
#define LL long long
using namespace std;
const int mo = 19930726;
const int MAXN = 1000005;
int N;
LL sum[MAXN], K;

int ksm(int a, int b)
{
	int res = 1;
	for (; b; b>>=1, a=1ll*a*a%mo)
		if(b&1) res = 1ll*res*a%mo;
	return res;
}

int ch[MAXN][26], fa[MAXN];
int cnt[MAXN], len[MAXN];
int last, n, ncnt;
char s[MAXN];
inline int NewNode(int l)
{
	return len[ncnt] = l, ncnt++;
}
void initpam()
{
	NewNode(0), NewNode(-1);
	s[0] = -1, fa[0] = 1;
}
inline int getfail(int x)
{
	while (s[n-len[x]-1] != s[n]) x = fa[x];
	return x;
}

void ins(int c)
{
	s[++n] = c;
	int cur = getfail(last);
	if (!ch[cur][c])
	{
		int now = NewNode(len[cur]+2);
		fa[now] = ch[getfail(fa[cur])][c];
		ch[cur][c] = now;
	}
	last = ch[cur][c];
	cnt[last] ++;
}

void calc()
{
	erp(i, ncnt-1, 0) cnt[fa[i]] += cnt[i];
	for (int i = 0; i<ncnt; ++i)
		if (len[i]>0&&(len[i]&1)) sum[len[i]] += cnt[i];
}

int main()
{
	freopen("rehearse.in","r",stdin);
	freopen("rehearse.out","w",stdout);
	char c;
	scanf("%d %lld\n", &N, &K);
	initpam();
	while (1)
	{
		c = getchar();
		if (c>='a'&&c<='z') ins(c-'a');
		else break;
	}
	calc();
	int ans = 1;
	for (int i = N; i; --i)
	{
		if (!(i&1)) continue;
		if (K<=sum[i]) { ans = 1ll*ans*ksm(i, K)%mo; K = 0; break; }
		else ans = 1ll*ans*ksm(i, sum[i])%mo, K -= sum[i];
	}
	if (K) puts("-1");
	else printf("%d\n", ans);
	return 0;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

Setup a simple host-to-host VPN with openvpn

http://www.eazynet.de/setup_a_simple_host_to_host_vpn_with_openvpnhttp://www.eazynet.de/setup_a_simp...

双连通缩点 POJ 1515

B - 双连通 Time Limit:1000MS    Memory Limit:10000KB    64bit IO Format:%I64d & %I64u SubmitStatus ...

poj 1515 Street Directions 【tarjan 求割边】【确定无向图边的方向 将其变成有向强连通图】

Street Directions Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 753 ...

NOIP2013华容道

NOIP2014华容道 说起来这道题还挺有难度的,我用了两个小时才把它AC,要是在赛场上的话。。。。这种题就果断放弃了, 下面步入正题 题目描述 Description 小 B 最近迷上了华容...

hdu3068之manacher算法+详解

最长回文 Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submi...

NOIP2015提高组第二题信息传递[图论]

这道题我一开始做感觉没有任何思路,但一开始就觉得是图论,就在思考用邻接矩阵还是链式前向星,可是邻接矩阵每判断一次就是$O(n^2)$,在次数最多是链式前向星空间复杂度达$O(n^2)$,直接爆空间,后...

10.9【noip动规专题考试】总结——众人AK,唯我爆零

(上不了p站我要死了)今天考试,看了前两道题之后就知道是dp专题了 然而却因为各种错误挂了。。。。有读题不仔细,也有基础不够扎实。题目总结 首先肯定是背包dp 发现买一次与买多次是不一样的效果。...

BZOJ 1008: [HNOI2008]越狱 计数原理

1008: [HNOI2008]越狱Time Limit: 1 Sec Memory Limit: 162 MB Submit: 10200 Solved: 4417 [Submit][Sta...

【bzoj4318】【OSU!】期望dp——维护多个期望值递推

[pixiv] https://www.pixiv.net/member_illust.php?mode=medium&illust_id=62369739 Description osu ...

【bzoj2160】拉拉队排练 manacher

先用manacher求出数组后,用一个差分数组计算每一个长度的字符串出现的次数,暴力计算就好了。 #include #include #include #include #include #...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)