并不对劲的bzoj4560:p3269:[JLOI2016]字符串覆盖

题目大意

(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
*/

一些感想

片手玩得好不是大佬,笛玩得好才是。——吹笛子时被车的片手流

原文地址:https://www.cnblogs.com/xzyf/p/10410845.html