标签:https 等价 ic卡 signed log har 方式 复杂 判断
考虑一个数组\(A\),其项数为\(n\)项。有\(m\)次询问,每次询问给定两个参数\(l\)和\(r\),要求求出\(A[l]+A[l+1]+...+A[r]\)。
怎么做呢?
那么前缀和究竟是什么呢?
在读入\(A\)数组的时候,我们预处理出一个前缀和数组\(S\),满足\(S[i]=A[1]+A[2]+...+A[i]\)。如何快速计算\(S\)?显然有\(S[i]=S[i-1]+A[i]\),\(S\)数组就可以这样线性递推出来了。
那么对于每次查询,就有\(A[l]+A[l+1]+...+A[r]=S[r]-S[l-1]\),实现了\(O(n)\)预处理,\(O(1)\)回答单次查询。
那么什么信息可以用这样的“前缀”方式维护呢?只要满足可加可减性也就是满足结合律的信息,都可以这样快速计算维护。例如:快速求区间和,快速求区间积,快速求区间异或和等等。
那么差分又是什么呢?差分是前缀和的逆运算。
在读入\(A\)数组的时候,我们预处理出一个差分数组\(T\),满足\(T[i]=A[i]-A[i-1]\)。差分数组就是原数组相邻两项的差。
考虑对\(n\)项的数组\(A\)给出\(m\)次询问,每次询问给定三个参数\(l\)和\(r\)和\(d\),要求对\(A[l],A[l+1],...,A[r]\)均执行\(+d\)操作。最后求出\(A\)数组的每一项的值。
利用差分运算,对\(A[l],A[l+1],...,A[r]\)均执行\(+d\)操作,就等价于对差分数据\(T\)进行\(T[l]+=d,T[r+1]-=d\)的操作。
那么直接维护差分数组\(T\),对差分数组做前缀和即可还原出原数组\(A\)。
和很多数据结构一样,前缀和与差分也可以推广到二维矩阵下。
二维前缀和:\(S[x][y]=S[x][y-1]+S[x-1][y]-S[x-1][y-1]+A[x][y]\)。就等价于一个容斥原理。
二维差分:对 \((x1,y1)\)~\((x2,y2)\)的 \(A\)数组\(+d\),等价于 \(T[x1][y1]+=d,T[x1][y2+1]?=d,T[x2+1][y1]?=d,T[x2+1][y2+1]+=d\)。知道\(T\)还原\(A\)做二维前缀和即可。
也可以将二维前缀和/二维差分,暴力拆成\(n\)个一维前缀和/一维差分进行运算。
P1115 最大子段和
最大子段和是一个经典问题,这里给出一种利用前缀和的线性做法。先预处理出一个前缀和数组\(S\)。
一个子段和\(A[l]+A[l+1]+...+A[r]\)就等价于\(S[r]-S[l-1]\)。要使得这个子段和最大,在维持\(S[r]\)为定值的同时,还要使得\(S[l-1]\)尽可能小。那么从左往右枚举右端点\(r\),并维护左侧的所有\(S\)中的最小值\(mins\)数组,两者相减即为此时的最大子段和,尝试更新答案。
#include<bits/stdc++.h>//P1115 最大子段和
using namespace std;
#define re register
#define ll long long
#define il inline
#define dou double
#define un unsigned
il int read()
{
char c=getchar();int x=0,f=1;
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
return x*f;
}
#define INF 114514114
#define clr(x) memset(x,0,sizeof(x))
#define N 200000+10
int n,ans=-INF;
int a[N],s[N],mins[N];
int main()
{
n=read();
for(re int i=1;i<=n;i++)mins[i]=INF;
for(re int i=1;i<=n;i++)
{
a[i]=read();
s[i]=s[i-1]+a[i];
mins[i]=min(mins[i-1],s[i]);
}
for(re int r=1;r<=n;r++)
ans=max(ans,s[r]-mins[r-1]);
cout<<ans<<endl;
return 0;
}
二维前缀和与二维差分的板子题。可以暴力拆成\(n\)个一维前缀和/一维差分维护。
#include<bits/stdc++.h>//P3397 地毯
using namespace std;
#define re register
#define ll long long
#define il inline
#define dou double
#define un unsigned
il int read()
{
char c=getchar();int x=0,f=1;
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
return x*f;
}
#define INF 114514114
#define clr(x) memset(x,0,sizeof(x))
#define N 1000+10
int n,m,a,b,c,d;
int tag[N][N],sum[N][N];
int main()
{
n=read();m=read();
for(re int i=1;i<=m;i++)
{
a=read();b=read();c=read();d=read();
for(re int j=a;j<=c;j++)
{
tag[j][b]++;
tag[j][d+1]--;
}
}
for(re int i=1;i<=n;i++)
{
for(re int j=1;j<=n;j++)
{
sum[i][j]=sum[i][j-1]+tag[i][j];
printf("%d ",sum[i][j]);
}
printf("\n");
}
return 0;
}
用差分维护每一段道路经过的次数,对于每一段道路,记经过次数为\(s\),那么它的最小花费就是\(min(as,bs+c)\)。即购买纸质票和购买IC卡的费用取最小值。注意每一段道路起点和终点的编号先后顺序可能需要调整。
#include<bits/stdc++.h>//P3406 海底高铁
using namespace std;
#define re register
#define ll long long
#define il inline
#define dou double
#define un unsigned
il int read()
{
char c=getchar();int x=0,f=1;
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
return x*f;
}
#define INF 114514114
#define clr(x) memset(x,0,sizeof(x))
#define N 100000+10
#define M 100000+10
ll n,m,ans;
ll p[M],a[N],b[N],c[N],tag[N],s[N];
int main()
{
n=read();m=read();
for(re int i=1;i<=m;i++)p[i]=read();
for(re int i=1;i<=n-1;i++)
{
a[i]=read();b[i]=read();c[i]=read();
}
for(re int i=1;i<=m-1;i++)
{
if(p[i]<p[i+1])
{
tag[p[i]]++;
tag[p[i+1]]--;
}
else
{
tag[p[i+1]]++;
tag[p[i]]--;
}
}
for(re int i=1;i<=n;i++)s[i]=s[i-1]+tag[i];
for(re int i=1;i<=n;i++)
ans=ans+min(a[i]*s[i],b[i]*s[i]+c[i]);
cout<<ans<<endl;
return 0;
}
二分答案\(x\),只需判断\(x\)为答案是否合法。
如何\(check\)?用差分维护借教室情况,前缀和还原之后比较,判断是否合法。
#include<bits/stdc++.h>
using namespace std;
#define re register
#define ll long long
#define il inline
#define dou double
#define un unsigned
il int read()
{
char c=getchar();int x=0,f=1;
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
return x*f;
}
#define N 1000000+10
#define M 1000000+10
#define clr(x) memset(x,0,sizeof(x))
int n,m;
int a[N],d[M],s[M],t[M],tag[N],sum[N];
il bool check(int x)
{
clr(tag);clr(sum);
for(re int i=1;i<=x;i++)
{
tag[s[i]]+=d[i];
tag[t[i]+1]-=d[i];
}
for(re int i=1;i<=n;i++)
{
sum[i]=sum[i-1]+tag[i];
if(sum[i]>a[i])
return false;
}
return true;
}
int main()
{
n=read();m=read();
for(re int i=1;i<=n;i++)
a[i]=read();
for(re int i=1;i<=m;i++)
{
d[i]=read();s[i]=read();t[i]=read();
}
int l=1,r=m;
if(check(m))
{
cout<<0<<endl;
return 0;
}
while(l<r)
{
int mid=(l+r)>>1;
if(check(mid))
l=mid+1;
else
r=mid;
}
cout<<-1<<endl<<l<<endl;
return 0;
}
标签:https 等价 ic卡 signed log har 方式 复杂 判断
原文地址:https://www.cnblogs.com/Hakurei-Reimu/p/11621447.html