其他分享
首页 > 其他分享> > [CF321E] Ciel and Gondolas

[CF321E] Ciel and Gondolas

作者:互联网

前言

决策单调性得好好学。

理解不深,瞎掰扯。

题目

洛谷

CF

讲解

经典题。

这道题的决策单调性可以由四边形不等式证明,形如 \(w(a,c)+w(b,d)\le w(a,d)+w(b,c) (a\le b<c\le d).\)

决策单调性可以用分治快速求解,我们考虑分 \(k'\) 块,显然其只和 \(k'-1\) 有关,所以我们一层一层解决问题。

具体的,我们每次二分一个中点,求出它的最优决策点,然后左右两个区间分别递归下去,顺便改一下决策区间,根据决策单调性即可满足复杂度。

时间复杂度 \(O(nk\log_n)\)。

代码

原谅我的工地英语
//12252024832524
#include <bits/stdc++.h>
#define TT template<typename T>
using namespace std;

typedef long long LL;
const int MAXN = 4005;
const int INF = 0x3f3f3f3f;
int n,m;
int pre[MAXN][MAXN],dp[MAXN][MAXN];

LL Read()
{
	LL x = 0,f = 1; char c = getchar();
	while(c > '9' || c < '0'){if(c == '-') f = -1;c = getchar();}
	while(c >= '0' && c <= '9'){x = (x*10) + (c^48);c = getchar();}
	return x * f;
}
TT void Put1(T x)
{
	if(x > 9) Put1(x/10);
	putchar(x%10^48);
}
TT void Put(T x,char c = -1)
{
	if(x < 0) putchar('-'),x = -x;
	Put1(x); if(c >= 0) putchar(c);
}
TT T Max(T x,T y){return x > y ? x : y;}
TT T Min(T x,T y){return x < y ? x : y;}
TT T Abs(T x){return x < 0 ? -x : x;}

int now;
int Sum(int l,int r){return (pre[r][r]-pre[l-1][r]-pre[r][l-1]+pre[l-1][l-1]) >> 1;}
void solve(int l,int r,int ql,int qr)//section to be solved,operation section
{
	if(l > r) return;
	int mid = (l+r) >> 1,ID = ql,MIN = INF;
	for(int i = ql,val;i <= qr && i <= mid;++ i)
		if((val = Sum(i,mid)+dp[i-1][now-1]) < MIN)
		{
			MIN = val;
			ID = i;
		}
	dp[mid][now] = MIN;
	solve(l,mid-1,ql,ID);
	solve(mid+1,r,ID,qr);
}

int main()
{
//	freopen(".in","r",stdin);
//	freopen(".out","w",stdout);
	n = Read(); m = Read();
	for(int i = 1;i <= n;++ i)
		for(int j = 1;j <= n;++ j)
			pre[i][j] = pre[i-1][j]+pre[i][j-1]-pre[i-1][j-1]+Read();
	memset(dp,INF,sizeof(dp));
	dp[0][0] = 0;
	for(now = 1;now <= m;++ now) solve(1,n,1,n);
	Put(dp[n][m],'\n');
	return 0;
}

标签:pre,return,Ciel,int,TT,决策,Gondolas,MAXN,CF321E
来源: https://www.cnblogs.com/PPLPPL/p/15552673.html