●UOJ 131 [NOI2015] 品酒大会

题链:

http://uoj.ac/problem/131

题解:

网上大多数的方法都是用并查集维护。
这里呢,给出另一种自己YY的解法(但实际上本质差不多吧):

后缀数组,RMQ,单调栈

1).预处理
1].首先对字符串后缀排序,得到 sa[i],rank[i],height[i]
   
2].然后维护出 L[i]:表示在后缀数组中,排名最小(记其排名为 L[i])的后缀与排名为 i的后缀的LCP>=hei[i]
    同理 R[i]:表示在后缀数组中,排名最大(记其排名为 R[i])的后缀与排名为 i的后缀的LCP>hei[i](注意这里没有‘=‘,避免重复)
    这两个数组可以用单调栈 O(n)维护出来。
   
3].把输入的权值数组(A数组)按照排好序的后缀的顺序重新排序,得到 B数组:B[i]=A[sa[i]];
    话句话说吧:B[i]表示 排名第i的后缀的首字母对应的权值
    然后对 B数组建两个 ST表,一个用于询问区间最大值,一个用于询问区间最小值。

2).贡献答案
记 num[i]表示相似值为 i的方案数(询问一),val[i]表示相似值为 i的最大美味度(询问二)。
    
不难发现,L[i],R[i]数组可以理解为 在后缀数组的 [L[i],R[i]]这个区间中,hei[i]为最小值。
那么这意味着,区间[L[i],i-1]中的任意一个后缀与区间[i,R[i]]中的任意一个后缀的 LCP==hei[i]
那么便可以贡献答案了:num[hei[i]]+=(i-L[i])*(R[i]-i+1)

同时查询 ST表,得到区间[L[i],i-1]的最大最小权值 max1,min1,以及区间[i,R[i]]最大最小权值 max2,min2。
那么:val[hei[i]]=max(val[hei[i]],max(max1*max2,min1*min2))

(因为有负数,所以要考虑最小值的相乘)
   
注意到以上贡献的值都是建立在 LCP(最长公共前缀)上的,
但还有 LCP-1,LCP-2,LCP-3......这些没有被贡献到,
所以就要最后反向遍历一下两个答案数组,用后面的去更新前面的:
    num[i]+=num[i+1]; val[i]=max(val[i],val[i+1]);

然后就没了。(第一道独立且一遍 AC的 NOI题目,虽然题有点水,但还是偷乐一下咯。)

代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#define MAXN 300500
#define filein(x) freopen(#x".in","r",stdin);
#define fileout(x) freopen(#x".out","w",stdout);
using namespace std;
char S[MAXN];
int A[MAXN],B[MAXN];
int sa[MAXN],rak[MAXN],hei[MAXN],log2[MAXN],L[MAXN],R[MAXN];
struct ST{
	int st[MAXN][20];
	void Make(int N,int *val,bool type){//0:小 1:大
		for(int i=0;i<N;i++) st[i][0]=val[i];
		for(int k=1;k<=log2[N];k++)
			for(int i=(1<<k)-1;i<N;i++)
				st[i][k]=type?max(st[i-(1<<(k-1))][k-1],st[i][k-1]):
							  min(st[i-(1<<(k-1))][k-1],st[i][k-1]);
	}
	int query(int l,int r,bool type){
		static int k;
		k=log2[r-l+1];
		return type?max(st[l+(1<<k)-1][k],st[r][k]):
					min(st[l+(1<<k)-1][k],st[r][k]);
	}
}BTMAX,BTMIN;
void build(int N,int M){
	static int cc[MAXN],ta[MAXN],tb[MAXN],*x,*y,h,p;
	x=ta; y=tb; h=0;
	for(int i=0;i<M;i++) cc[i]=0;
	for(int i=0;i<N;i++) cc[x[i]=S[i]]++;
	for(int i=1;i<M;i++) cc[i]+=cc[i-1];
	for(int i=N-1;i>=0;i--) sa[--cc[x[i]]]=i;
	for(int k=1;p=0,k<N;k<<=1){
		for(int i=N-k;i<N;i++) y[p++]=i;
		for(int i=0;i<N;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
		for(int i=0;i<M;i++) cc[i]=0;
		for(int i=0;i<N;i++) cc[x[y[i]]]++;
		for(int i=1;i<M;i++) cc[i]+=cc[i-1];
		for(int i=N-1;i>=0;i--) sa[--cc[x[y[i]]]]=y[i];
		swap(x,y); y[N]=-1; x[sa[0]]=0; M=1;
		for(int i=1;i<N;i++)
			x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?M-1:M++;
		if(M>=N) break;
	}
	for(int i=0;i<N;i++) rak[sa[i]]=i;//puts(S+sa[i]);
	for(int i=0,j;i<N;i++){
		if(h) h--;
		if(rak[i]){
			j=sa[rak[i]-1];
			while(S[i+h]==S[j+h]) h++;
		}
		hei[rak[i]]=h;
	}
}
void preLR(int N){
	static int stk[MAXN],stp[MAXN],top;
	top=0; stp[top]=0;
	for(int i=0;i<N;i++){
		while(top&&stk[top]>=hei[i]) top--;
		L[i]=stp[top]; top++;
		stk[top]=hei[i]; stp[top]=i;
	}
	top=N; stp[top]=N;
	for(int i=N-1;i>=0;i--){
		while(top<N&&stk[top]>hei[i]) top++;
		R[i]=stp[top]-1; top--;
		stk[top]=hei[i]; stp[top]=i;
	}
	/*for(int i=0;i<N;i++) printf("%d ",L[i]);printf("\n");
	for(int i=0;i<N;i++) printf("%d ",R[i]);printf("\n");*/
}
void getans(int N){
	static long long num[MAXN],val[MAXN],min1,min2,max1,max2;
	memset(val,0xcc,sizeof(val));
	for(int i=1;i<N;i++){
		num[hei[i]]+=1ll*(i-L[i])*(R[i]-i+1);
		min1=BTMIN.query(L[i],i-1,0);
		min2=BTMIN.query(i,R[i],0);
		max1=BTMAX.query(L[i],i-1,1);
		max2=BTMAX.query(i,R[i],1);
		val[hei[i]]=max(val[hei[i]],max(1ll*min1*min2,1ll*max1*max2));
	}
	for(int i=N-2;i>=0;i--){
		num[i]+=num[i+1];
		val[i]=max(val[i],val[i+1]);
	}
	for(int i=0;i<N;i++){
		if(!num[i]) val[i]=0;
		printf("%lld %lld\n",num[i],val[i]);
	}
}
int main()
{
	//filein(ex_savour); fileout(ex_savour);
	int N; log2[1]=0;
	for(int i=2;i<=300000;i++) log2[i]=log2[i>>1]+1;
	scanf("%d",&N); scanf("%s",S);
	for(int i=0;i<N;i++) scanf("%d",&A[i]);
	build(N,300);
	for(int i=0;i<N;i++) B[i]=A[sa[i]];
	BTMIN.Make(N,B,0);
	BTMAX.Make(N,B,1);
	preLR(N); getans(N);
	return 0;
}
时间: 2024-10-11 16:56:14

●UOJ 131 [NOI2015] 品酒大会的相关文章

[UOJ#131][BZOJ4199][NOI2015]品酒大会 后缀数组 + 并查集

[UOJ#131][BZOJ4199][NOI2015]品酒大会 试题描述 一年一度的“幻影阁夏日品酒大会”隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发“首席品酒家”和“首席猎手”两个奖项,吸引了众多品酒师参加. 在大会的晚餐上,调酒师 Rainbow 调制了 n杯鸡尾酒.这 n杯鸡尾酒排成一行,其中第 i杯酒 (1≤i≤n ) 被贴上了一个标签 si ,每个标签都是 26 个小写英文字母之一.设 Str(l,r)表示第 l杯酒到第 r 杯酒的 r−l+1 个标签顺次连接构成的字

【BZOJ4199】[Noi2015]品酒大会 后缀数组+并查集

[BZOJ4199][Noi2015]品酒大会 题面:http://www.lydsy.com/JudgeOnline/wttl/thread.php?tid=2144 题解:听说能用SAM?SA默默水过~ 本题的实现还是非常简单的,先求出height数组,然后两杯酒'r'相似就等价于二者中间的height都>=r,于是我们将height排序,从大到小扔进去,那么所有连续的区间中的点对就都是相似的了.维护连续区间可以用并查集.统计点对个数需要维护size,统计最大乘积需要维护最(次)大(小)值,

【uoj131】 NOI2015—品酒大会

http://uoj.ac/problem/131 (题目链接) 题意 给出一个字符串,每个后缀有一个权值${a_i}$,这些后缀两两之间存在公共前缀.问能够组成长度从0~n-1的公共前缀的后缀的方案数以及他们权值的最大乘积. Solution 听LCF说这是水题,就来做了.. lyp学长说SAM构出来后就两个东西:在自动机上dp,在后缀树上搞事情. 于是我们很容易想到在后缀树上dp(感觉在自动机上dp没什么卵用),于是我们把串反过来见后缀自动机,然后建出后缀树,在上面树形dp就可以了. 假设当

BZOJ 4199 [Noi2015]品酒大会:后缀数组 + 并查集

题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=4199 题意: 给你一个长度为n的字符串s,和一个长为n的数组v. 对于每个整数r∈[0,n-1]: (1)问你有多少对后缀(suffix(i), suffix(j)),满足LCP(suffix(i), suffix(j)) >= r (2)输出mul[r] = max(v[i]*v[j]),其中i,j满足(1)的条件 题解: 先考虑第(1)问. 由于LCP只受连续的一段height中最小

BZOJ 4199: [Noi2015]品酒大会( 后缀数组 + 并查集 )

求出后缀数组后, 对height排序, 从大到小来处理(r相似必定是0~r-1相似), 并查集维护. 复杂度O(NlogN + Nalpha(N)) ----------------------------------------------------------------------------------- #include<cstdio> #include<cstring> #include<algorithm> using namespace std; ty

字符串(后缀数组||SAM):NOI2015 品酒大会

这道题可以转化为计数类问题. 若使用后缀数组,那么答案就是所有位置二元组(i,j)的lcp对0~lcp答案段的贡献.然后发现若一个二元组有x的贡献,那么对x-1有同样的贡献,考虑先求出lcp(max)的答案,再传给lcp(max-1)等等,复杂度是O(N)的. 若用SAM,那么需要求的答案在x与fa[x]的转移之间,因为后缀自动机中每个点所代表的出现位置和出现次数是相等的,可以合并,事实上还是用到了后缀数组的叠加的原理,从叶子节点一路传上去…………………………………………… 还有要注意INF要足

BZOJ4199: [Noi2015]品酒大会

后缀数组height排序后并查集合并 也就是height较大的合并不影响较小的 num[i]=num[i+1]  ans[i]=ans[i+1] 合并时,num+=sz[x]*sz[y],ans=max(mn[x]*mn[y],mx[x]*mx[y],ans) 这种思路适应于求点对,还可以考虑启发式合并 后缀数组还有的常见思路就是二分,height分组 1 #include<bits/stdc++.h> 2 using namespace std; 3 #define N 300005 4 #

Bzoj4199:[NOI2015]品酒大会

题面 Bzoj4199 Sol 后缀数组 显然的暴力就是求\(LCP\)+差分 \(40\)分 # include <bits/stdc++.h> # define RG register # define IL inline # define Fill(a, b) memset(a, b, sizeof(a)) using namespace std; typedef long long ll; const int _(3e5 + 5); IL int Input(){ RG int x =

bzoj4199: [Noi2015]品酒大会 (并查集 &amp;&amp; 后缀数组)

据说用后缀自动机 + dp也能做 然而并不会 后缀数组的做法呢 就是先建个后缀数组,求出height值,此时如果直接找,复杂度是n ^ 2的,肯定会超时. 但是height大的值是不会对小的产生影响的,所以可以按height大小,从大到小合并两个区间,用并查集维护就可以了 代码如下 1 #include <cstdio> 2 #include <algorithm> 3 #include <cstring> 4 using namespace std; 5 const