题目大意
(T)((Tleq10))组询问
每组询问给出一个字符串(A)((|A|leq10^4)),(n)((nleq4))个(A)的子串(B_1,B_2,B_3,...,B_n)((forall i in[1,n],|B_i|leq10^3))
如果(|B_i|=r-l+1)且(B_i)的每一个字符依次与(A_l,A_{l+1},...,A_{r})相等,那么区间([l,r])每个位置都被覆盖了一次
求如果必须把(B_1,...,B_n)全部用上,那么字符串(A)至少有几个位置被覆盖不少于一次?至多有几个位置被覆盖不少于一次?
题解
至少:
如果一个串(B_i)是(B_j)的子串,那么让(B_i)覆盖被(B_j)覆盖的部分会更优,所以先去掉“是别的串的子串的串”
设(f(S,j))表示已经用过的子串集合为(S),目前覆盖的最靠右的位置为(j)
因为已经去掉了“是别的串的子串的串”,所以以(j)为结尾的(B)串只有一个,假设它是(B_x)
有两种情况:1.在选取(B_x)前(此时已用过的子串集合为(Sigoplus{ x})),被覆盖的区域与(B_x)有交;2.在选取(B_x)前,被覆盖的区域与(B_x)没有交
1.(f(S,j)=min{ f(Sigoplus{ x},k)+j-k mid j-|B_x|+1leq kleq j})因为已经去重,所以集合(Sigoplus{ x})里最靠后的子串的起点一定不会比(j-|B_x|+1)更靠右,加入(B_x)后新覆盖的区域为(j-k)
因为没有“是别的串的子串的串”,所以当(j)右移时(j-|B_x|+1)也右移,用单调队列维护
2.(f(S,j)=min{ f(Sigoplus{ x},k) mid 0leq kleq j-|B_x|+1}+|B_x|),维护前缀最小值就行
时间复杂度(Theta(2^n*|A|*n))或(Theta(2^n*|A|))
至多:
发现并不能去重
暴力枚举(n)个串的顺序,暴力枚举每个串与之前的部分是否有交
有交就让交的部分尽可能小,没有交就尽可能把这个串往前放
时间复杂度大概(Theta(n!*2^n*n))
代码
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<iomanip>
#include<iostream>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<vector>
#define rep(i,x,y) for(register int i=(x);i<=(y);++i)
#define dwn(i,x,y) for(register int i=(x);i>=(y);--i)
#define maxn 10010
#define maxm 1010
#define pii pair<int,int>
#define fi first
#define se second
#define mp(x,y) make_pair(x,y)
#define maxs (1<<5)
#define chmn(x,y) x=min(x,y)
using namespace std;
int read()
{
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)&&ch!='-')ch=getchar();
if(ch=='-')f=-1,ch=getchar();
while(isdigit(ch))x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
return x*f;
}
void write(int x)
{
if(x==0){putchar('0');return;}
int f=0;char ch[20];
if(x<0)putchar('-'),x=-x;
while(x)ch[++f]=x%10+'0',x/=10;
while(f)putchar(ch[f--]);
return;
}
int T,pos[5][maxn],n,lens,lent[5],fa[5][maxm],ord[20],ord2[20],vis[5],f[maxs][maxn],maxans,minans,cntnd,hsh[5],ext[5][maxn],pref[maxs][maxn];
char t[5][maxm],s[maxn];
pii a[20];
struct deque
{
int w[maxn],tim[maxn],hd,tl;
void reset(){hd=1,tl=0;}
void push(int x,int y){while(hd<=tl&&w[tl]>=x)tl--;w[++tl]=x,tim[tl]=y;}
void del(int y){while(hd<=tl&&tim[hd]<y)hd++;}
int empty(){return hd>tl?1:0;}
}q[maxs];
bool cmp1(int x,int y){return lent[x]<lent[y];}
bool cmp2(int x,int y){return a[x].fi==a[y].fi?a[x].se<a[y].se:a[x].fi<a[y].fi;}
int calmx()
{
rep(i,1,cntnd)ord2[i]=i;
sort(ord2+1,ord2+cntnd+1,cmp2);
int now=0,lenth=0;
rep(i,1,cntnd)
{
now+=a[ord2[i]].se;
if(i==cntnd||a[ord2[i+1]].fi!=a[ord2[i]].fi){if(now)lenth+=a[ord2[i+1]].fi-a[ord2[i]].fi;}
}
return lenth;
}
int getpla(int step,int maxp)
{
if(step==n+1){return calmx();}
int now=ord[step],res=0;
int l=1,r=pos[now][0],ans=0;
while(l<=r)
{
int mid=(l+r)>>1;
if(pos[now][mid]-lent[now]+1<=maxp){l=mid+1;ans=max(ans,mid);}
else r=mid-1;
}
if(ans!=0&&pos[now][ans]>=maxp)a[++cntnd]=mp(pos[now][ans]+1,-1),a[++cntnd]=mp(pos[now][ans]-lent[now]+1,1),res=getpla(step+1,pos[now][ans]),cntnd-=2;
if(ans+1!=pos[now][0]+1&&pos[now][ans+1]>=maxp)
a[++cntnd]=mp(pos[now][ans+1]+1,-1),a[++cntnd]=mp(pos[now][ans+1]-lent[now]+1,1),
res=max(res,getpla(step+1,pos[now][ans+1])),cntnd-=2;
return res;
}
void getord(int step)
{
if(step==n+1){cntnd=0;maxans=max(maxans,getpla(1,0));return;}
rep(i,1,n)if(!vis[i])vis[i]=1,ord[step]=i,getord(step+1),vis[i]=0;
return;
}
int main()
{
T=read();
while(T--)
{
scanf("%s",s+1);lens=strlen(s+1);maxans=0;
n=read();memset(ext,0,sizeof(ext));
rep(i,1,n)
{
pos[i][0]=vis[i]=0,scanf("%s",t[i]+1),lent[i]=strlen(t[i]+1);
fa[i][0]=-1,fa[i][1]=0;
rep(j,2,lent[i])
{
int k=fa[i][j-1];
while(k&&t[i][k+1]!=t[i][j])k=fa[i][k];
if(t[i][k+1]==t[i][j])fa[i][j]=k+1;
else fa[i][j]=0;
}
int k=0;
rep(j,1,lens)
{
while(k&&(k==lent[i]||t[i][k+1]!=s[j]))k=fa[i][k];
if(t[i][k+1]==s[j]){k++;if(k==lent[i])pos[i][++pos[i][0]]=j,ext[i][j]=1;}
else k=0;
}
}
getord(1);
rep(i,1,n)ord[i]=i,vis[i]=0;
sort(ord+1,ord+n+1,cmp1);
rep(i,1,n)
rep(j,i+1,n)
{
int fl=0,k=0;
rep(l,1,lent[ord[j]])
{
while(k&&(t[ord[i]][k+1]!=t[ord[j]][l]))k=fa[ord[i]][k];
if(t[ord[i]][k+1]==t[ord[j]][l]){k++;if(k==lent[ord[i]]){fl=1;break;}}
else k=0;
}
vis[ord[i]]|=fl;
}
cntnd=0;
rep(i,1,n)if(!vis[i])hsh[++cntnd]=i;
int fulls=(1<<cntnd)-1;
rep(i,0,fulls)q[i].reset();
memset(f,0x7f,sizeof(f)),memset(pref,0x7f,sizeof(pref));
int inf=minans=f[0][0];f[0][0]=pref[0][0]=0;q[0].push(0,0);
rep(j,1,lens)rep(s,0,fulls)
{
int num=0;pref[s][j]=pref[s][j-1];
if(!s)continue;
rep(i,1,cntnd)
{
if(ext[hsh[i]][j]&&((1<<(i-1))&s))
{
num++;
int pres=s^(1<<(i-1));
q[pres].del(j-lent[hsh[i]]);
if(!q[pres].empty())f[s][j]=min(f[s][j],q[pres].w[q[pres].hd]+j);
if(j-lent[hsh[i]]>=0)f[s][j]=min(f[s][j],pref[pres][j-lent[hsh[i]]]+lent[hsh[i]]);
pref[s][j]=min(pref[s][j],f[s][j]);if(f[s][j]!=inf)q[s].push(f[s][j]-j,j);
}
}
}
rep(j,1,lens)chmn(minans,f[fulls][j]);
write(minans),putchar(' '),write(maxans),putchar('
');
}
return 0;
}
/*
2
hello
4
he
l
l
o
abacaba
4
ab
ba
a
c
*/
一些感想
片手玩得好不是大佬,笛玩得好才是。——吹笛子时被车的片手流