P3089 [USACO13NOV] Pogo-Cow S

合集 - 题解(30)

1.[ABC363G] Dynamic Scheduling 与 P4511 [CTSC2015] 日程管理2024-07-262.P9058 [Ynoi2004] rpmtdq 与 P9678 [ICPC2022 Jinan R] Tree Distance2024-07-293.P1081 [NOIP2012 提高组] 开车旅行2024-07-304.P1973 [NOI2011] NOI 嘉年华2024-07-305.P1398 [NOI2013] 书法家2024-07-316.P3957 [NOIP2017 普及组] 跳房子2024-07-317.P5017 [NOIP2018 普及组] 摆渡车2024-08-018.P5665 [CSP-S2019] 划分2024-08-019.P4544 [USACO10NOV] Buying Feed G2024-08-02

10.P3089 [USACO13NOV] Pogo-Cow S2024-08-02

11.P3572 [POI2014] PTA-Little Bird2024-08-0312.P6764 [APIO2020] 粉刷墙壁2024-08-0313.P2831 [NOIP2016 提高组] 愤怒的小鸟2024-08-0414.P3959 [NOIP2017 提高组] 宝藏2024-08-0515.P2150 [NOI2015] 寿司晚宴2024-08-0516.P8037 [COCI2015-2016#7] Prokletnik2024-08-1217.CF1943C Tree Compass2024-08-1318.P1084 [NOIP2012 提高组] 疫情控制2024-08-1419.P6805 [CEOI2020] 春季大扫除2024-08-1620.P4689 [Ynoi2016] 这是我自己的发明 与 P5268 [SNOI2017] 一个简单的询问02024-08-1821.CF1264D1/2 Beautiful Bracket Sequence (easy/hard version)2024-08-2122.CF924D Contact ATC2024-08-2223.P9640 [SNCPC2019] Digit Mode2024-08-2324.P2757 [国家集训队] 等差子序列 与 CF452F Permutation2024-08-2625.P4423 [BJWC2011] 最小三角形 与 SP7209 CLOSEST - Closest Triplet2024-08-2926.P2825 [HEOI2016/TJOI2016] 游戏 与 P10945 Place the Robots2024-08-2927.P8304 [CoE R4 D] 01 串2024-08-3028.P3320 [SDOI2015] 寻宝游戏 与 P10930 异象石 与 CF176E Archaeology2024-08-3029.P9108 [PA2020] Malowanie płotu2024-09-0330.CF773D Perishable Roads2024-09-03

收起

思路：

首先要先将目标点按照坐标从小到大排序。

考虑动态规划算法，定义 $f_{i,j}$ 表示从 $j\to i$ 的最大得分：

$$f_{i,j}=p_i+\underset{k=i+1}{\overset{n}{max}}[x_k-x_i\ge x_i-x_j]f_{k,i}$$

朴素实现是 $O(N^3)$ 的但是数据过水，可以冲过去。

考虑优化，注意一下 $[x_k-x_i\ge x_i-x_j]$ 的条件，即 $[x_k\ge 2x_i-x_j]$，则当 $j$ 单增时，$2x_i-x_j$ 单减，那么我们可以走指针维护 $k$，维护中统计 $f_{k,i}$ 的最大值即可。

时间复杂度为 $O(N^2)$。

注意还可以往左走，状态方程类似，就不多说了。

完整代码：

#include<bits/stdc++.h>
#define Add(x,y) (x+y>=mod)?(x+y-mod):(x+y)
#define lowbit(x) x&(-x)
#define pi pair<ll,ll>
#define pii pair<ll,pair<ll,ll>>
#define iip pair<pair<ll,ll>,ll>
#define ppii pair<pair<ll,ll>,pair<ll,ll>>
#define fi first
#define se second
#define full(l,r,x) for(auto it=l;it!=r;it++) (*it)=x
#define Full(a) memset(a,0,sizeof(a))
#define open(s1,s2) freopen(s1,"r",stdin),freopen(s2,"w",stdout);
using namespace std;
typedef double db;
typedef unsigned long long ull;
typedef long long ll;
bool Begin;
const ll N=1010; 
inline ll read(){
    ll x=0,f=1;
    char c=getchar();
    while(c<'0'||c>'9'){
        if(c=='-')
          f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9'){
        x=(x<<1)+(x<<3)+(c^48);
        c=getchar();
    }
    return x*f;
}
inline void write(ll x){
	if(x<0){
		putchar('-');
		x=-x;
	}
	if(x>9)
	  write(x/10);
	putchar(x%10+'0');
}
struct Node{
	ll x,w;
	bool operator<(const Node&rhs)const{
		return x<rhs.x;
	}
}A[N];
ll n,ans,s;
ll a[N],b[N];
ll f[N][N],g[N][N];
bool End;
int main(){
	n=read();
	for(int i=1;i<=n;i++)
	  A[i]={read(),read()};
	sort(A+1,A+n+1);
	for(int i=1;i<=n;i++){
		a[i]=A[i].x;
		b[i]=A[i].w;
	}
	for(int i=1;i<=n;i++)
	  f[n][i]=b[n];
	for(int i=n-1;i>=1;i--){
		ll Max=0;
		for(int k=n,j=1;j<=i;j++){
			while(a[k]>=2ll*a[i]-a[j]&&k>=i+1){
				Max=max(Max,f[k][i]);
				k--;
			}
			f[i][j]=Max+b[i];
		}
	}
	g[1][1]=b[1];
	for(int i=1;i<=n;i++){
		ll Max=0;
		for(int j=1,k=n;k>=i+1;k--){
			while(a[j]<=2ll*a[i]-a[k]&&j<i){
				Max=max(Max,g[i][j]);
				j++;
			}
			g[k][i]=max({g[k][i],b[i]+b[k],Max+b[k]});	
		}
	}
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=i;j++)
	    ans=max({ans,f[i][j],g[i][j]});
	write(ans);
	cerr<<'\n'<<abs(&Begin-&End)/1048576<<"MB";
	return 0;
}