agc036D Negative Cycle

agc036D Negative Cycle

• 有一个$N$个点的有向图，编号$[1,n]$，首先有$n-1$条边，第$i$条为$(i,i+1)$，边权为0。
• 然后对于所有有序对$(i,j)(i!=j)$，如果$i<j$，边权为$-1$，否则边权为$1$,共加入$n\ast (n-1)$条边，每一个有序对都有一个贡献。
• 要求将后面新加入的$n\ast (n-1)$条边中的若干边删去，使得图上没有负环，求最小的贡献。
• $1\le N\le 500$

Solution

• 想了很久一个负环在以上条件上的形态，根据形态来DP，但是状态实在是不好划分。
• 由于这是一个有顺序的图，我们可以考虑一个差分约束的模型，即从$1$开始走最短路，设$d_x$为到达$x$的最短路长度，那么对于一条边$(i,j)$，有$d_i+w[i,j]\ge d_j$。
• 那么不难发现，负环一定不能满足上面的条件，反之则没有负环。
• 由于$d_i\ge d_{i+1}$，可以发现$c_i=d_i-d_{i+1}$，有$0\le c_i\le 1$。
• 考虑将$n-1$个$c_i$放$0/1$，并将不满足这个序列的边删去，对于$(i,j)(i<j)$，有${\sum }_{k=i}^{j-1}{c}_i\ge 1$，对于$(i,j)(i>j)$，有${\sum }_{k=j}^{i-1}{c}_i\le 1$。
• $n^3$的简单DP，记录最后两个$1$的位置，预处理二维前缀和即可简单转移。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define maxn 505
#define ll long long 
using namespace std;
	
int n,i,j,k,s[maxn][maxn];
ll a[maxn][maxn],b[maxn][maxn],f[maxn][maxn];

ll Suma(int x,int xx,int y,int yy){
	if (x>xx||y>yy) return 0;
	return a[xx][yy]-a[x-1][yy]-a[xx][y-1]+a[x-1][y-1];
}

ll Sumb(int l,int r){
	return b[r][r]-b[r][l-1]-b[l-1][r]+b[l-1][l-1];
}
	
int main(){
	freopen("ceshi.in","r",stdin);
	scanf("%d",&n);
	for(i=1;i<=n;i++){
		for(j=1;j<i;j++) scanf("%d",&a[i][j]);
		for(j=i+1;j<=n;j++) scanf("%d",&a[i][j]);
		for(j=1;j<=n;j++) s[i][j]=a[i][j];
	}
	for(i=1;i<=n;i++) for(j=1;j<=n;j++) 
		b[i][j]=a[i][j]*(i<j)+b[i-1][j]+b[i][j-1]-b[i-1][j-1];
	for(i=1;i<=n;i++) for(j=1;j<=n;j++)
		a[i][j]+=a[i-1][j]+a[i][j-1]-a[i-1][j-1];
	memset(f,127,sizeof(f)),f[0][0]=0;
	for(i=0;i<n;i++) for(j=0;j<=i;j++) if (f[i][j]<1e18){
		for(k=i+1;k<=n;k++)
			f[k][i]=min(f[k][i],f[i][j]+Sumb(i+1,k)+Suma(i+1,k,1,j));
	}
	ll ans=1e18;
	for(i=0;i<=n;i++) ans=min(ans,f[n][i]);
	printf("%lld\n",ans);
}