• [NOI2015]品酒大会


    原题意:有多少个串的(lcp(i, j) >= r)(其中r为1~n中每一个数)

    我们先不管第二问,只看第一问

    第一次转化

    首先不难发现一个非常好的性质:对于一个r相似的两杯酒,他们肯定也是r-1相似,r-2相似……

    于是,我们考录倒序枚举,於是问题转化成了:有多少个串的(lcp(i, j) == r)(其中r为1~n中每一个数)

    第二次转化

    看到lcp,我们不难想到后缀数组,(lcp(i, j) == r) 又等价于(min(he[i + 1], he[i + 2], ……, he[j]))

    于是我们把问题重新考虑,就转化成了:求出he数组,然后问有多少个数对,满足(i - j)的he的最小值恰好等于r

    第三次转化

    我们把he数组降序排一遍序,然后按照顺序插入

    然后我们可以把问题理解成:对于每次插入,我可以连接连续两堆数,然后这两堆数中每一个数值都比r大,所以我们可以把第一堆数的个数和第二堆数的个数相乘,得到的就是这两堆数对答案的贡献

    举个例子:

    假设我现在的数列长这样(还未被插入的数(即比r小的数)为*):(5 4 * 4 5 6 7 * * 4 * 8)

    然后我现在要在第三个位置插入一个3

    不难发现,一共有((2 + 1) * (4 + 1) - 1 = 14)种方案(+1是因为3本身也可以算进来, -1是因为([3, 3])不能算)

    所以这个合并的过程类似于并查集,我们只需要维护一个size就可以求出第一问了

    对于第二问,我们只需要维护一个最小值和最大值即可(维护最小值是因为有复数)

    代码如下

    #include<bits/stdc++.h>
    using namespace std;
    #define il inline
    #define re register
    #define debug printf("Now is Line : %d
    ",__LINE__)
    #define file(a) freopen(#a".in","r",stdin);freopen(#a".out","w",stdout)
    #define int long long
    #define inf 1234567890000000000
    il int read() {
        re int x = 0, f = 1; re char c = getchar();
        while(c < '0' || c > '9') { if(c == '-') f = -1; c = getchar();}
        while(c >= '0' && c <= '9') x = x * 10 + c - 48, c = getchar();
        return x * f;
    }
    #define rep(i, s, t) for(re int i = s; i <= t; ++ i)
    #define drep(i, s, t) for(re int i = t; i >= s; -- i)
    #define _ 300005
    int n, m, sa[_], rk[_], tp[_], a[_], he[_], val[_], ans1[_], ans2[_], now, num, b[_];
    char c[_];
    vector<int> q[_];
    struct set {int size, mi, ma, fa;}e[_];
    il void Qsort() {
    	rep(i, 1, m) a[i] = 0;
    	rep(i, 1, n) ++ a[rk[i]];
    	rep(i, 1, m) a[i] += a[i - 1];
    	drep(i, 1, n) sa[a[rk[tp[i]]] --] = tp[i];
    }
    il void get_sort() {
    	for(re int w = 1, p = 0; p < n && w <= n; m = p, p = 0, w <<= 1) {
    		rep(i, n - w + 1, n) tp[++ p] = i;
    		rep(i, 1, n) if(sa[i] > w) tp[++ p] = sa[i] - w;
    		Qsort(), swap(tp, rk), rk[sa[1]] = p = 1;
    		rep(i, 2, n) rk[sa[i]] = (tp[sa[i]] == tp[sa[i - 1]] && tp[sa[i] + w] == tp[sa[i - 1] + w]) ? p : ++ p;
    	}
    }
    il void get_height() {
    	int j = 0;
    	rep(i, 1, n) {
    		if(rk[i] == 1) continue;
    		if(j) -- j;
    		while(c[i + j] == c[sa[rk[i] - 1] + j]) ++ j;
    		he[rk[i]] = j;
    	}
    }
    il int find(int x) {
    	while(x != e[x].fa) x = e[x].fa = e[e[x].fa].fa;
    	return x;
    }
    il void merge(int a, int b) {
    	int x = find(a), y = find(b);
    	now += e[x].size * e[y].size, num = max(num, max(e[x].ma * e[y].ma, e[x].mi * e[y].mi));
    	e[y].fa = x, e[x].size += e[y].size, e[x].ma = max(e[x].ma, e[y].ma), e[x].mi = min(e[x].mi, e[y].mi);
    }
    signed main() {
    	n = read(), scanf("%s", c + 1), m = 26, num = -inf;
    	rep(i, 1, n) b[i] = read(), rk[i] = c[i] - 'a' + 1, tp[i] = i;
    	Qsort(), get_sort(), get_height();
    	rep(i, 1, n) e[i].fa = i, e[i].size = 1, e[i].ma = e[i].mi = b[sa[i]], q[he[i]].push_back(i);
    	drep(i, 0, n - 1) {
    		int pax = q[i].size();
    		rep(j, 0, pax - 1) merge(q[i][j], q[i][j] - 1);
    		if(now) ans1[i] = now, ans2[i] = num;
    	}
    	rep(i, 0, n - 1) printf("%lld %lld
    ", ans1[i], ans2[i]);
    	return 0;
    }
    
  • 相关阅读:
    DateUtils
    Java静态绑定与动态绑定
    Mysql中实现递归查询
    架构一、核心概念
    Spring cron 表达式
    MySql点点滴滴(一)之可视化工具介绍
    java中注解的使用与实例
    3、第一个Python程序
    CSS命名
    如何在Notepad++ 中成功地安装Emmet 插件
  • 原文地址:https://www.cnblogs.com/bcoier/p/11441727.html
Copyright © 2020-2023  润新知