Description

\(n \le 2*10^5\)

给定限制序列 \(A\)

求满足 \(P_i\le A_i\) 的所有排列中

逆序对个数的和

Solution

考虑知道一个 \(A\) 序列时怎么计算排列个数

记 \(C[i]\) 表示 \(A\ge i\) 的个数

然后依次决定 \(n, n-1, \cdots 1\) 填在哪里

填 \(k\) 时有 \(C[k]\) 个可选位置, 其中 \(n-k\) 被占用

因此总排列个数为 \(\prod_{k=1}^n C[k]-(n-k)\)

不难证明无解时该式返回0

用期望线性性去统计答案

考虑 \(i < j, P_i>P_j\)

当 \(A_i\le A_j\) 时, 可知 \(P_j\le A_i\), 把 \(A_j\) 改为 \(A_i\), 此时的所有排列中有一半是满足逆序对的

当 \(A_i\gt A_j\) 时, 转为统计 \(P_i<P_j\) 的方案数, 用总的去减这个就好了

注意到修改时把 \((A_j, A_i]\) 的 \(C\) 都减了一, 定义这个是 \(C'\)

按权值从大到小枚举 \(A_j\)

我们要求的是 \(suff[A_i+1] * C'[A_j+1..A_i] * pref[A_j] * sign(j<i)\)

按标号维护线段树, 第一项维护在线段树里, 第二项通过打标记维护

Code

#include <bits/stdc++.h>
using namespace std;
#define ri rd<int>
#define rep(i, a, b) for (int i = (a), _ = (b); i <= _; ++i)
#define per(i, a, b) for (int i = (a), _ = (b); i >= _; --i)
#define For(i, a, b) for (int i = (a), _ = (b); i < _; ++i)
const int maxN = 2e5 + 7;
typedef long long LL;
const LL O = 1e9 + 7; template<class T> T rd() {
bool f = 1; char c = getchar(); for (; !isdigit(c); c = getchar()) if (c == '-') f = 0;
T x = 0; for (; isdigit(c); c = getchar()) x = x * 10 + c - '0'; return f ? x : -x;
} int n;
struct Node {
int v, i;
inline bool operator < (const Node &y) const {
return v > y.v;
}
}a[maxN];
LL c[maxN], d[maxN];
LL C[maxN], sufC[maxN]; namespace Seg {
const int maxN = ::maxN << 1;
struct Node;
typedef Node* node;
struct Node {
node lc, rc;
int l, mid, r;
int sz;
LL sum, tag;
bool havtag;
Node(int _sz = 0, LL _sum = 0) : sz(_sz), sum(_sum) { havtag = false; } Node operator + (const Node &v) const {
return Node(sz + v.sz, (sum + v.sum) % O);
} void totag(LL d) {
(sum *= d) %= O;
if (!havtag) havtag = true, tag = 1;
(tag *= d) %= O;
} void pushdown() {
if (havtag) {
lc->totag(tag);
rc->totag(tag);
havtag = false;
}
}
}*rt, pool[maxN], *tpool = pool; void build(node &x, int l, int r) {
x = tpool++;
x->l = l, x->r = r, x->mid = (l + r) >> 1;
if (l == r) return;
build(x->lc, l, x->mid);
build(x->rc, x->mid+1, r);
} void init() {
build(rt, 1, n);
} void ins(node x, int to, LL d) {
x->sz++;
(x->sum += d) %= O;
if (x->l == x->r) return;
x->pushdown();
if (to <= x->mid) ins(x->lc, to, d);
else ins(x->rc, to, d);
} void ins(int x, LL d) {
ins(rt, x, d);
} Node get(node x, int l, int r) {
if (l <= x->l && x->r <= r) return *x;
x->pushdown();
if (r <= x->mid) return get(x->lc, l, r);
if (x->mid < l) return get(x->rc, l, r);
return get(x->lc, l, x->mid) + get(x->rc, x->mid+1, r);
} Node get(int l, int r) {
if (l > r) return Node();
return get(rt, l, r);
}
} int main() {
#ifndef ONLINE_JUDGE
freopen("a.in", "r", stdin);
#endif n = ri();
rep (i, 1, n) {
c[a[i].v = ri()]++;
a[i].i = i;
} per (i, n, 1) c[i] += c[i+1];
rep (i, 1, n) {
c[i] = c[i] - (n-i);
d[i] = c[i] - 1;
}
sufC[n] = 1; per (i, n-1, 0) sufC[i] = sufC[i+1] * c[i+1] % O;
C[0] = 1; rep (i, 1, n) C[i] = C[i-1] * c[i] % O; Seg::init();
sort(a+1, a+n+1);
LL ans = 0, cnt = 0;
for (int v = n, i = 1; v; --v) {
LL res = 0;
for (; i <= n && a[i].v == v; ++i) {
Seg::Node tp = Seg::get(1, a[i].i-1);
res -= tp.sum;
cnt += tp.sz;
tp = Seg::get(a[i].i+1, n);
res += tp.sum;
Seg::ins(a[i].i, sufC[v]);
}
ans += res % O * C[v] % O;
Seg::rt->totag(d[v]);
}
ans = ans % O * (O+1) / 2 % O;
ans += cnt % O * C[n] % O;
printf("%lld\n", (ans % O + O) % O); return 0;
}

AGC023E - Inversion的更多相关文章

  1. HDU 1394 Minimum Inversion Number ( 树状数组求逆序数 )

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1394 Minimum Inversion Number                         ...

  2. 控制反转Inversion of Control (IoC) 与 依赖注入Dependency Injection (DI)

    控制反转和依赖注入 控制反转和依赖注入是两个密不可分的方法用来分离你应用程序中的依赖性.控制反转Inversion of Control (IoC) 意味着一个对象不会新创建一个对象并依赖着它来完成工 ...

  3. HDU 1394 Minimum Inversion Number(最小逆序数 线段树)

    Minimum Inversion Number [题目链接]Minimum Inversion Number [题目类型]最小逆序数 线段树 &题意: 求一个数列经过n次变换得到的数列其中的 ...

  4. 依赖倒置原则(Dependency Inversion Principle)

    很多软件工程师都多少在处理 "Bad Design"时有一些痛苦的经历.如果发现这些 "Bad Design" 的始作俑者就是我们自己时,那感觉就更糟糕了.那么 ...

  5. HDU 1394 Minimum Inversion Number(最小逆序数/暴力 线段树 树状数组 归并排序)

    题目链接: 传送门 Minimum Inversion Number Time Limit: 1000MS     Memory Limit: 32768 K Description The inve ...

  6. Inversion Sequence(csu 1555)

    Description For sequence i1, i2, i3, … , iN, we set aj to be the number of members in the sequence w ...

  7. ACM: 强化训练-Inversion Sequence-线段树 or STL·vector

    Inversion Sequence Time Limit:2000MS     Memory Limit:262144KB     64bit IO Format:%lld & %llu D ...

  8. ACM Minimum Inversion Number 解题报告 -线段树

    C - Minimum Inversion Number Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d &a ...

  9. HDU-Minimum Inversion Number(最小逆序数)

    Problem Description The inversion number of a given number sequence a1, a2, ..., an is the number of ...

随机推荐

  1. FreeBSD--如何最有效率的安装软件

    Freebsd 下如何最有效率的安装软件 From <https://www.cnblogs.com/apexchu/p/4131821.html> FreeBSD的默认下载工具是fetc ...

  2. [Codeforces958C2]Encryption (medium)(区间DP)

    Description 题目链接 Solution 显然的区间DP,正常想法f[i][j]表示前i个数分成j块,每次在i前找一个k使得balala,然而常规打法会超时 我们发现,对于i前面的所有点,他 ...

  3. Linux命令、权限

    一.新建用户natasha,uid为1000,gid为555,备注信息为“master”: groupadd -g 555 natasha useradd -u 1000 -g 555 -c mast ...

  4. Android 自定义圆形图片 CircleImageView

    1.效果预览 1.1.布局中写自定义圆形图片的路径即可 1.2.然后看一看图片效果 1.3.原图是这样的 @mipmap/ic_launcher 2.使用过程 2.1.CircleImageView源 ...

  5. Java语言基础---转义符

    转义符 转义符使用“\”表示.常用转义符如下: 1.‘\n’回车 2.‘\t’制表位字符,一个表示向右跳8-10个字符 3.‘\\’表示’\’ 4.‘\’’表示单引号 5.‘\’’’表示双引号 6.‘ ...

  6. Java语言基础---变量与数据类型

    变量的作用域 java用一对大括号“{}”作为语句块的范围,称为作用域.作用域中的变量不能重复定义:离开作用域,变量所分配的内存空间将被JVM所收回. 基本数据类型的包装类 java为基础数据类型提供 ...

  7. MethodTrace 生成的trace文件为空

    今天我准备生成一个trace文件,看看程序卡在哪里. 一般: Debug.startMethodTracing("yuge"); Debug.stopMethodTracing() ...

  8. CSS 一些基础知识(优先级、行内元素的一些属性、font-size单位) 怎样不加载图片

    CSS大小写不敏感 选择器优先级如下所示: 在属性后面使用 !important 会覆盖页面内任何位置定义的元素样式. 作为style属性写在元素内的样式 id选择器 类选择器 标签选择器 通配符选择 ...

  9. Android学习记录(3)—Android中ContentProvider的基本原理学习总结

    一.ContentProvider简介        当应用继承ContentProvider类,并重写该类用于提供数据和存储数据的方法,就可以向其他应用共享其数据.虽然使用其他方法也可以对外共享数据 ...

  10. iOS 中 AFNetworking HTTPS 的使用

    由于我们公司由HTTP转HTTPS,出现了一系列问题特此记录下. 一.HTTPS 二.App Transport Security 三.iOS 中用HTTPS 注意的问题 四.使用 AFNetwork ...