(洛谷官方题单)前缀和与差分

1. P1115 最大子段和
解题思路：一维最大前缀和的模板题。
对第 i 个数来说，如果前面的 i-1个数的和小于 0，那么加上前面 i-1 个数构成的前缀和不如只有第 i 个数，所以如果前 i-1 个数的和小于 0 则舍弃。

#include <iostream>
using namespace std;
int ma=-99999,sum=0;
int main(){
	int n,now;
	cin >> n;
	while(n--){
		cin >> now;
		sum+=now;
		if(ma<sum){
			ma=sum;
		}
		if(sum<0){
			sum=0;
		}
	}
	cout << ma;
	return 0;
} 

2. P1719 最大加权矩形
解题思路：看完一维的我们来看二维的。
二维矩阵压缩成一维的即可，设矩阵为3*3的一个方阵，1-3行分别标记为a，b，c，此时1-3行可以组合压缩为 a,ab,abc,b,bc,c 三种情况
对 a 这一行我们求一维最大前缀和即可
对 ab 这两行我们先将每一列对应的元素加在一起：a1+b1=n1,a2+b2=n2,a3+b3=n3, 此时将ab两行压缩成一行了，求一维最大前缀和即可。
其它情况以此类推，以下为AC代码：

#include <iostream>
#include <algorithm>
using namespace std;
const int maxn=121;
int mp[maxn][maxn],now[maxn],n,nmax;
int pd(){
	int nsum=0,sum=now[1];
	for(int i=1;i<=n;i++){
		nsum+=now[i];
		if(nsum<0){
			nsum=0;
		}
		if(nsum>sum){
			sum=nsum;
		}
	}
	return sum;
}
int main(){
	cin >> n;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			cin >> mp[i][j];
		}
	}
	nmax=mp[1][1];
	for(int i=1;i<=n;i++){
		for(int j=i;j<=n;j++){
			if(i==j){
				for(int k=1;k<=n;k++){
					now[k]=mp[i][k];
				}
			}else{
				for(int k=1;k<=n;k++){
					now[k]+=mp[j][k];
				}
			}
			nmax=max(nmax,pd());
		}
	}
	cout << nmax;
	return 0;
}

3. P2004 领地选择
解题思路：虽然是前缀和的题目，我们还是祭出我们的DP大法。左上角 (1,1)，右下角 (X,Y) 的矩形的带权矩阵和的值为 dp[x,y]=dp[x-1][y]+dp[x][y-1]-dp[x-1][y-1]+val[x][y]; 左上角 (x1,y1), 右下角 (x2,y2) 的矩阵的带权矩阵和为 dp[x2][y2]+dp[x1-1][y1-1]-dp[x1-1][y2]-dp[x2][y1-1]. 最后遍历一遍左上角为 (i,j) 找最大值即可。AC代码如下：

#include <iostream>
using namespace std;
const int maxn=1001;
typedef long long ll;
long long mp[maxn][maxn],dp[maxn][maxn]={0};
ll nmax=0,nsum;
int main(){
	int n,m,c,x,y;
	cin >> n >> m >> c;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			cin >> mp[i][j];
		}
	}
	for(int i=1;i<=c;i++){
		for(int j=1;j<=c;j++){
			nmax+=mp[i][j];
		}
	}
	x=1;
	y=1;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			dp[i][j]=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1]+mp[i][j];
		}
	}
	for(int i=1;i<=n-c+1;i++){
		for(int j=1;j<=m-c;j++){
			nsum=dp[i+c-1][j+c-1]+dp[i-1][j-1]-dp[i-1][j+c-1]-dp[i+c-1][j-1];
			if(nsum>nmax){
				nmax=nsum;
				x=i;
				y=j;
			}
		}
	}
	cout << x << " " << y;
	return 0;
}

4. P5638 【CSGRound2】光骓者的荣耀
解题思路：DP（没想到吧，还是DP）
首先确定状态 dp[i][0] 表示到城市 i 且不使用传送的时候的最短距离，dp[i][1] 表示到城市 i 且传送依次的最短距离。
状态转移方程：

if(i>k){
	dp[i][1]=min(dp[i-k][0],dp[i-1][1]+val[i]);
}else{
	dp[i][0]=dp[i-1][0]+val[i];
}

AC代码如下:

#include <iostream>
#include <algorithm>
using namespace std;
const int maxn=1000010;
typedef long long ll;
ll a[maxn],dp[maxn][2];
int main(){
	int n,k;
	cin >> n >> k;
	for(int i=1;i<=n-1;i++){
		cin >> a[i];
	}
	for(int i=2;i<=n;i++){
		if(i-k>0){
			dp[i][1]=dp[i-k][0];
			dp[i][1]=min(dp[i-1][1]+a[i-1],dp[i][1]);
		}
		dp[i][0]=dp[i-1][0]+a[i-1];
	}
	cout << dp[n][1];
	return 0;
} 

5. P3406 海底高铁
解题思路：差分，这个知识点很重要。
这道题里面如果我们从城市 1 到城市 3 则需要走 1-2，2-3 这两条路，即以 1 开头有一条路, 以 2 开头有一条路。如果我们标记 1 为 1，标记 3 为 -1，此时我们将前缀和相加，并用 x[i] 存储, 我们惊奇的发现 x[1]=1, x[2]=1, x[3]=0, 恰好就是城市 i 到城市 i+1 的次数。对每一段出差都这么处理，可以大大简化计算的过程。最后我们判断每一段路走过的次数及买卡和买票的花费，选择花费小的即可。
AC代码如下：

#include <iostream>
#include <queue>
#include <algorithm>
using namespace std;
const int maxn=100010;
typedef long long ll;
ll n,m,x,y,z,p1,p2,now,a[maxn],b[maxn],c[maxn],d[maxn]={0};
ll nsum=0,sum=0;
queue<ll> q;
int main(){
	cin >> n >> m;
	for(int i=0;i<m;i++){
		cin >> now;
		q.push(now);
	}
	for(int i=1;i<n;i++){
		cin >> a[i] >> b[i] >> c[i];
	}
	x=q.front();
	q.pop();
	while(!q.empty()){
		y=q.front();
		q.pop();
		p1=min(x,y);
		p2=max(x,y);
		//cout << p1 << " " << p2 << endl;
		d[p1]++;
		d[p2]--;
		x=y;
	}
	for(int i=1;i<=n;i++){
		d[i]+=d[i-1];
		//cout << d[i] << " ";
	}
	for(int i=1;i<=n-1;i++){
		sum+=min(a[i]*d[i],(b[i]*d[i]+c[i]));
	}
	cout << sum;
	return 0;
}

6. P3397 地毯
解题思路：此题和第五题类似，但是实际上只需要暴力模拟即可。

#include <iostream>
using namespace std;
int mp[1001][1001]={0};
int main(){
	int n,m,x1,x2,y1,y2;
	cin >> n >> m;
	for(int k=0;k<m;k++){
		cin >> x1 >> y1 >> x2 >> y2;
		for(int i=x1;i<=x2;i++){
			for(int j=y1;j<=y2;j++){
				mp[i][j]++;
			}
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(j!=1){
				cout << " ";
			}
			cout << mp[i][j];
		}
		cout << endl;
	}
	return 0;
}