【bzoj 4392】 Counting Haybales 【Usaco2015 dec】

Description

Farmer John is trying to hire contractors to help rearrange his farm, but so far all of them have quit when they saw the complicated sequence of instructions FJ wanted them to follow. Left to complete the project by himself, he realizes that indeed, he has made the project perhaps more complicated than necessary. Please help him follow his instructions to complete the farm upgrade.

FJ's farm consists of NN fields in a row, conveniently numbered 1…N1…N. In each field there can be any number of haybales. Farmer John's instructions contain three types of entries:

1) Given a contiguous interval of fields, add a new haybale to each field.

2) Given a contiguous interval of fields, determine the minimum number of haybales in a field within that interval.

3) Given a contiguous interval of fields, count the total number of haybales inside that interval.

农夫约翰打算重修他的农场。他有 N 块土地,连续排列成一行,标号为 1…N。在每块土地上有任意数量的草堆。他可以发出三种指令:
1) 对一个连续区间的土地,每块土地增加相同数量的草堆。
2) 对一个连续区间的土地,输出其中最少的草堆数量。
3) 对一个连续区间的土地,输出草堆数量总数。

第一行两个正整数,N (1≤N≤200,000) 和 Q (1≤Q≤100,000)。
下一行是N个非负整数,最大100,000,表示每块土地上有多少个草堆。
以下Q行,每行一单个大写字母开头(M,P或S),空格后跟随两个正整数 A 和 B (1≤A≤B≤N), 或者三个正整数 A, B, 和 C (1≤A≤B≤N; 1≤C≤100,000)。当且仅当第一个字母是 P 时,是三个正整数。
当该字母是M,输出区间A…B的最小草堆数。
当该字母是P,在区间A…B,每块土地增加C堆草。
当该字母是M,输出区间A…B的草堆数之和。

每行一个数字,用于响应'M' 或 'S' 命令。

 

Input

The first line contains two positive integers, N (1≤N≤200,000) and Q (1≤Q≤100,000).

The next line contains N nonnegative integers, each at most 100,000, indicating how many haybales are initially in each field.

Each of the next Q lines contains a single uppercase letter, either M, P or S, followed by either two positive integers AA and BB (1≤A≤B≤N), or three positive integers AA, BB, and CC (1≤A≤B≤N; 1≤C≤100,000). There will be three positive integers if and only if the uppercase letter is P.

If the letter is M, print the minimum number of haybales in the interval of fields from A…B.

If the letter is P, put C new haybales in each field in the interval of fields from A…B.

If the letter is S, print the total number of haybales found within interval of fields from A…B.

Output

A line in the output should appear in response to every 'M' or 'S' entry in FJ's instructions.

Sample Input

4 5
3 1 2 4
M 3 4
S 1 3
P 2 3 1
M 3 4
S 1 3

Sample Output

2
6
3
8

这道题是一个线段树模板题,下面是程序:

#include<stdio.h>
#include<iostream>
#define ll long long
using namespace std;
const int N=200005;
struct treenode{
	int lc,rc;
	ll sum,c,mn;
}t[N<<2];
int k=1;
ll min(ll a,ll b){
	return a<b?a:b;
}
inline void read(ll &s){
	s=0;
	char c=getchar();
	while(c<'0'||c>'9'){
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		s=(s<<3)+(s<<1)+c-'0';
		c=getchar();
	}
}
inline void build(int rt,int l,int r){
	if(l==r){
		read(t[rt].mn);
		t[rt].sum=t[rt].mn;
		t[rt].c=t[rt].lc=t[rt].rc=0;
		return;
	}
	int m=l+r>>1;
	build(t[rt].lc=++k,l,m);
	build(t[rt].rc=++k,m+1,r);
	t[rt].sum=t[t[rt].lc].sum+t[t[rt].rc].sum;
	t[rt].mn=min(t[t[rt].lc].mn,t[t[rt].rc].mn);
}
inline void add(int rt,int l,int r,int a,int b,ll &w){
	t[rt].sum+=(b-a+1)*w;
	if(l==a&&r==b){
		t[rt].mn+=w;
		t[rt].c+=w;
		return;
	}
	int m=l+r>>1;
	if(t[rt].c){
		t[t[rt].lc].c+=t[rt].c;
		t[t[rt].rc].c+=t[rt].c;
		t[t[rt].lc].sum+=(m-l+1)*t[rt].c;
		t[t[rt].rc].sum+=(r-m)*t[rt].c;
		t[t[rt].lc].mn+=t[rt].c;
		t[t[rt].rc].mn+=t[rt].c;
		t[rt].c=0;
	}
	if(b<=m){
		add(t[rt].lc,l,m,a,b,w);
	}
	else{
		if(a>m){
			add(t[rt].rc,m+1,r,a,b,w);
		}
		else{
			add(t[rt].lc,l,m,a,m,w);
			add(t[rt].rc,m+1,r,m+1,b,w);
		}
	}
	t[rt].sum=t[t[rt].lc].sum+t[t[rt].rc].sum;
	t[rt].mn=min(t[t[rt].lc].mn,t[t[rt].rc].mn);
}
inline ll askmn(int rt,int l,int r,int a,int b){
	if(l==a&&r==b){
		return t[rt].mn;
	}
	int m=l+r>>1;
	if(t[rt].c){
		t[t[rt].lc].c+=t[rt].c;
		t[t[rt].rc].c+=t[rt].c;
		t[t[rt].lc].sum+=(m-l+1)*t[rt].c;
		t[t[rt].rc].sum+=(r-m)*t[rt].c;
		t[t[rt].lc].mn+=t[rt].c;
		t[t[rt].rc].mn+=t[rt].c;
		t[rt].c=0;
	}
	if(b<=m){
		return askmn(t[rt].lc,l,m,a,b);
	}
	if(a>m){
		return askmn(t[rt].rc,m+1,r,a,b);
	}
	return min(askmn(t[rt].lc,l,m,a,m),askmn(t[rt].rc,m+1,r,m+1,b));
}
inline ll asksum(int rt,int l,int r,int a,int b){
	if(l==a&&r==b){
		return t[rt].sum;
	}
	int m=l+r>>1;
	if(t[rt].c){
		t[t[rt].lc].c+=t[rt].c;
		t[t[rt].rc].c+=t[rt].c;
		t[t[rt].lc].sum+=(m-l+1)*t[rt].c;
		t[t[rt].rc].sum+=(r-m)*t[rt].c;
		t[t[rt].lc].mn+=t[rt].c;
		t[t[rt].rc].mn+=t[rt].c;
		t[rt].c=0;
	}
	if(b<=m){
		return asksum(t[rt].lc,l,m,a,b);
	}
	if(a>m){
		return asksum(t[rt].rc,m+1,r,a,b);
	}
	return asksum(t[rt].lc,l,m,a,m)+asksum(t[rt].rc,m+1,r,m+1,b);
}
inline char ch(){
	char c=getchar();
	while(c!='M'&&c!='S'&&c!='P'){
		c=getchar();
	}
	return c;
}
inline void out(ll n){
	if(n>9){
		out(n/10);
	}
	putchar(n%10+'0');
}
int main(){
	int n,m,l,r;
	ll w;
	char c;
	scanf("%d%d",&n,&m);
	build(1,1,n);
	while(m--){
		c=ch();
		scanf("%d%d",&l,&r);
		switch(c){
			case 'S':
				out(asksum(1,1,n,l,r));
				putchar('\n');
				break;
			case 'M':
				out(askmn(1,1,n,l,r));
				putchar('\n');
				break;
			case 'P':
				read(w);
				add(1,1,n,l,r,w);
				break;
		}
	}
	return 0;
}

 

题目描述 牛牛和她的朋友们正在玩一个有趣的游戏,他们需要构建一个有 $n$ 个节点的无向图,每个节点都有一个唯一的编号并且编号从 $1$ 到 $n$。他们需要从节点 $1$ 到节点 $n$ 找到一条最短路径,其中路径长度是经过的边权的和。为了让游戏更有趣,他们决定在图上添加一些额外的边,这些边的权值都是 $x$。他们想知道,如果他们添加的边数尽可能少,最短路径的长度最多会增加多少。 输入格式 第一行包含两个正整数 $n$ 和 $m$,表示节点数和边数。 接下来 $m$ 行,每行包含三个整数 $u_i,v_i,w_i$,表示一条无向边 $(u_i,v_i)$,权值为 $w_i$。 输出格式 输出一个整数,表示最短路径的长度最多会增加多少。 数据范围 $2 \leq n \leq 200$ $1 \leq m \leq n(n-1)/2$ $1 \leq w_i \leq 10^6$ 输入样例 #1: 4 4 1 2 2 2 3 3 3 4 4 4 1 5 输出样例 #1: 5 输入样例 #2: 4 3 1 2 1 2 3 2 3 4 3 输出样例 #2: 2 算法 (BFS+最短路) $O(n^3)$ 我们把问题转化一下,假设原图中没有添加边,所求的就是点 $1$ 到点 $n$ 的最短路,并且我们已经求出了这个最短路的长度 $dis$。 接下来我们从小到大枚举边权 $x$,每次将 $x$ 加入图中,然后再次求解点 $1$ 到点 $n$ 的最短路 $dis'$,那么增加的最短路长度就是 $dis'-dis$。 我们发现,每次加入一个边都需要重新求解最短路。如果我们使用 Dijkstra 算法的话,每次加入一条边需要 $O(m\log m)$ 的时间复杂度,总的时间复杂度就是 $O(m^2\log m)$,无法通过本题。因此我们需要使用更优秀的算法。 观察到 $n$ 的范围比较小,我们可以考虑使用 BFS 求解最短路。如果边权均为 $1$,那么 BFS 可以在 $O(m)$ 的时间复杂度内求解最短路。那么如果我们只是加入了一条边的话,我们可以将边权为 $x$ 的边看做 $x$ 条边的组合,每次加入该边时,我们就在原始图上添加 $x$ 条边,边权均为 $1$。这样,我们就可以使用一次 BFS 求解最短路了。 但是,我们不得不考虑加入多条边的情况。如果我们还是将边权为 $x$ 的边看做 $x$ 条边的组合,那么我们就需要加入 $x$ 条边,而不是一条边。这样,我们就不能使用 BFS 了。 但是,我们可以使用 Floyd 算法。事实上,我们每次加入边时,只有边权等于 $x$ 的边会发生变化。因此,如果我们枚举边权 $x$ 时,每次只需要将边权等于 $x$ 的边加入图中,然后使用 Floyd 算法重新计算最短路即可。由于 Floyd 算法的时间复杂度为 $O(n^3)$,因此总的时间复杂度为 $O(n^4)$。 时间复杂度 $O(n^4)$ 空间复杂度 $O(n^2)$ C++ 代码 注意点:Floyd算法计算任意两点之间的最短路径,只需要在之前的路径基础上加入新的边构成的新路径进行更新即可。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值