CodeForces 351B Jeff and Furik 概率DP 逆序对

题目大意：

就是现在给出一个1~n的排列, Jeff和Furik分别轮流进行操作, Jeff先手, Jeff会选择相邻的两个数p[i], p[i + 1]交换位置, 然后轮到Furik, Furiki每次都会抛一个硬币, 出现正面就在序列中选取相邻的满足p[i] > p[i + 1]的两个数交换, 出现反面则选取任意一个相邻的满足p[i] < p[i + 1]的一对数进行交换, 操作时当这个序列变成递增序列的时候操作结束, 假设Jeff每一步都操作都足够好(使得接下来剩余的操作次数最少)那么问一共需要的操作步数的期望是多少

大致思路：

首先Jeff操作要足够好的话就会在每一次操作的时候尽量的交换p[i] > p[i + 1]的数的位置, 这样剩余的期望步数才会最少, 那么不难找出递推关系

具体递推关系见代码注释

代码如下：

Result  :  Happy New Year!     Memory  :  100 KB     Time  :  62 ms

/*
 * Author: Gatevin
 * Created Time:  2015/1/6 19:44:21
 * File Name: Kotomi.cpp
 */
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#include<iomanip>
using namespace std;
const double eps(1e-8);
typedef long long lint;

#define lson l, mid, rt << 1
#define rson mid + 1, r, rt << 1 | 1

/*
 * 首先这个问题涉及到逆序数的概念
 * 在一个排列中，如果一对数的前后位置与大小顺序相反，即前面的数大于后面的数，
 * 那么它们就称为一个逆序。一个排列中逆序的总数就称为这个排列的逆序数。
 * 那么我们视一个一个升序序列为逆序数为0的序列的话
 * 用E[i]表示当前序列的逆序数为i时到达目标状态的概率, 有E[0] = 0; E[1] = 1(Jeff先手);
 * 有余每次操作时, Jeff都会使得逆序数减少1, 而Furik会有0.5的可能性使得逆序数+1, 0.5的可能性-1
 * 则E[i] = 1 + 1 + E[i - 2]*0.5 + E[i - 1 + 1]*0.5 (i >= 2)
 * 化简得E[i] = E[i - 2] + 4
 * 所以E[i] = 2*i (i是偶数)
 *     E[i] = 2*i - 1 (i是奇数)
 */

/*
 * 那么问题变为求一个长度为n的排列中的逆序对的个数
 * 本题 n <= 3000可以采用O(n^2)的暴力求解
 * 不过还是采用线段树来求逆序对的个数了, 复杂度O(nlogn)算是了解一下逆序对的线段树写法
 * 线段树中rt节点表示l到r中被标记出现的数的个数, 每次插入一个数更新节点值并查询 [p[i] + 1, n] 之间出现过的数的个数
 */


int sum[3000 << 2];
int n, p;

void pushUp(int rt)
{
    sum[rt] = sum[rt << 1] + sum[rt << 1 | 1];
}

void build(int l, int r, int rt)
{
    if(l == r)
    {
        sum[rt] = 0;
        return;
    }
    int mid = (l + r) >> 1;
    build(lson);
    build(rson);
    pushUp(rt);
    return;
}

void update(int l, int r, int rt, int num)
{
    if(l == r)
    {
        sum[rt] = 1;
        return;
    }
    int mid = (l + r) >> 1;
    if(mid >= num) update(lson, num);
    if(mid + 1 <= num) update(rson, num);
    pushUp(rt);
    return;
}

int query(int l, int r, int rt, int L, int R)
{
    if(l >= L && r <= R)
        return sum[rt];
    int mid = (l + r) >> 1;
    int ret = 0;
    if(mid >= L) ret += query(lson, L, R);
    if(mid + 1 <= R) ret += query(rson, L, R);
    return ret;
}

int main()
{
    scanf("%d", &n);
    build(1, n, 1);
    int ans = 0;
    for(int i = 1; i <= n; i++)
    {
        scanf("%d", &p);
        update(1, n, 1, p);
        ans += query(1, n, 1, p + 1, n);
    }
    printf("%.6f\n", 2.*ans - (ans & 1));
    return 0;
}