【学习笔记】[AGC036D] Negative Cycle

难以评价。。。

首先,没有负环等价于差分约束有解。这谁想得到啊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; }

__EOF__

本文作者仰望星空的蚂蚁
本文链接https://www.cnblogs.com/cqbzly/p/17530041.html
关于博主:评论和私信会在第一时间回复。或者直接私信我。
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   仰望星空的蚂蚁  阅读(7)  评论(0编辑  收藏  举报  
相关博文:
阅读排行:
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· DeepSeek 开源周回顾「GitHub 热点速览」
点击右上角即可分享
微信分享提示