【USACO 2013 February Gold】旅行线路

Description

贝西经营着一家旅行社，一天贝西带着几队游客沿着亚马逊河旅行，河的两边分布着一些景点，每个景点都对应着一个观赏值。景点间由一些穿过河流的道路相连(位于河流同一侧的景点间没有直接道路相连)，贝西想要设计游客的旅行线路，使得该线路经过的景点的总的观赏值最大。 

但是，贝西可能同时带着几个旅行团，贝西希望它给安排的旅行线路不能相交。 
两条线路 (a <-> x) 和 (b <-> y) 出现如下情况之一就算相交: 
if(a < b and y < x) or (b < a and x < y) or (a = b and x = y). 

帮助贝西找出最佳旅行线路，贝西可以从任意景点开始，在任意景点结束。

Input

第一行，三个空格间隔的整数N (1 <= N <= 40,000), M (1 <=M <= 40,000), and R (0 <= R <= 100,000) 。N表示河岸左边的景点数，M表示河岸右边的景点数，R表示道路的条数 
接下来N行，每行一个整数，表示河岸左边对应景点的观赏值。 
接下来M行，每行一个整数，表示河岸右边对应景点的观赏值。 
接下来R行，每行两个空格间隔的整数X,Y，表示左岸的X号景点与右岸的Y号景点间有道路直接相连。 

(0 <= 观赏值<= 40,000),

Output

一个整数，表示所求最大观赏值。

Sample Input

3 2 4
1
1
5
2
2
1 1
2 1
3 1
2 2

Sample Output

8

【分析】

DP

仔细分析题目后我们发现，我们走的路线一定是在L方向严格递增，R方向严格递增。我们自然想到以点作为状态，但这样时空复杂度都是O(N^2)级别的，不能承受。于是我们考虑以边作为事件来解决。

我们先对边排序，左端点由大到小。这时候被这条边连接的两个点便可以状态转移。我们设DPA[i]表示左边以i结尾的最大收益，DPB[i]表示右边以i结尾的最大收益。

得到状态转移方程：

存在边<i,j>时

    DPA[i]=max{A[i]+DPB[j]} 

    DPB[j]=max{B[j]+DPA[i]} 

转移的时候拿两个临时变量存一下就行了

【代码】

/**********************
    ID:Ciocio
	LANG:C++
	DATE:2013-12-24
	TASK:route
**********************/
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>
#include <functional>

using namespace std;

#define MAXN 40005
#define MAXE 100005

typedef long long LL;
typedef pair<int,int> pii;
int N,M,R;
LL A[MAXN],B[MAXN],DPA[MAXN],DPB[MAXN];
pii E[MAXE];

template <class T>
void _read(T &x)
{
	char tt=getchar();
	while(tt<'0'||'9'<tt) tt=getchar();
	for(x=0;'0'<=tt&&tt<='9';x=(x<<3)+(x<<1)+tt-'0',tt=getchar());
}

void _init()
{
	_read(N);_read(M);_read(R);
	for(int i=1;i<=N;i++) _read(A[i]);
	for(int i=1;i<=M;i++) _read(B[i]);
	for(int i=1;i<=R;i++) _read(E[i].first),_read(E[i].second);
	sort(E+1,E+R+1);
}

void _solve()
{
	memcpy(DPA,A,sizeof(A));
	memcpy(DPB,B,sizeof(B));
	for(int i=R;i;i--)
	{
		int u=E[i].first,v=E[i].second;
		LL va=A[u]+DPB[v];
		LL vb=B[v]+DPA[u];
		DPA[u]=max(DPA[u],va);
		DPB[v]=max(DPB[v],vb);
	}
	LL ans=max(*max_element(DPA+1,DPA+N+1),*max_element(DPB+1,DPB+M+1));
	cout<<ans<<endl;
}

int main()
{
	_init();
	_solve();
	return 0;
}