【BZOJ 4199】 [Noi2015]品酒大会
【链接】h在这里写链接
【题意】
标志了每一杯酒;
以及n个数字,表示每一杯酒的美味度ai。
两杯酒(i,j)称为r相似
当且仅当s的后缀i和后缀j的lcp>=r
然后问你r相似的酒里面任意两个酒的ai*aj的最大值.
以及r相似的酒的对数 (i<j)
【题解】
如果连续出现了一段height值大于k的。
就说明这一段后缀的lcp至少为k
->肯定是k相似的酒。
那么我们就从height值最大的考虑。
单个height单个height地考虑。
这个height如果为k;
那么因为是从大到小的。
所以之前的height肯定>=k;
那么如果我们此时把i和i-1并在一起。
我们并起来的这两个集合里的任意两个后缀肯定都是k相似的。
->因为它中间会有一个height[i]=k,那么根据lcp的计算规则。
两个集合里的后缀任意两两组合都是lcp==k的。
然后我们再把k+1相似的答案累加一下。就能求出所有的答案了。
并查集维护一下最大值、最小值就好.
这种单个ai*aj的问题
感觉都能用(a1+a2..) * (a3 + a4...)的形式来搞
求出height数组;
按照height数组从大到小排;
->记录每个height是排名第几的后缀的height;
然后从最大的height开始逆序枚举;
答案累加上相似值为i+1的对数,最大值和之前的一样。
把这个height[i]的i和i-1合并在一起。
对数递增两个集合的大小,最大值尝试更改一下。
这个集合是排名的集合。
->排名第i一开始单独一个集合,大小为1,最大值、最小值都为val[Sa[i]];
【错的次数】
【反思】
【代码】
#include<bits/stdc++.h>
using namespace std;
#define ll long long const int N = 3e5;
const int MAX_CHAR = 255;//每个数字的最大值。
char s[N + 10];//如果是数字,就写成int s[N+10]就好,从0开始存
int Sa[N + 10], T1[N + 10], T2[N + 10], C[N+10];
int Height[N + 10], Rank[N + 10];
ll a[N+10],sum[N+10],ans[N+10];
ll cnt[N+10],zd[N+10],zx[N+10];
int f[N+10];
pair <int,int> temp[N+10]; void build_Sa(int n, int m) {
int i, *x = T1, *y = T2;
for (i = 0; i<m; i++) C[i] = 0;
for (i = 0; i<n; i++) C[x[i] = s[i]]++;
for (i = 1; i<m; i++) C[i] += C[i - 1];
for (i = n - 1; i >= 0; i--) Sa[--C[x[i]]] = i;
for (int k = 1; k <= n; k <<= 1)
{
int p = 0;
for (i = n - k; i<n; i++) y[p++] = i;
for (i = 0; i<n; i++) if (Sa[i] >= k) y[p++] = Sa[i] - k;
for (i = 0; i<m; i++) C[i] = 0;
for (i = 0; i<n; i++) C[x[y[i]]]++;
for (i = 1; i<m; i++) C[i] += C[i - 1];
for (i = n - 1; i >= 0; i--) Sa[--C[x[y[i]]]] = y[i];
swap(x, y);
p = 1; x[Sa[0]] = 0;
for (i = 1; i<n; i++)
x[Sa[i]] = y[Sa[i - 1]] == y[Sa[i]] && y[Sa[i - 1] + k] == y[Sa[i] + k] ? p - 1 : p++;
if (p >= n) break;
m = p;
}
} void getHeight(int n)
{
int i, j, k = 0;
for (i = 1; i <= n; i++) Rank[Sa[i]] = i;
for (i = 0; i<n; i++) {
if (k) k--;
j = Sa[Rank[i] - 1];
while (s[i + k] == s[j + k]) k++;
Height[Rank[i]] = k;
}
} int ff(int x)
{
if (f[x]==x) return x;
else
return f[x] = ff(f[x]);
} int main() {
//freopen("F:\\rush.txt","r",stdin);
int n;
scanf("%d",&n);
scanf("%s", s);
for (int i = 0;i < n;i++) scanf("%lld",&a[i]); s[n] = 0;
build_Sa(n + 1, MAX_CHAR);//注意调用n+1
getHeight(n); for (int i = 2;i <= n;i++)
{
temp[i-1].first = Height[i];
temp[i-1].second = i;
} for (int i = 1;i <= n;i++)
{
f[i] = i;
cnt[i] = 1;
zd[i] = zx[i] = a[Sa[i]];
} sort(temp+1,temp+1+n-1);
reverse(temp+1,temp+1+n-1); for (int i = 0;i <= n-1;i++)
ans[i] = -1e18-10; int j = 1;
for (int r = temp[1].first;r >= 0;r--)
{
sum[r]=sum[r+1],ans[r] = ans[r+1];
for ( ;j <= n-1 && temp[j].first==r;j++)
{
int x = ff(temp[j].second),y = (temp[j].second-1);
sum[r] += cnt[x]*cnt[y];
ans[r] = max(ans[r],zd[x]*zd[y]);
ans[r] = max(ans[r],zx[x]*zx[y]);
f[y] = x;
cnt[x] += cnt[y];
zd[x] = max(zd[x],zd[y]);
zx[x] = min(zx[x],zx[y]);
}
} for (int i = 0;i <= n-1;i++)
printf("%lld %lld\n",sum[i],sum[i]==0?0:ans[i]);
return 0;
}
【BZOJ 4199】 [Noi2015]品酒大会的更多相关文章
- BZOJ 4199: [Noi2015]品酒大会 [后缀数组 带权并查集]
4199: [Noi2015]品酒大会 UOJ:http://uoj.ac/problem/131 一年一度的“幻影阁夏日品酒大会”隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发“首席品 ...
- bzoj 4199: [Noi2015]品酒大会
Description 一年一度的"幻影阁夏日品酒大会"隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发"首席品酒家"和"首席猎手&quo ...
- 【刷题】BZOJ 4199 [Noi2015]品酒大会
Description 一年一度的"幻影阁夏日品酒大会"隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发"首席品酒家"和"首席猎手&quo ...
- BZOJ.4199.[NOI2015]品酒大会(后缀自动机 树形DP)
BZOJ 洛谷 后缀数组做法. 洛谷上SAM比SA慢...BZOJ SAM却能快近一倍... 只考虑求极长相同子串,即所有后缀之间的LCP. 而后缀的LCP在后缀树的LCA处.同差异这道题,在每个点处 ...
- BZOJ.4199.[NOI2015]品酒大会(后缀数组 单调栈)
BZOJ 洛谷 后缀自动机做法. 洛谷上SAM比SA慢...BZOJ SAM却能快近一倍... 显然只需要考虑极长的相同子串的贡献,然后求后缀和/后缀\(\max\)就可以了. 对于相同子串,我们能想 ...
- BZOJ 4199 [Noi2015]品酒大会:后缀数组 + 并查集
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=4199 题意: 给你一个长度为n的字符串s,和一个长为n的数组v. 对于每个整数r∈[0,n ...
- BZOJ 4199: [Noi2015]品酒大会( 后缀数组 + 并查集 )
求出后缀数组后, 对height排序, 从大到小来处理(r相似必定是0~r-1相似), 并查集维护. 复杂度O(NlogN + Nalpha(N)) ------------------------- ...
- [BZOJ]4199: [Noi2015]品酒大会(后缀数组+笛卡尔树)
Time Limit: 10 Sec Memory Limit: 512 MB Description Input Output Sample Input 10 ponoiiipoi 2 1 4 7 ...
- bzoj 4199: [Noi2015]品酒大会 后缀树
题目大意: 给定一个长为n的字符串,每个下标有一个权\(w_i\),定义下标\(i,j\)是r相似的仅当\(r \leq LCP(suf(i),suf(j))\)且这个相似的权为\(w_i,w_j\) ...
- uoj 131/bzoj 4199 [NOI2015]品酒大会 后缀树+树d
题目大意 见uoj131 分析 题目的提示还是很明显的 \(r\)相似就就代表了\(0...r-1\)相似 建出后缀树我们能dfs算出答案 再后缀和更新一下即可 注意 细节挺多的,但数据很良心 不然我 ...
随机推荐
- 如何使用jquery判断一个元素是否含有一个指定的类(class)
如何使用jquery判断一个元素是否含有一个指定的类(class) 一.总结 一句话总结:可以用hasClass方法(专用)和is方法 1.is(expr|obj|ele|fn)的方法几个参数表示什么 ...
- js数组操作(增、删、改、查)
数组是js中非常常用的一个对象,它有一些经典的操作,今天零度就为大家介绍介绍. 首先,声明一个数组的方法有两种: var arr = []; 或者 var arr = new Array(); 一般情 ...
- K短路 spfa + A*
#include <stdio.h> #include <string.h> #include <queue> #include <algorithm> ...
- Impala储存与分区
不多说,直接上干货! hive的元数据存储在/user/hadoop/warehouse Impala的内部表也在/user/hadoop/warehouse. 那两者怎么区分,看前面的第一列. 下面 ...
- Server.MapPath()的用法
http://blog.csdn.net/qiuhaifeng_csu/article/details/19416407 Server.MapPath(string path)作用是返回与Web服务器 ...
- 洛谷 P1801 黑匣子_NOI导刊2010提高(06)(未完)
P1801 黑匣子_NOI导刊2010提高(06) 题目描述 Black Box是一种原始的数据库.它可以储存一个整数数组,还有一个特别的变量i.最开始的时候Black Box是空的.而i等于0.这个 ...
- [Python] List & Object spread in Python
def myfunc(x, y, z): print(x, y, z) tuple_vec = (, , ) dict_vec = {, , } >>> myfunc(*tuple_ ...
- [Android] 图像各种处理系列文章合集
这是我近期在做Android随手拍一个项目的各种网上关于图片处理的资料,曾经学过数字图像处理都是用C++写的,以下的资料个人觉得是很优秀的各种集合,还有一方面它是在线笔记,希望对大家有所帮助吧 ...
- stm8开发环境配置及测试
需要准备的软件,硬件,IAR for stm8 (EWSTM8).stm8s标准固件库.ST-LINK.STM8s003f3核心板 安装IAR(其中包括st-link的驱动), 到这个网址下载stm8 ...
- Day5下午解题报告1
预计分数:100+60+30=190 实际分数:100+60+30=190 终于有一道无脑T1了哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈 ...