【codeforces 749E】 Inversions After Shuffle
http://codeforces.com/problemset/problem/749/E (题目链接)
题意
给出一个1~n的排列,从中等概率的选取一个连续段,设其长度为l。对连续段重新进行等概率的全排列,求排列后整个原序列的逆序对的期望个数。
Solution
考虑对于每一对数${(a_i,a_j),i<j}$算贡献。
1.连续段包含${a_i,a_j}$
不妨设${a_i<a_j}$,则只有当排列后${a_j}$再${a_i}$前面才会对答案有贡献(${a_i>a_j}$的情况同理),连续段长度为${l}$。
于是满足${a_i}$在${a_j}$前面的排列数为${P_l^{l-2}}$,概率:${\frac{P_l^{l-2}}{P_l^l}=\frac{1}{2}}$。
满足包含${a_i}$和${a_j}$的连续段有${i*(n-j+1)}$个,其概率为:${\frac{2*i*(n-j+1)}{n*(n+1)}}$。
所以其期望等于两个概率相乘:
$${q_{i,j}=\frac{i*(n-j+1)}{n*(n+1)}}$$
这是${O(n^2)}$的,考虑优化。总期望:
$${Q=\sum_{i=1}^n \sum_{j=i+1}^n q_{i,j}}$$
$${Q=\sum_{i=1}^n \sum_{j=i+1}^n \frac{i*(n-j+1)}{n*(n+1)}}$$
发现${(n-j+1)}$是连续的,于是就变成了:
$${Q=\sum_{i=1}^n \frac {i*(n-i)*(n-i+1)} {2*n*(n+1)}}$$
这样复杂度就是${O(n)}$的了。
2.连续段不同时包含${a_i,a_j}$
如果${a_i<a_j}$,那么因为不被连续段同时包含,它们不会有机会改变相对位置,所以不会对答案做出贡献。考虑${a_i>a_j}$的情况。
那么连续段可能取到的区间有:${[1,j-1],[i+1,n]}$。考虑到区间${[i+1,j-1]}$被算了2次,容斥一下,所以区间的概率:
$${P_{i,j}=\frac {(j-1)*j+(n-i)*(n-i+1)-(j-i-1)*(j-i)} {n*(n+1)}}$$
$${P_{i,j}=\frac {(n^2+n)-(2*i+2*n*i)+2*i*j} {n*(n+1)}}$$
这个${P_{i,j}}$怎么快速求解呢,考虑逆序对这个东西。
$${Q=\sum_{i=1}^n \sum_{j=i+1}^n \frac {(n^2+n)-(2*i+2*n*i)+2*i*j} {n*(n+1)}}$$
设满足${a_j<a_i,j>i}$的${a_j}$的个数为${x}$,显然${x}$我们可以通过树状数组用求逆序对的方法${O(nlogn)}$的求出来,则:
$${Q=\sum_{i=1}^n \frac {x*((n^2+n)-(2*i+2*n*i)) + \sum_{j=i+1}^n 2*i*j} {n*(n+1)}}$$
那么现在${\sum_{j=i+1}^n 2*i*j}$怎么求呢。把${2*i}$提出去,那么就成了${2*i*\sum_{j=i+1}^n j}$我们用${y}$记录满足${a_j<a_i,j>i}$的${a_j}$的位置的和,也就是${\sum_{j=i+1}^n j}$,那么显然这个东西我们也是可以通过树状数组用求逆序对的方法${O(nlogn)}$的算出来的。则:
$${Q=\sum_{i=1}^n \frac {x*((n^2+n)-(2*i+2*n*i)) + 2*i*y} {n*(n+1)}}$$
于是问题就${O(nlogn)}$的解决了。
细节
mdzz不晓得哪里爆掉了还是精度问题,调了2天,最后莫名AC。。。
代码
// codeforces 749E
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define LL long long
#define inf 1<<30
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std; const int maxn=100010;
LL c[maxn],s[maxn],n;
int a[maxn];
long double ans; int lowbit(int x) {
return x&-x;
}
void add(LL *c,int x,LL val) {
for (int i=x;i<=n;i+=lowbit(i)) c[i]+=val;
}
LL query(LL *c,int x) {
LL res=0;
for (int i=x;i;i-=lowbit(i)) res+=c[i];
return res;
} void solve1() { //区间包含
long double Q=0;
for (LL i=1;i<=n;i++)
Q+=(long double)(i*(n-i)*(n-i+1))/2/n/(n+1);
ans+=Q;
}
void solve2() { //区间不包含
long double Q=0;
for (int i=n;i>=1;i--) {
LL x=query(c,a[i]-1);
Q-=(long double)(x*((2*i+2*n*i)-(n*n+n)))/n/(n+1);
Q+=(long double)(2*i)/n/(n+1)*query(s,a[i]-1);
add(c,a[i],1);
add(s,a[i],i);
}
ans+=Q;
}
int main() {
scanf("%lld",&n);
for (int i=1;i<=n;i++) scanf("%d",&a[i]);
solve1();
solve2();
printf("%.20Lf",ans);
return 0;
}
贴一个暴力,供参考:
// codeforces 749E
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define LL long long
#define inf 1<<30
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std; const int maxn=100010;
LL c[maxn],s[maxn],n;
int a[maxn];
long double ans; int main() {
freopen("aaa.in","r",stdin);freopen("ccc.out","w",stdout);
scanf("%lld",&n);
for (int i=1;i<=n;i++) scanf("%d",&a[i]);
for (LL i=1;i<=n;i++)
ans+=(long double)(i*(n-i)*(n-i+1))/(2*n*(n+1));
long double res=0;
for (LL i=n;i>=1;i--) {
for (LL j=i+1;j<=n;j++)
if (a[i]>a[j]) res+=(long double)((j-1)*j+(n-i)*(n-i+1)-(j-i-1)*(j-i))/(n*(n+1));
}
ans+=res;
printf("%.20Lf",ans);
return 0;
}
【codeforces 749E】 Inversions After Shuffle的更多相关文章
- 【codeforces 415D】Mashmokh and ACM(普通dp)
[codeforces 415D]Mashmokh and ACM 题意:美丽数列定义:对于数列中的每一个i都满足:arr[i+1]%arr[i]==0 输入n,k(1<=n,k<=200 ...
- 【codeforces 707E】Garlands
[题目链接]:http://codeforces.com/contest/707/problem/E [题意] 给你一个n*m的方阵; 里面有k个联通块; 这k个联通块,每个连通块里面都是灯; 给你q ...
- 【codeforces 707C】Pythagorean Triples
[题目链接]:http://codeforces.com/contest/707/problem/C [题意] 给你一个数字n; 问你这个数字是不是某个三角形的一条边; 如果是让你输出另外两条边的大小 ...
- 【codeforces 709D】Recover the String
[题目链接]:http://codeforces.com/problemset/problem/709/D [题意] 给你一个序列; 给出01子列和10子列和00子列以及11子列的个数; 然后让你输出 ...
- 【codeforces 709B】Checkpoints
[题目链接]:http://codeforces.com/contest/709/problem/B [题意] 让你从起点开始走过n-1个点(至少n-1个) 问你最少走多远; [题解] 肯定不多走啊; ...
- 【codeforces 709C】Letters Cyclic Shift
[题目链接]:http://codeforces.com/contest/709/problem/C [题意] 让你改变一个字符串的子集(连续的一段); ->这一段的每个字符的字母都变成之前的一 ...
- 【Codeforces 429D】 Tricky Function
[题目链接] http://codeforces.com/problemset/problem/429/D [算法] 令Si = A1 + A2 + ... + Ai(A的前缀和) 则g(i,j) = ...
- 【Codeforces 670C】 Cinema
[题目链接] http://codeforces.com/contest/670/problem/C [算法] 离散化 [代码] #include<bits/stdc++.h> using ...
- 【codeforces 515D】Drazil and Tiles
[题目链接]:http://codeforces.com/contest/515/problem/D [题意] 给你一个n*m的格子; 然后让你用1*2的长方形去填格子的空缺; 如果有填满的方案且方案 ...
随机推荐
- Qt——组件位置随窗口变化
当我们用Qt Designer设计界面时,有时会面临这样一个问题:需要在窗口指定位置放置组件,并且当窗口位置大小改变时,该组件相对其父对象的位置是不变的,如下面两幅图所示 ,首先看上面这幅图,注意bu ...
- 乱码引起的CSS失效原理,解决技巧。
由于一个中文是两个字符组成,在编码不一致的情况下会引发字符的“重新”组合,(半个汉字的编码字符与后面的字符组合生成新的“文字”)引发原本的结束符合“变异”,从而导致找不到结束符号,使得后面的CSS就会 ...
- SharePoint 2013 入门教程
以下文章是自己在学习SharePoint的过程中,不断积累和总结的博文,现在总结一个目录,分享给大家.这个博客也是自己从SharePoint入门,到一个SharePoint开发的成长记录,里面记录的都 ...
- 安装cocoapods遇到两大坑-Ruby版本升级和Podfile的配置
今天安装cocoapods #移除原有ruby源 $ gem sources --remove https://rubygems.org/ #使用可用的淘宝网 $ gem sources -a htt ...
- Android 面试题--Activity
1.什么是 Activity?Activity是Android组件中最基本也是最为常见用的四大组件(Activity,Service服务,Content Provider内容提供,BroadcastR ...
- Android中webView和网页的交互
Android中webView和网页的交互 Android中webView跟网页的交互式通过JavaScript进行的.具体步骤: 1.创建JavaScript,在点击的时候调用JavaScript ...
- tuple放入dict中
tuple放入dict中是否可以正常运行 # 将tuple放入dict中 a = ('AI','Kobe','Yao') b = ('AI',['Kobe','Yao']) dict1 = {'a': ...
- SQL Server 2008 R2:快速清除日志文件的方法
本例,快速清理“students”数据库的日志,清理后日志文件不足1M. USE [master] GO ALTER DATABASE students SET RECOVERY SIMPLE WIT ...
- uboot的配置流程分析
简单介绍一下uboot的基本配置流程.和绝大多数源码编译安装一样,uboot在执行make之前需要执行make XXXconfig来配置相关信息,而且uboot本身是针对多种平台的bootloader ...
- HDU 1856 Brave Game(巴什博奕)
十年前读大学的时候,中国每年都要从国外引进一些电影大片,其中有一部电影就叫<勇敢者的游戏>(英文名称:Zathura),一直到现在,我依然对于电影中的部分电脑特技印象深刻. 今天,大家选择 ...