【学习笔记】[AGC036D] Negative Cycle

21 篇文章 0 订阅

首先,没有负环等价于差分约束有解。这谁想得到啊qwq

那么记变量 { x i } \{x_i\} {xi},相当于满足若干不等式。假设我们得到了一组合法的 { x i } \{x_i\} {xi},那么显然 x 1 ≥ x 2 ≥ . . . ≥ x n x_1\ge x_2\ge...\ge x_n x1x2...xn,记 q i = x i − x i + 1 q_i=x_i-x_{i+1} qi=xixi+1,如果 i < j i<j i<j有负边那么 x i − x j ≥ 1 x_i-x_j\ge 1 xixj1,等价于 [ i , j ) [i,j) [i,j)不全为 0 0 0,如果 i < j i<j i<j有正边那么 x i − x j ≤ 1 x_i-x_j\le 1 xixj1,显然我们可以看出 q i ∈ { 0 , 1 } q_i\in \{0,1\} qi{0,1}。那么正边的情况也等价于 [ i , j ) [i,j) [i,j)最多只有一个 1 1 1

那么我们可以对 { q i } \{q_i\} {qi}进行 d p dp dp,去计算删边的花费。

复杂度 O ( n 3 ) O(n^3) O(n3)

#include<bits/stdc++.h>
#define ll long long
#define inf 0x3f3f3f3f3f3f3f3f
#define pb push_back
#define pii pair<int,int>
#define fi first
#define se second
using namespace std;
int n;
ll a[505][505],b[505][505],f[505][505],g[505][505],res(inf);
void chmin(ll &x,ll y){
	x=min(x,y);
}
signed main(){
	cin>>n;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(i!=j){
				if(i<j)cin>>a[j-1][i];
				else cin>>b[i-1][j];
			}
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=i;j++){
			a[i][j]+=a[i][j-1],b[i][j]+=b[i][j-1];
		}
	}
	memset(f,0x3f,sizeof f),f[0][0]=0;
	for(int i=1;i<n;i++){
		memset(g,0x3f,sizeof g);
		g[0][0]=f[0][0]+a[i][i];
		g[0][i]=f[0][0];
		for(int j=1;j<i;j++){
			chmin(g[0][j],f[0][j]+a[i][i]-a[i][j]);
			chmin(g[j][i],f[0][j]+b[i][j]);
		}
		for(int j=1;j<i;j++){
			for(int k=j+1;k<i;k++){
				if(f[j][k]==inf)continue;
				chmin(g[j][k],f[j][k]+(a[i][i]-a[i][k])+b[i][j]);
				chmin(g[k][i],f[j][k]+b[i][k]);
			}
		}memcpy(f,g,sizeof g);
	}for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			res=min(res,f[i][j]);
		}
	}cout<<res;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值