poj 3263 Tallest Cow ， 差分学习日记

差分学习日记：（来自算法竞赛进阶指南）

差分的定义：对于一个数列A， 其差分数列B满足： B【1】=A【1】 ， B【n】=A【n】-A【n-1】 。差分经常用于区间操作，跟前缀和一样，主要用来减少计算量，将区间操作转化为两个端点的操作。

比如说： 给一个长度为n的数列，每次可以选择一个区间 【l，r】，并且将这个区间内的数字+1/-1 ， 问至少需要操作多少次才能使得数列中所有数都一样，此时数列共有几种不同的情况。

那么此题就可以转化为差分来做，造一个差分数列，使得b1=a1，bn=an-an-1， bn+1=0.

那么每一次操作就成了数列b里面挑选两个数使得一个数+1，另外一个数-1，使得数列b2.b3 …… bn变为0。 

至于为什么b1=a1， 结束之后b1的值就是数列的值。

对于中间的数，我们需要找到所有一正一负的两个数，一个-1，一个+1

假如全部都是同号的话，再找中间的数一个+1，一个-1就纯属浪费操作，因此，我们需要找开头结尾，b1/bn+1 和 bi （ 2<=i<=n) 进行操作

设b2~bn的正数总和为p，负数总和为q，则需要进行 min(p,q)+abs(p-q)=max(p,q)次 ， 变化过程中一共可以产生p-q+1种不同的b1值。

 

回到poj 3263

题目：http://poj.org/problem?id=3263

做法：

其实这道题不难想到一开始全部置为最高，然后根据出现的l,r可以互相看得见，将l，r中间的数都-1 ， 注意l，r只能出现过一次，重复出现不操作。此时便是最优解。

这道题数据不强，暴力都可以过，下面我就给出暴力代码和差分代码。

暴力：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
const int maxn = 1e4 + 7;
typedef pair<int, int> P;
map<P, bool> M;
int num[maxn];
int main() {
	int n, i, h, r;
	scanf("%d %d %d %d", &n, &i, &h, &r);
	fill(num + 1, num + 1 + n, h);
	for (int i = 1; i <= r; i++) {
		int fro, to;
		scanf("%d %d", &fro, &to);

		if (fro > to) swap(fro, to);
		if (M.count(P(fro, to)) != 0) continue;
		for (int j = fro + 1; j < to; j++) {
			num[j]--;
		}
		M[P(fro, to)] = 1;
	}
	for (int i = 1; i <= n; i++) printf("%d\n", num[i]);
	return 0;
}

 

差分：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
typedef pair<int, int> P;
const int maxn = 1e4 + 7;
map<P, bool> M;
int d[maxn];
int main() {
	int n, i, h, r;
	cin >> n >> i >> h >> r;
	for (int i = 1; i <= r; i++) {
		int fro, to;
		scanf("%d %d", &fro, &to);
		if (M.count(P(fro, to)) != 0) continue;
		if (fro > to) swap(fro, to);
		//if (abs(fro - to) <= 1) continue;
		d[fro + 1]--;
		d[to]++;
		P p;
		p.first = fro;
		p.second = to;
		M[p] = 1;
	}
	int flag = 0;
	for (int i = 1; i <= n; i++) {
		flag += d[i];
		printf("%d\n", h + flag);
	}
	return 0;
}