传送门:http://acm.hdu.edu.cn/showproblem.php?pid=5306
三种操作:
0 x y t: 使这个区间的值变成min(t,a[i])。
1 x y: 区间最大值
2 x y: 区间和
1和2都是常见操作,就是这个0操作很难搞定。参考了吉如一论文。感觉非常的强。
区间的max可以看做是区间标记,我们可以维护3个值,区间最大值,区间次大值,区间最大值的个数。然后我们每次做0操作的时候,就会有3种情况。
把t作为要更新的值。
1:t>=区间最大值,那么就是无效操作。
2:区间次大值 < t <区间最大值,因为我们已经求得了最大值的个数,所以我们可以直接更新这段的sum。
3:其他情况。继续搜下去,直到搜到前两种情况为止。
然后我们可以把区间的最大值看做是这个区间的标记,因为每次更新的时候只会更新到(区间次大值 < t < 区间最大值)的情况,而且因为没有加法操作,所以max从上到下一定是非递增的,那么我们在更新的时候就可以从每次pushdown的时候判断一下子节点和当前节点的最大值的情况,然后更新一下子节点的sum和max值。复杂度证明并不会,但是论文里说是O(mlogn)的。
2000+ms AC
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <string>
#include <iostream>
#include <vector>
#include <map>
#include <set>
#include <stack>
#include <queue>
#include <ctime>
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
#define pb push_back
#define mp make_pair
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define calm (l+r)>>1
const int INF = 2139062143;
const int maxn=1000010;
int n,m;
struct Seg{
ll sum[maxn<<2];
int mx[maxn<<2],mxnum[maxn<<2],sc[maxn<<2];
inline void pushup(int rt){
sum[rt]=sum[rt<<1]+sum[rt<<1|1];
mx[rt]=max(mx[rt<<1],mx[rt<<1|1]);
mxnum[rt]=0;
sc[rt]=max(sc[rt<<1],sc[rt<<1|1]);
if(mx[rt]==mx[rt<<1])mxnum[rt]+=mxnum[rt<<1];
else sc[rt]=max(sc[rt],mx[rt<<1]);
if(mx[rt]==mx[rt<<1|1])mxnum[rt]+=mxnum[rt<<1|1];
else sc[rt]=max(sc[rt],mx[rt<<1|1]);
}
inline void pushdown(int rt){
if(mx[rt]<mx[rt<<1]){
sum[rt<<1]-=(ll)mxnum[rt<<1]*(mx[rt<<1]-mx[rt]);
mx[rt<<1]=mx[rt];
}
if(mx[rt]<mx[rt<<1|1]){
sum[rt<<1|1]-=(ll)mxnum[rt<<1|1]*(mx[rt<<1|1]-mx[rt]);
mx[rt<<1|1]=mx[rt];
}
}
void build(int l,int r,int rt){
if(l==r){
scanf("%I64d",&sum[rt]);
mx[rt]=sum[rt];
mxnum[rt]=1;
sc[rt]=-1;
return;
}
int m=calm;
build(lson);build(rson);
pushup(rt);
}
void update(int L,int R,int v,int l,int r,int rt){
if(v>=mx[rt])return;
if(L<=l&&r<=R&&sc[rt]<v){
sum[rt]-=(ll)mxnum[rt]*(mx[rt]-v);
mx[rt]=v;
return;
}
pushdown(rt);
int m=calm;
if(L<=m)update(L,R,v,lson);
if(R>m)update(L,R,v,rson);
pushup(rt);
}
int getMax(int L,int R,int l,int r,int rt){
if(L<=l&&r<=R){
return mx[rt];
}
pushdown(rt);
int m=calm,ans=0;
if(L<=m)ans=getMax(L,R,lson);
if(R>m)ans=max(ans,getMax(L,R,rson));
return ans;
}
ll query(int L,int R,int l,int r,int rt){
if(L<=l&&r<=R){
return sum[rt];
}
pushdown(rt);
int m=calm;
ll ans=0;
if(L<=m)ans+=query(L,R,lson);
if(R>m)ans+=query(L,R,rson);
return ans;
}
}tree;
int main(){
//freopen("D://input.txt","r",stdin);
int T;scanf("%d",&T);
while(T--){
scanf("%d%d",&n,&m);
tree.build(1,n,1);
while(m--){
int op,l,r;
scanf("%d%d%d",&op,&l,&r);
if(op==0){
int x;scanf("%d",&x);
tree.update(l,r,x,1,n,1);
}
else if(op==1){
printf("%d\n",tree.getMax(l,r,1,n,1));
}
else{
printf("%I64d\n",tree.query(l,r,1,n,1));
}
}
}
return 0;
}