Why Did the Cow Cross the Road I P

W h y   D i d   t h e   C o w   C r o s s   t h e   R o a d   I   P Why\ Did\ the\ Cow\ Cross\ the\ Road\ I\ P Why Did the Cow Cross the Road I P

题目链接: l u o g u   P 3656 luogu\ P3656 luogu P3656 \ j z o j   6684 jzoj\ 6684 jzoj 6684

题目

为什么牛要过马路?我们可能永远都不知道全部原因,但肯定的是,农民约翰的牛经常穿过马路。事实上,他们经常过马路,当他们的路径交叉时,他们经常碰到对方,这是约翰希望改进的情况。

农夫约翰将奶牛分成了 N N NN个不同的品种 ( 1 ≤ N ≤ 100000 ) (1\le N\le 100000) 1N100000,他的每一个田地都是专门为一种特定的品种放牧; 例如,专门用于品种 12 12 12奶牛的田地只能用于品种 12 12 12的奶牛,而不能用于任何其他品种的奶牛。一条长长的路穿过他的农场。道路一侧有 N N N个田地(每个品种一个),道路另一侧也有 N N N个田地(每个品种也有一个)。所以,当一头牛穿过马路,它会穿过特定品种的两个田地。

如果农民约翰更仔细的研究他的计划,他会在道路两旁订购同品种田地,所以每个品种的两个田野将直接穿过彼此的道路。这样就可以让奶牛随意过马路,而不会有不同品种的牛碰撞在一起。唉,但现在道路两边的田地可能不一样,所以农民约翰明白,可能会有一对品种 ( a , b ) (a,b) (a,b)是交叉的当且仅当是一对不同的品种,且道路上一边田地品种 a a a上的任何路线都与道路另一边田地品种 b b b的任何路线相交 。

农民约翰希望尽量减少品种的交叉对数。由于物流方面的原因,他表示,他可以在道路的一边移动牛,所以那边的田地可以“循环转移”。也就是说,对于一个给定的常值 k k k来说,循环转移是指每头牧场里的牛都会向前移动 k k k个田野,而最前面的 k k k只奶牛会移动到最后,因为向前移动的奶牛已经填补了前 k k k个田地。比如,如果道路一侧的田地移动前的排列为 3 , 7 , 1 , 2 , 5 , 4 , 6 3,7,1,2,5,4,6 3,7,1,2,5,4,6,假设 k k k 2 2 2,新排列将为 4 , 6 , 3 , 7 , 1 , 2 , 5 4,6,3,7 ,1,2,5 4,6,3,7,1,2,5.请编程确定在道路一侧的田野适当循环移动后存在的品种交叉对的最小值。

输入

一行包含 N ( 1 ≤ N ≤ 100 , 000 ) N (1 \leq N \leq 100,000) N(1N100,000),接下来 N N N行按顺序描述了小路一旁田地的品种的 I D ID ID号,每一个 I D ID ID号是一个在 1... N 1...N 1...N之间的整数。倒数 N N N行描述了小路另一旁田地的品种的 I D ID ID号。

输出

循环移动道路侧的田野后,请输出可能的品种交叉对的最小数量(两边均可移动)

样例输入

5
5
4
1
3
2
1
3
2
5
4

样例输出

0

思路

这道题要用树状数组。

我们先让 p l [ b [ i ] ] = i pl[b[i]] = i pl[b[i]]=i p a r [ i ] = p l [ a [ i ] ] par[i] = pl[a[i]] par[i]=pl[a[i]],那么 p a r [ i ] par[i] par[i]就代表 a [ i ] a[i] a[i]在序列 b b b中的位置。
那么我们就可以通过树状数组求出 p a r [ i ] par[i] par[i]的逆序对,而它的逆序对就是当前的交叉对。
那我们只要不断的把最后面的数移到最前面,然后维护逆序对个数,就可以了。

至于怎么维护,就是把原来的逆序对 − ( n − p a r [ i ] ) + ( p a r [ i ] − 1 ) - (n - par[i]) + (par[i] - 1) (npar[i])+(par[i]1)
(一个是把后面拿掉所扣的,另一个是放在前面所加的)

不过还有一个就是,因为 a a a b b b两个序列都可以动,所以我们要分别试让 a a a序列动和让 b b b序列动的情况。
(不过好像也不用,但我没试,不知道能不能只动一个)

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#define ll long long

using namespace std;

int n, a[100001], b[100001], pl[100001], par[100001], ni[100001];
ll ans;

int get(int now) {//单点求值
	int answer = 0;
	for (; now <= n; now += now & -now)
		answer += ni[now];
	return answer;
}

void build(int now, int add) {//单点加值
	for (; now; now -= now & -now)
		ni[now] += add;
}

ll work() {
	memset(ni, 0, sizeof(ni));//初始化
	ll an = 0;
	for (int i = 1; i <= n; i++) {//求出当前逆序对个数
		an += get(par[i]);
		build(par[i], 1);
	}
	ll getnew = an;
	for (int i = n; i >= 0; i--) {//移数
		getnew = getnew - n + par[i] * 2 - 1;//得到新的逆序对的数量
		an = min(getnew, an);//找到逆序对最少的情况
	}
	return an;
}

int main() {
//	freopen("mincross.in", "r", stdin);
//	freopen("mincross.out", "w", stdout);
	
	scanf("%d", &n);//读入
	for (int i = 1; i <= n; i++) scanf("%d", &a[i]);//读入
	for (int i = 1; i <= n; i++) scanf("%d", &b[i]);//读入
	
	for (int i = 1; i <= n; i++) pl[b[i]] = i;
	for (int i = 1; i <= n; i++) par[i] = pl[a[i]];//初始化
	ans = work();
	
	for (int i = 1; i <= n; i++) pl[a[i]] = i;
	for (int i = 1; i <= n; i++) par[i] = pl[b[i]];//反过来也要一次(因为b也可以动)
	ans = min(ans, work());
	
	printf("%lld", ans);//输出
	
//	fclose(stdin);
//	fclose(stdout);
	
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值