待学知识点:
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;
}