【训练】2018 ICPC Regional Nanjing

待学知识点：

wqs二分

*Pollard Rho

字符串问题（KMP、exKMP、*PAM、SA、SAM）

待复习知识点：

KMP

网络流

A.Adrien and Austin

签到。容易发现先手只要能把石头分成相同数量的两组就必胜，因为接下来可以复制对方操作。

B.Tournament

Unsolved

C.Cherry and Chocolate

Unsolved

D.Country Meow

最小球覆盖。

三分做法比较好写，但不会证明。

退火不懂。

E.Eva and Euro coins

Unsolved

F.Frank

Unsolved

G.Pyramid

通过差分应用得到这种纯找规律的东东如果差分x次后是等差数列，则其是最高次为x+1的多项式。

打表得到系数即可。

打表就把点扔进去n^3判断即可。

H.Huge Discount

Unsloved

I.Magic Potion

训练时没仔细想蒙了个两次二分图匹配，冲了一发就过了。

实际上也是一个最大流板子题，直接建1号源点和2号源点全部连向左半图，再从1号源点连一条权值为k的边到2号源点即可。

原理上两者好像差不多。

J.Prime Game

我的做法是先预处理质数，然后用质数去标记1e6以内的数，如果标记到序列中的数就把序列下标放入优先队列，然后稍微统计一下答案即可。

但是不太会算时间复杂度，类比埃氏筛的话应该是O(NloglogN)的。

直接暴力Pollard Rho应该也可以？

K.Kangaroo Puzzle

Unsolved

L.Lagrange the Chef

Unsolved

M.Mediocre String Problem

由于要求描述从T串中找出的某个串一定与S串找出的那个串部分前缀回文，因此只需要枚举S的每个点，然后二分+hash跟T匹配即可，然后就只需要知道每个点为左端点的回文串个数就行了，这个在manacher的时候差分一下就可以了。

顺带，本题其他做法所需算法（包括但不限于exKMP、SA(M)、PAM）我都不会。。。。。

有空先学一下SA和KMP吧，听上去好入门一些，顺便给各类自动机铺路。

Code:(注意二分写法，以后都用l<=r,l=mid+1,r=mid-1吧）

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<vector>
#include<map>
using namespace std;
const long long maxn=1e6+5;
const long long b1=13331;
const long long b2=13531;
const long long mod1=998244353;
const long long mod2=1e9+7;
inline long long read()
{
	long long x=0,f=1;char c=getchar();
	while(c<'0'||c>'9')
	{
		if(c=='-') f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9')
	{
		x=x*10+(c-'0');
		c=getchar();
	}
	return x*f;
}
char s[maxn],t[maxn],cp[maxn<<1];
long long len1,len2;
long long p[maxn<<1],cf[maxn],d[maxn],sum[maxn];
long long h1[maxn],h2[maxn],h3[maxn],h4[maxn];
long long jz1[maxn],jz2[maxn];
long long las;
inline void init()
{
	cp[0]='~';
	for(long long i=1;i<=len1;i++)
	{
		cp[i*2-1]='#';
		cp[i*2]=s[i];
	}
	cp[len1*2+1]='#';
}
inline void manacher()
{	
	long long pos=0,maxr=0;
	for(long long i=1;i<=len1*2+1;i++)
	{
		if(i<maxr) p[i]=min(p[pos*2-i],maxr-i);
		while(cp[i+p[i]+1]==cp[i-p[i]-1]) p[i]++;
		if(i+p[i]>maxr) 
		{
			maxr=i+p[i];pos=i;
		}
	}
	for(long long i=1;i<=len1*2+1;i++)
	{
		if(!p[i]) continue;
		cf[(i-p[i]+1)/2]++;
		if(i%2==0) cf[i/2]--;
		else cf[(i+1)/2]--;
	}
}
int main()
{
	
	jz1[0]=jz2[0]=1;
	for(long long i=1;i<=1000000;i++)
	{
		jz1[i]=(jz1[i-1]*b1)%mod1;
		jz2[i]=(jz2[i-1]*b2)%mod2;
	}
	scanf("%s%s",s+1,t+1);
	len1=strlen(s+1);len2=strlen(t+1);d[len1+1]=-1;
	init();
	manacher();
	for(long long i=1;i<=len1;i++) d[i]=d[i-1]+cf[i];
	reverse(t+1,t+len2+1);
	for(long long i=len1;i>=1;i--)
	{
		h1[i]=(h1[i+1]*b1+s[i])%mod1;
		h2[i]=(h2[i+1]*b2+s[i])%mod2;
	}
	for(long long i=len2;i>=1;i--)
	{
		h3[i]=(h3[i+1]*b1+t[i])%mod1;
		h4[i]=(h4[i+1]*b2+t[i])%mod2;
	}
	for(long long i=len1;i>=1;i--)
	{
		if(s[i]!=t[len2]) continue;
		long long l=0,r=i;
		long long ans=0;
		while(l<=r)
		{
			long long mid=(l+r)/2;
			if(len2-mid<1)
			{
				r=mid-1;
				continue;
			}
			if((((h1[i-mid]-h1[i+1]*jz1[mid+1]%mod1+mod1)%mod1)==h3[len2-mid])&&(((h2[i-mid]-h2[i+1]*jz2[mid+1]%mod2+mod2)%mod2)==h4[len2-mid]))
			{
				l=mid+1;ans=mid;
			}
			else r=mid-1;
		}
		las+=((ans+1)*(d[i+1]+1));
	}
	cout<<las<<endl;
}