[USACO 2010 Nov G]Cow Photographs（思维 + 树状数组/dp）

题目链接

好题！！！md看了半天题解，绝大部分讲的太抽象了…服了…

抽象一下题意：

给定一个n的排列。每次操作可以交换相邻两个数。
问将序列变成一个：i,i+1,i+2,…,n,1,2,…,i−1 形式的序列最少操作次数。

分析

参考题解点这里

首先想一下把序列变成 [1 2 … n] ，即i=1的情况怎么做？
由于每次只能交换相邻的数对，因此答案就是原序列中的逆序对个数；
求逆序对个数，这个简单啊，从前往后扫一遍树状数组维护一下就行了。

然后我们考虑如何从i变成i+1。
每次，我们考虑将当前序列中最小的数变成最大的数加1，再次求这个序列的逆序对个数，就是i+1的答案。
他们之间的关系，就是最小的数后面的数的个数减去前面的数的个数。

具体例子：
考虑变成 [2 3 … n 1] 的情况要怎么做？
原序列 ：[… 1 …]
上个序列：[1 2 3 … n]
当前序列：[2 3 … n 1]
考虑 当前序列 和 上个序列 的答案关系：
由于1是最小数，在上个序列中，它跑到了最前面，对答案产生的贡献就是原序列中1前面的数的个数cnt1；
此时1不在当前序列中的第一个数了，因此我们要撤销这些贡献，即 -cnt1；
由于1是最小数，在当前序列中跑到了最后一个，那么对答案产生的贡献就是原序列中1前面的数的个数cnt2，因此要+cnt2；

以此类推，我们从小到大遍历每个数作为最小数，然后每次算出贡献更新答案；
证明看题解！

代码

#include<iostream>
#include<queue>
#include<cstring>
#include<vector>
#include<stdio.h>
#include<map>
#include<algorithm>
#include<deque>
#include<stack>
#include<set>
// #include <unordered_map>
#include<math.h>
#include<string.h>
#define IOS ios::sync_with_stdio(false),cin.tie(0);
using namespace std;
 
#define pb push_back
#define coutl cout<<"------------"<<endl;
#define fi first
#define se second
  
#define ire(x) scanf("%d",&x)
#define iire(a,b) scanf("%d %d",&a,&b)
#define lre(x) scanf("%lld",&x)
#define llre(a,b) scanf("%lld %lld",&a,&b)
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define endl "\n"
#define PI acos(-1.0)

typedef long long ll;
typedef unsigned long long ull;
      
typedef pair<int, int> PII;
typedef pair<double, int> PDI;
typedef pair<ll, ll> PLL;
typedef pair<double, double> PDD;
typedef pair<double, pair<int, double> > PDID;
typedef pair<char, char> PCC;
typedef pair<char, pair<int, int> > PCII;
typedef pair<int, pair<int, int> > PIII;
typedef pair<int, pair<int, pair<int, int> > > PIIII;
typedef pair<ll, pair<int, int> > PLII;
 
const int maxn = 1e6 + 7;
const int N = 2010 + 7;
const int M = 1e6 + 7;
const int mod = 998244353;
const int inf = 0x3f3f3f3f;
// const ll INF = 0x3f3f3f3f3f3f3f3f;
const double pi = acos(-1);
const double eps = 1e-8;
  
ll gcd(ll a,ll b) {return b==0 ? a : gcd(b,a%b);}
ll lcm(ll a,ll b) {return a*b / gcd(a,b);}
ll qmi(ll a,ll b,ll p) {ll ans = 1; while(b) { if(b & 1) ans = ans * a % p; a = a * a % p; b >>= 1; } return ans;}
int lowbit(int x) {return x & (-x);}

int a[maxn];
int idx[maxn];
int tr[maxn];
int n;

void update(int x,int v)
{
	while(x <= n)
	{
		tr[x] += v;
		x += lowbit(x);
	}
}

int query(int x)
{
	int ans = 0;
	while(x)
	{
		ans += tr[x];
		x -= lowbit(x);
	}
	return ans;
}

int main()
{
	ire(n);
	
	ll last = 0;
	for(int i=1;i<=n;i++)
	{
		ire(a[i]);
		idx[a[i]] = i;
		
		last += query(n)-query(a[i]);
		update(a[i],1);
	}
	
	ll ans = last;
	for(int i=1;i<n;i++)
	{
		last = last - (idx[i]-1) + (n-idx[i]);
		
		ans = min(ans,last);
	}
	
	cout<<ans<<endl;
	
	return 0;
}