差分&前缀和

前缀和嘛，在算法里是一种比较重要的预处理工具。。。。。。。

好吧，我编不下去了，好像也没有了吧 [doge]

差分的话，可以说是前缀和的逆运算，往往要配合前缀和一起使用，来解决一些区间问题。

最近我也是遇到了好多的差分问题，所以准备总结一下水一篇 。这里的话，我们主要介绍一下一维和二维的前缀和和差分。好了话不多说了。

前缀和;

一维的话，就没什么说的了，

我们来看一下二维：

 我们先来看看某个点 [ 2,3 ] 的前缀和怎么求：

可以看到，如果我们想求，我们就要加上，从[1,1] 到 [ 2,3 ] 全部点的和；

重点就是我们怎么加呢；

我们在计算[ 2,3 ] 前缀和的时候，[ 2,2 ] 和 [ 1,3 ] 就一定是已经算好了的；

那我们可以发现的，[ 2,2 ] 和 [1 , 3] 的前缀和是包含在[2 , 3] 内部的，并且重复包含[ 1,2 ]

于是乎：(s:前缀和数组 , a:原数组)

好，相信大家已经学会了吧；

那我们再来看看前缀和区间怎么求；（比如我们要算红色方块的区间前缀和 [6,4],[4,3]）

我们发现同理，那就剪掉多余的呗；

但是我们发现，剪的过程中，[3,2] 这个多余的部分，被减掉了两次，所以最后我们还要加上这个部分；

于是乎；

;

好了，那我们可以来整题了；

 领地选择 - 洛谷

这题就不用解释了吧【doge】

#include<bits/stdc++.h>
using namespace std;

const int N = 1010;
typedef long long LL; 
int n,m,c;

int a[N][N];

LL s[N][N];
int main()
{
	
	cin>>n>>m>>c;
	for(int i =1 ;i <= n;i ++)
		for(int j =1 ;j <=m ;j ++)
			cin>>a[i][j];
	
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			s[i][j] = s[i-1][j] + s[i][j-1]- s[i-1][j-1] + a[i][j];
	
	LL res = -1e8;
	int x,y;
	for(int i =1 ;i <= n -c + 1 ;i ++)
		for(int j = 1 ;j <= m -c + 1 ;j ++)
		{
			int x1,y1,x2,y2;
			x1= i;
			y1 = j;
			x2 = i + c -1;
			y2 = j + c -1;
			
			LL mx = s[x2][y2]- s[x1-1][y2] - s[x2][y1-1] + s[x1-1][y1-1];
			if(mx > res)
			{
				res = mx;
				x = i,y = j;
			}
		}
	cout<<x<<" "<<y<<endl; 
	return 0;
}

差分;

逆过程嘛；a[i] -= a[i-1];

那我们想要得到原数组的话，只需要 a[i] += a[i-1] (前缀和);

上面说了，差分最常用的用法就是，维护某个区间，使每次操作对该区间的贡献相同；

说人话就是: 给该区间上的数都加上一个值；

 直接看图；

如何我们想给 [3,5] 加上1；

我们考虑前缀和，前缀和中第一个数的值，会对后面所有数都产生一个贡献，那我们在3上加上1，这个1将对后面全部的数都产生一个贡献1，但是我只想要加到5啊，也就是说我们不想让贡献加到6 上，那我们在6的位置上减1 即可，这个贡献也就从6开始不在对后面的数产生贡献了。

 

我们如法炮制;

我们想在 [4,3] -- [4,6] 上加上1;

我们在[4,3] 上 加 1 将会对该点右下角所有的点都作出贡献；

因此我们减掉其他部分的贡献；

我们在[4,5] 上加上 -1 导致灰色框内不在受到[4,3] 加1的贡献；在[7,3] 上加上 -1 导致深红色框内不在受到[4,3] 加1的贡献；

但是，又出现了重叠，我们多减了一部分，所以我们在[7,5] 上加上 1 将多减掉的加回来；

我们在来整个差分的题吧;

海底高铁 - 洛谷

#include<bits/stdc++.h>
using namespace std;

const int N  = 100010;
typedef long long LL;


LL n,m;
queue<int> p;

LL a[N],b[N],c[N];

LL s[N];
LL res;
int main()
{
	scanf("%d%d",&n,&m);
	
	for(int i =1;i <= m;i ++)
	{
		int c;
		scanf("%d",&c);
		p.push(c);
	} 
	
	for(int i =1;i <= n - 1 ;i ++)	scanf("%d%d%d",&a[i],&b[i],&c[i]);
		
	int st,ed;
	st = p.front(),p.pop();
	while(!p.empty())
	{
		ed = p.front(), p.pop();
		
		if(st > ed)
		{
			s[ed] ++;
			s[st] --;	
		} 
		else 
		{
			s[st] ++;
			s[ed] --;
		}
		st = ed;
	}
	
	for(int i =1 ;i <= n - 1 ;i ++) s[i] += s[i - 1];
	
	
	for(int i =1 ;i <= n-1 ; i ++)
	{
		if(s[i] * a[i] > s[i]* b[i] + c[i]) res += s[i]* b[i] + c[i];
		else res += s[i] * a[i];
	}
	
	printf("%lld\n",res);
	return 0;
} 

 收工吧；