加油鸭!
比赛自闭,补题吃力,新算法不深,基础算法不透。
写个计划,博客打卡,督促自己。
一周一专题,一至五按计划写题,六七补题。
题目来自 kuangbin带我飞系列。
欢迎大家找我讨论。
当刷题日记写,水题直接贴,经典题难题单独写博客。
每日更新。
第一周:并查集(鸽)
这周我没了,后期补上。
第二周:KMP & 扩展KMP & Manacher
周一:[ KMP ]
A - Number Sequence
#include<cstdio>
const int Mn=1e6+17;
void preKMP(int x[],int m,int kmpNext[]){
int i=0,j=-1;kmpNext[0]=-1;
while(i<m){
while(-1!=j&&x[i]!=x[j]) j=kmpNext[j];
if(x[++i]==x[++j])kmpNext[j]=kmpNext[j];
else kmpNext[i]=j;
}
}
int next[Mn];
int KMP(int x[],int m,int y[],int n){ //x模式串,y主串
int i=0,j=0,ans=-1;
preKMP(x,m,next);
while(i<n){
while(-1!=j&&y[i]!=x[j]) j=next[j];
++i,++j;
if(j>=m){ans=i-m+1;break;}
}
return ans;
}
int a[1000017],b[10017];
int main(){
int T;scanf("%d",&T);
while(T--){
int n,m;scanf("%d%d",&n,&m);
for(int i=0;i<n;++i)scanf("%d",&a[i]);
for(int j=0;j<m;++j)scanf("%d",&b[j]);
printf("%d\n",KMP(b,m,a,n));
}
}
B - Oulipo
#include<cstdio>
#include<cstring>
const int Mn=1e5+17;
void preKMP(char x[],int m,int kmpNext[]){
int i=0,j=-1;kmpNext[0]=-1;
while(i<m){
while(-1!=j&&x[i]!=x[j]) j=kmpNext[j];
if(x[++i]==x[++j])kmpNext[j]=kmpNext[j];
else kmpNext[i]=j;
}
}
int next[Mn];
int KMP(char x[],int m,char y[],int n){ //x模式串,y主串
int i=0,j=0,ans=0;
preKMP(x,m,next);
while(i<n){
while(-1!=j&&y[i]!=x[j]) j=next[j];
++i,++j;
if(j>=m){++ans;j=next[j];}
}
return ans;
}
char a[10017],b[1000017];
int main(){
int T;scanf("%d",&T);
while(T--){
scanf("%s%s",a,b);
printf("%d\n",KMP(a,strlen(a),b,strlen(b)));
}
return 0;
}
C - 剪花布条
#include<cstdio>
#include<cstring>
const int Mn=1e5+17;
void preKMP(char x[],int m,int kmpNext[]){
int i=0,j=-1;kmpNext[0]=-1;
while(i<m){
while(-1!=j&&x[i]!=x[j]) j=kmpNext[j];
if(x[++i]==x[++j])kmpNext[j]=kmpNext[j];
else kmpNext[i]=j;
}
}
int next[Mn];
int KMP(char x[],int m,char y[],int n){ //x模式串,y主串
int i=0,j=0,ans=0;
preKMP(x,m,next);
while(i<n){
while(-1!=j&&y[i]!=x[j]) j=next[j];
++i,++j;
if(j>=m){++ans;j=0;}
}
return ans;
}
char a[10017],b[1000017];
int main(){
while(~scanf("%s",a)){
if(a[0]=='#') break;
scanf("%s",b);
printf("%d\n",KMP(b,strlen(b),a,strlen(a)));
}
return 0;
}
D - Cyclic Nacklace
next[ ]应用。
藏在这里。
E - Period (×)
F - Power Strings (×)
G - Seek the Name, Seek the Fame (×)
周三:[ KMP & Manacher ]
A - How many
最小表示法。
藏在这里。
B - Period II
这题我看了一个中午都没明白case4,后来9f巨巨教会我了(%%%)。
也是
n
e
x
t
[
]
next[]
next[] 数组的应用。
之前有篇博客写过 len−next[i]
为此字符串的最小循环节,这道题求(太难文字叙述了,略),同样用到这个知识点,把
n
e
x
t
[
]
回
退
next[]回退
next[]回退 next[next[next[len-next[len]]]]
,大概长这样,就可以求出所有循环节了。
其实我没弄太懂是为什么。
我弄懂了再写博客吧。
目前思路是这样
n
e
s
t
[
i
]
nest[i]
nest[i]存的不是失配后,跳转匹配的位置嘛,说明
i
i
i 前面部分子串是和串串的某部分匹配成功的。
index初始化为最小循环节index=len−next[i]
……
这题要注意不要忽略除了最小循环节以外的循环节。
还有就是,用容器存答案会超时。emmm,不都是nlogn嘛,也不懂为什么。
注意:next数组不要命名为next[],因为next是标识符,会报错。
#include<cstdio>
#include<set>
#include<cstring>
#include<algorithm>
using namespace std;
const int Mn=1e6+17;
void kmp_pre(char x[],int m,int nextt[]){
int i=0,j=-1;nextt[0]=-1;
while(i<m){
while(-1!=j&&x[i]!=x[j]) j=nextt[j];
nextt[++i]=++j;
}
// for(int i=0;i<=m;++i)printf("%2d ",i);printf("\n");
// for(int i=0;i<m;++i)printf(" %c ",x[i]);printf("\n");
// for(int i=0;i<=m;++i)printf("%2d ",nextt[i]);printf("\n");
}
int nextt[Mn],res[Mn];
// set<int> ans;
char a[Mn];
int main(){
int T;scanf("%d",&T);
for(int tt=1;tt<=T;++tt){
// ans.clear();
scanf("%s",a);
int len=strlen(a);
// preKMP(a,len,nextt); //这个是假的next[]数组
kmp_pre(a,len,nextt);
int index=len,tail=1;res[0]=len;
// ans.insert(len);
while(nextt[index]>0){
// ans.insert(len-nextt[index]);
res[tail++]=len-nextt[index];
index=nextt[index];
}
sort(res,res+tail);
printf("Case #%d: %d\n",tt,tail);
// printf("Case #%d: %d\n",tt,(int)ans.size());
// set<int>::iterator it;
// for(auto i: ans) printf("%d ",i);printf("\n");
for(int i=0;i<tail-1;++i)printf("%d ",res[i]);printf("%d\n",res[tail-1]);
}
return 0;
}
C - Teacher YYF
恶心打表题,略过。
D - Best Reward
马拉车。
藏在这里。
E - Finding Palindromes(需补)
没什么思路。需要学字典树。
周四:[ Manacher ]
A - Palindrome
manacher板子题。
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int Mn=1000017;
const int INF=0x3f3f3f3f;
char s[Mn],Ma[Mn*2];
int Mp[2*Mn];
int Manacher(char s[],int len){
int maxx=-INF;
int l=2; Ma[0]='$';Ma[1]='#';
for(int i=0;i<len;++i) Ma[l++]=s[i],Ma[l++]='#';
Ma[l]=0;
int mx=0,id=0;
for(int i=0;i<l;++i){
Mp[i]=mx>i?min(Mp[2*id-i],mx-i):1;
while(Ma[i+Mp[i]]==Ma[i-Mp[i]]) Mp[i]++;
if(i+Mp[i]>mx)mx=i+Mp[i],id=i;
}
for(int i=1;i<l;++i) maxx=max(maxx,Mp[i]);
return maxx-1;
}
int main(){
int Case=1;
while(~scanf("%s",s)){
if(s[0]=='E') break;
printf("Case %d: %d\n",Case++,Manacher(s,strlen(s)));
}
return 0;
}
B - 吉哥系列故事——完美队形II
manacher板子题,求单峰回文(中点最大值)。Manacher里多一个判断。
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int Mn=1000017;
const int INF=0x3f3f3f3f;
int s[Mn],Ma[Mn*2];
int Mp[2*Mn];
int Manacher(int s[],int len){
int maxx=-INF;
int l=2; Ma[0]=-2;Ma[1]=-1;
for(int i=0;i<len;++i) Ma[l++]=s[i],Ma[l++]=-1;
Ma[l]=-3;
int mx=0,id=0;
for(int i=0;i<l;++i){
Mp[i]=mx>i?min(Mp[2*id-i],mx-i):1;
while(Ma[i+Mp[i]]==Ma[i-Mp[i]]&&Ma[i-Mp[i]]<=Ma[i-Mp[i]+2]) Mp[i]++; //多加一个判断
if(i+Mp[i]>mx)mx=i+Mp[i],id=i;
}
for(int i=1;i<l;++i) maxx=max(maxx,Mp[i]);
return maxx-1;
}
int main(){
int T,n;scanf("%d",&T);
while(T--){
scanf("%d",&n);
for(int i=0;i<n;++i)scanf("%d",&s[i]);
printf("%d\n",Manacher(s,n));
}
return 0;
}
C - Girls’ research
manacher板子题。多了个字母转化。输入很不好搞,我wa了3发。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int Mn=1000017;
const int INF=0x3f3f3f3f;
char s[Mn],Ma[Mn*2],ans[Mn];
int Mp[2*Mn];
int Manacher(char s[],int len){
int index=0,maxx=-INF;
int l=2; Ma[0]='$';Ma[1]='#';
for(int i=0;i<len;++i) Ma[l++]=s[i],Ma[l++]='#';
Ma[l]=0;
int mx=0,id=0;
for(int i=0;i<l;++i){
Mp[i]=mx>i?min(Mp[2*id-i],mx-i):1;
while(Ma[i+Mp[i]]==Ma[i-Mp[i]]) Mp[i]++;
if(i+Mp[i]>mx)mx=i+Mp[i],id=i;
}
for(int i=1;i<l;++i){
if(maxx<Mp[i]){
index=i;
maxx=Mp[i];
}
}
return index;
}
int main(){
char w;
while(cin>>w>>s){
// while(~scanf("%c",&w)){ //注释部分为错误代码(因为getchar)
// getchar();
// scanf("%s",s);
memset(Mp,0,sizeof(Mp));
int index=Manacher(s,strlen(s)),tail=0,p1=0,p2=0;
bool flag=true;
if(Mp[index]<=2) {puts("No solution!");continue;}
for(int i=index-Mp[index]+1;i<index+Mp[index];++i){
if(Ma[i]!='#'){
ans[tail++]=Ma[i]-(w-'a');
p2=i;
if(flag) {p1=i;flag=false;}
}
}
printf("%d %d\n",p1/2-1,p2/2-1);
for(int i=0;i<tail;++i)printf("%c",(ans[i]<'a')?(ans[i]+26):ans[i]);printf("\n");
}
return 0;
}
周五
A - 最长回文
manacher板中板题。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int Mn=110017;
const int INF=0x3f3f3f3f;
char s[Mn],Ma[Mn*2];
int Mp[2*Mn];
int Manacher(char s[],int len){
int maxx=-INF;
int l=2; Ma[0]='$';Ma[1]='#';
for(int i=0;i<len;++i) Ma[l++]=s[i],Ma[l++]='#';
Ma[l]=0;
int mx=0,id=0;
for(int i=0;i<l;++i){
Mp[i]=mx>i?min(Mp[2*id-i],mx-i):1;
while(Ma[i+Mp[i]]==Ma[i-Mp[i]]) Mp[i]++;
if(i+Mp[i]>mx)mx=i+Mp[i],id=i;
}
for(int i=1;i<l;++i)maxx=max(maxx,Mp[i]);
return maxx-1;
}
int main(){
while(~scanf("%s",s)){
memset(Mp,0,sizeof(Mp));
printf("%d\n",Manacher(s,strlen(s)));
}
return 0;
}
B - Wow! Such Doge!
挺傻的一题。(要不是因为它是西安邀请赛的题,我就跳过了)
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int Mn=110017;
char s[Mn];
int main(){
int cnt=0;
while(~scanf("%s",s)){
int len=strlen(s);
for(int i=0;i<len-3;++i){
if((s[i]=='d'||s[i]=='D')&&(s[i+1]=='o'||s[i+1]=='O')&&(s[i+2]=='g'||s[i+2]=='G')&&(s[i+3]=='e'||s[i+3]=='E'))
++cnt;
}
}
printf("%d\n",cnt);
return 0;
}
C - Theme Section
在一个大串串,中间找一段子串,满足和串串的前缀与后缀都相同,求这个子串最长长度(不能与前缀后缀重合)。
嘤嘤嘤,是个水题,但我对next数组理解不深,debug了好久。
这个题加深了我对next的理解,感觉前面那题可以写博客了。
#include<cstdio>
#include<cstring>
const int Mn=1e6+17;
char s[Mn];
int nxt[Mn];
void kmp_pre(char x[],int m,int nxt[]){
int i=0,j=-1;nxt[0]=-1;
while(i<m){
while(-1!=j&&x[i]!=x[j]) j=nxt[j];
nxt[++i]=++j;
}
// for(int i=0;i<m;++i)printf(" %c ",s[i]);printf("\n");
// for(int i=0;i<=m;++i)printf("%2d ",nxt[i]);printf("\n");
}
int KMP(char x[],int m){
for(int i=nxt[m];i!=-1;i=nxt[i]){
// printf("i=%d,m=%d,next[%d]=%d\n",i,m,m,nxt[m]);
for(int j=i+1;j<=m-i;++j){
// printf("i=%d,j=%d,Next[%d]=%d\n",i,j,j,nxt[j]);
if(nxt[j]==i) return i;
}
}
return 0;
}
int main(){
int T;scanf("%d",&T);
while(T--){
scanf("%s",s);
kmp_pre(s,strlen(s),nxt);
printf("%d\n",KMP(s,strlen(s)));
}
return 0;
}
第三周:AC自动机
周一
A - Keywords Search
板子。
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
struct Trie{
int nxt[500010][26],fail[500010],end[500010];
int root,L;
int newnode(){
for(int i=0;i<26;++i) nxt[L][i]=-1;
end[L++]=0;
return L-1;
}
void init(){
L=0;
root=newnode();
}
void insert(char buf[]){
int len=strlen(buf);
int now=root;
for(int i=0;i<len;++i){
if(nxt[now][buf[i]-'a']==-1)
nxt[now][buf[i]-'a']=newnode();
now=nxt[now][buf[i]-'a'];
}
end[now]++;
}
void build(){
queue<int> Q;
fail[root]=root;
for(int i=0;i<26;++i)
if(nxt[root][i]==-1) nxt[root][i]=root;
else {
fail[nxt[root][i]]=root;
Q.push(nxt[root][i]);
}
while(!Q.empty()){
int now=Q.front();
Q.pop();
for(int i=0;i<26;++i)
if(nxt[now][i]==-1) nxt[now][i]=nxt[fail[now]][i];
else {
fail[nxt[now][i]]=nxt[fail[now]][i];
Q.push(nxt[now][i]);
}
}
}
int query(char buf[]){
int len=strlen(buf);
int now=root;
int res=0;
for(int i=0;i<len;++i){
now=nxt[now][buf[i]-'a'];
int temp=now;
while(temp!=root){
res+=end[temp];
end[temp]=0;
temp=fail[temp];
}
}
return res;
}
void debug(){
printf("id fail end chi ");
for(int i=0;i<26;++i)printf(" %c",'a'+i);printf("\n");
for(int i=0;i<L;++i){
printf("id=%3d,fail=%3d,end=%3d,chi=[",i,fail[i],end[i]);
for(int j=0;j<26;++j) if(nxt[i][j]!=0)printf("%2d",nxt[i][j]); else printf(" _");printf("]\n");
}
}
};
char buf[1000010];
Trie ac;
int main(){
int T;scanf("%d",&T);
while(T--){
int n;scanf("%d",&n);
ac.init();
for(int i=0;i<n;++i){
scanf("%s",buf);
ac.insert(buf);
}
ac.build();
scanf("%s",buf);
printf("%d\n",ac.query(buf));
}
return 0;
}
B - 病毒侵袭
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
const int Mn=200*500+17;
const int ALTN=128;
struct Trie{
int nxt[Mn][ALTN+7];
int fail[Mn];
int end[Mn];
int root,L;
int newnode(){
for(int i=0;i<ALTN;++i) nxt[L][i]=-1;
end[L++]=-1;
return L-1;
}
void init(){
L=0;
root=newnode();
}
void insert(char buf[],int id){
int len=strlen(buf);
int now=root;
for(int i=0;i<len;++i){
if(nxt[now][(int)buf[i]]==-1)
nxt[now][(int)buf[i]]=newnode();
now=nxt[now][(int)buf[i]];
}
end[now]=id;
}
void build(){
queue<int> Q;
fail[root]=root;
for(int i=0;i<ALTN;++i)
if(nxt[root][i]==-1) nxt[root][i]=root;
else {
fail[nxt[root][i]]=root;
Q.push(nxt[root][i]);
}
while(!Q.empty()){
int now=Q.front();Q.pop();
for(int i=0;i<ALTN;++i)
if(nxt[now][i]==-1)
nxt[now][i]=nxt[fail[now]][i];
else {
fail[nxt[now][i]]=nxt[fail[now]][i];
Q.push(nxt[now][i]);
}
}
}
int query(char buf[],int id,int tn){
int len=strlen(buf);
int now=root;
bool book[517],ans=false; for(int i=1;i<=tn;++i)book[i]=false;
for(int i=0;i<len;++i){
now=nxt[now][(int)buf[i]];
int temp=now;
while(temp!=root){
if(end[temp]!=-1) {book[end[temp]]=true;ans=true;}
temp=fail[temp];
}
}
if(!ans) return 0;
printf("web %d:",id);for(int i=1;i<=tn;++i)if(book[i]) printf(" %d",i);printf("\n");
return 1;
}
// void debug(){
// printf("id fail end chi ");
// for(int i=97;i<123;++i)printf(" %c",i);printf("\n");
// for(int i=0;i<L;++i){
// printf("id=%3d,fail=%3d,end=%3d,chi=[",i,fail[i],end[i]);
// for(int j=97;j<123;++j) if(nxt[i][j]!=0)printf("%2d",nxt[i][j]); else printf(" _");printf("]\n");
// }
// }
};
char buf[10017];
Trie ac;
int main(){
int n,m;scanf("%d",&n);
ac.init();
for(int i=1;i<=n;++i){
scanf("%s",buf);
ac.insert(buf,i);
}
//ac.debug();
ac.build();
scanf("%d",&m);
int cnt=0;
for(int i=1;i<=m;++i){
scanf("%s",buf);
cnt+=ac.query(buf,i,n);
}
printf("total: %d\n",cnt);
return 0;
}
C - 病毒侵袭持续中
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
const int Mn=1000*50+100;
const int ALTN=128;
char buf[2000017],s[1017][57];
struct Trie{
int nxt[Mn][ALTN+7];
int fail[Mn];
int end[Mn];
int root,L;
int newnode(){
for(int i=0;i<ALTN;++i) nxt[L][i]=-1;
end[L++]=-1;
return L-1;
}
void init(){
L=0;
root=newnode();
}
void insert(char buf[],int id){
int len=strlen(buf);
int now=root;
for(int i=0;i<len;++i){
if(nxt[now][(int)buf[i]]==-1)
nxt[now][(int)buf[i]]=newnode();
now=nxt[now][(int)buf[i]];
}
end[now]=id;
}
void build(){
queue<int> Q;
fail[root]=root;
for(int i=0;i<ALTN;++i)
if(nxt[root][i]==-1) nxt[root][i]=root;
else {
fail[nxt[root][i]]=root;
Q.push(nxt[root][i]);
}
while(!Q.empty()){
int now=Q.front();Q.pop();
for(int i=0;i<ALTN;++i)
if(nxt[now][i]==-1)
nxt[now][i]=nxt[fail[now]][i];
else {
fail[nxt[now][i]]=nxt[fail[now]][i];
Q.push(nxt[now][i]);
}
}
}
void query(char buf[],int tn){
int len=strlen(buf);
int now=root;
int book[1017]; for(int i=1;i<=tn;++i)book[i]=0;
for(int i=0;i<len;++i){
now=nxt[now][(int)buf[i]];
int temp=now;
while(temp!=root){
if(end[temp]!=-1) book[end[temp]]++;
temp=fail[temp];
}
}
for(int i=1;i<=tn;++i){if(book[i]) printf("%s: %d\n",s[i],book[i]);}
}
// void debug(){
// printf("id fail end chi ");
// for(int i=97;i<123;++i)printf(" %c",i);printf("\n");
// for(int i=0;i<L;++i){
// printf("id=%3d,fail=%3d,end=%3d,chi=[",i,fail[i],end[i]);
// for(int j=97;j<123;++j) if(nxt[i][j]!=0)printf("%2d",nxt[i][j]); else printf(" _");printf("]\n");
// }
// }
};
Trie ac;
int main(){
int n;
while(~scanf("%d",&n)){
ac.init();
for(int i=1;i<=n;++i){
scanf("%s",s[i]);
ac.insert(s[i],i);
}
//ac.debug();
ac.build();
scanf("%s",buf);
ac.query(buf,n);
}
return 0;
}
D - DNA Sequence
AC自动机+矩阵快速幂。
藏在这里。