04-前缀和 差分

1、一维前缀和

AcWing 795. 前缀和

输入一个长度为 n的整数序列。

接下来再输入 m个询问，每个询问输入一对 l,r。

对于每个询问，输出原序列中从第 l个数到第 r个数的和。

【输入格式】

第一行包含两个整数 n和 m。

第二行包含 n个整数，表示整数数列。

接下来 m行，每行包含两个整数 l 和 r，表示一个询问的区间范围。

【输出格式】

共 m行，每行输出一个询问的结果。

数据范围

1≤l≤r≤n,

1≤n,m≤100000,

−1000≤数列中元素的值≤1000

输入样例：

5 3

2 1 3 6 4

1 2

1 3

2 4

输出样例：

3

6

10

<代码>

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+6;
int n,m,a[N];
int main()
{
    cin>>n>>m;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        a[i]+=a[i-1];
    }
    for(int i=1;i<=m;i++)
    {
        int l,r;
        scanf("%d%d",&l,&r);
        printf("%d\n",a[r]-a[l-1]);
    }
    return 0;
}

2、二维前缀和

AcWing 796. 子矩阵的和

输入一个 n 行 m 列的整数矩阵，再输入 q 个询问，每个询问包含四个整数 x1,y1,x2,y2，表示一个子矩阵的左上角坐标和右下角坐标。

对于每个询问输出子矩阵中所有数的和。

【输入格式】

第一行包含三个整数 n，m，q。

接下来 n行，每行包含 m个整数，表示整数矩阵。

接下来 q行，每行包含四个整数 x1,y1,x2,y2，表示一组询问。

【输出格式】

共 q行，每行输出一个询问的结果。

数据范围

1≤n,m≤1000,
1≤q≤200000,
1≤x1≤x2≤n,
1≤y1≤y2≤m,
−1000≤矩阵内元素的值≤1000

输入样例：

3 4 3

1 7 2 4

3 6 2 8

2 1 2 3

1 1 2 2

2 1 3 4

1 3 3 4

输出样例：

17

27

21

<代码>

#include<bits/stdc++.h>
using namespace std;
const int N=1e3+5;
int a[N][N],n,m,q;
int main()
{
    cin>>n>>m>>q;
    for(int i=1;i<=n;i++)
      for(int j=1;j<=m;j++)
      {
          scanf("%d",&a[i][j]);
          a[i][j]+=a[i-1][j]+a[i][j-1]-a[i-1][j-1];
      }  
    while(q--)
    {
        int x,y,x2,y2;
        scanf("%d%d%d%d",&x,&y,&x2,&y2);
        printf("%d\n",a[x2][y2]-a[x-1][y2]-a[x2][y-1]+a[x-1][y-1]);
    }
    return 0;
}

3、贪心+前缀和

1224：最大子矩阵

【题目描述】

已知矩阵的大小定义为矩阵中所有元素的和。给定一个矩阵，你的任务是找到最大的非空(大小至少是1×1)子矩阵。

比如，如下4×4的矩阵

的最大子矩阵是

这个子矩阵的大小是15。

【输入】

输入是一个N×N的矩阵。输入的第一行给出N(0<N<=100)。再后面的若干行中，依次(首先从左到右给出第一行的N个整数，再从左到右给出第二行的N个整数……)给出矩阵中的N2个整数，整数之间由空白字符分隔(空格或者空行)。已知矩阵中整数的范围都在[−127,127]。

【输出】

输出最大子矩阵的大小。

【输入样例】

4

 0 -2 -7  0

 9  2 -6  2

-4  1 -4  1

-1  8  0 -2

【输出样例】

15

<代码1>二维前缀和

#include<bits/stdc++.h>
using namespace std;
int n,t[101][101],sum[101][101],ans;
int  main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=n;j++)
	    cin>>t[i][j];
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=n;j++)
	    sum[i][j]=sum[i][j-1]+sum[i-1][j]-sum[i-1][j-1]+t[i][j];
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=n;j++)
	    for(int k=i;k<=n;k++)
	      for(int l=j;l<=n;l++)
	        ans=max(ans,sum[k][l]-sum[i-1][l]-sum[k][j-1]+sum[i-1][j-1]);
	cout<<ans;
	return 0;
}

<代码2>贪心+一维前缀和

#include<bits/stdc++.h>
using namespace std;
const int N=110;
int a[N][N],n,m,ans=-2e9,sum;
int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=n;j++)
	  {
		scanf("%d",&a[i][j]);
		a[i][j]+=a[i][j-1];//行前缀和
	  }
    for(int l=1;l<=n;l++)
		for(int r=l;r<=n;r++)
		{
			sum=0;
			for(int i=1;i<=n;i++)//行
			{
                int tmp=a[i][r]-a[i][l-1];
                sum+=tmp;
				ans=max(ans,sum);
				if(sum<0)sum=0;
			}
		}

	cout<<ans<<endl;
	return 0;
}

4、一维差分

AcWing 797. 差分

输入一个长度为 n的整数序列。

接下来输入 m个操作，每个操作包含三个整数 l,r,c，表示将序列中 [l,r] 之间的每个数加上 c。

请你输出进行完所有操作后的序列。

【输入格式】

第一行包含两个整数 n和 m。

第二行包含 n个整数，表示整数序列。

接下来 m行，每行包含三个整数 l，r，c，表示一个操作。

【输出格式】

共一行，包含 n个整数，表示最终序列。

数据范围

1≤n,m≤100000,
1≤l≤r≤n,
−1000≤c≤1000,
−1000≤整数序列中元素的值≤1000

输入样例：

6 3

1 2 2 1 2 1

1 3 1

3 5 1

1 6 1

输出样例：

3 4 5 3 4 2

<代码>

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+5;
int a[N],b[N],n,m,l,r,c;
int main()
{
    cin>>n>>m;
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&l,&r,&c);
        b[l]+=c;
        b[r+1]-=c;
    }
    for(int i=1;i<=n;i++)
    {
        b[i]+=b[i-1];
        printf("%d ",a[i]+b[i]);
    }
    return 0;
}

5、二维差分

AcWing 798. 差分矩阵

输入一个 n 行 m 列的整数矩阵，再输入 q 个操作，每个操作包含五个整数 x1,y1,x2,y2,c，其中 (x1,y1) 和 (x2,y2)表示一个子矩阵的左上角坐标和右下角坐标。

每个操作都要将选中的子矩阵中的每个元素的值加上 c。

请你将进行完所有操作后的矩阵输出。

【输入格式】

第一行包含整数 n,m,q。

接下来 n行，每行包含 m个整数，表示整数矩阵。

接下来 q行，每行包含 5 个整数 x1,y1,x2,y2,c，表示一个操作。

【输出格式】

共 n行，每行 m个整数，表示所有操作进行完毕后的最终矩阵。

数据范围

1≤n,m≤1000,
1≤q≤100000,
1≤x1≤x2≤n,
1≤y1≤y2≤m,
−1000≤c≤1000,
−1000≤矩阵内元素的值≤1000

输入样例：

3 4 3

1 2 2 1

3 2 2 1

1 1 1 1

1 1 2 2 1

1 3 2 3 2

3 1 3 4 1

输出样例：

2 3 4 1

4 3 4 1

2 2 2 2

<代码>

#include<bits/stdc++.h>
using namespace std;
const int N=1e3+5;
int a[N][N],b[N][N],n,m,q;
int main()
{
    cin>>n>>m>>q;
    for(int i=1;i<=n;i++)
      for(int j=1;j<=m;j++)
        scanf("%d",&a[i][j]);
    while(q--)
    {
        int x,y,x2,y2,c;
        scanf("%d%d%d%d%d",&x,&y,&x2,&y2,&c);
        b[x][y]+=c,b[x2+1][y2+1]+=c;
        b[x][y2+1]-=c,b[x2+1][y]-=c;
    }
    for(int i=1;i<=n;i++)
    {     
        for(int j=1;j<=m;j++)
          {
              b[i][j]+=b[i-1][j]+b[i][j-1]-b[i-1][j-1];
              printf("%d ",a[i][j]+b[i][j]);
          }
        printf("\n");
    }
    return 0;
}