【高手训练】【最短路】负环 记录

设$dis[k][i][j]$表示从点$i$走到点$j$走了$2^k$条边的最短路径长度。

$asd[i][j]$表示二分到当前答案时点$i$到点$j$的最短路。

$dis[k][i][j]$用$FloydO(n^3\log n)$预处理

对于任意一个正整数$x$，它显然可以被表示成$2$的整数次幂的和。

所以可以从大到小枚举$k$，每次用$dis[k][i][j]$去更新$asd[i][j]$。表示从点$i$到点$j$的路径上常识一次走$2^k$条边

$qwe[i][j]=min(qwe[i][j],min(dis[z][i][k]+asd[k][j],dis[z][k][j]+asd[i][k]))$

此处需要另外开一个数组$qwe$来储存更新的结果，如果直接在原$asd$数组上更新会导致$asd[i][j]$又被刚更新完的$asd[i][z],asd[z][j]$更新，这样的话一次就不止尝试走$2^k$步了。

更新后枚举原图的某一条边，如果存在$w_i+qwe[u_i][v_i]<0$，则原图存在负环。

当前答案肯定大于最终答案，否则要把$asd$数组赋成更新后的$qwe$，把答案加上$2^k$

但是此做法是有点问题的，从点$i$到点$j$走$k$条边的最短路不一定大于等于走$k+1$条边的最短路，即不满足单调性

对此有一个巧妙的方法就是给每个点加一条长度为$0$的自环即可，这样对答案没有影响

复杂度$O(n^3\log n)$

Code:

#include<bits/stdc++.h>
using namespace std;
#define f1(a,b,c) for(int c=a;c<=b;c++)
#define f2(a,b,c) for(int c=a;c>=b;c--)
#define f3(a,b,c) for(int c=a;c;c=b)
#define so1(a,n) sort(a+1,a+n+1,mycmp);
#define so2(a,n) sort(a+1,a+n+1);
#define ll long long
#define itn int 
#define ubt int 
#define mp make_pair
#define pii pair<int,int>
#define pll pair<ll,ll>
const int twx=300+100;
const int inf=0x3f3f3f3f;
ll read()
{
    ll sum=0;
    ll flag=1;
    char c=getchar();
    while(c<'0'||c>'9')
    {
        if(c=='-')
        {
            flag=-1;
        }
        c=getchar();
    }
    while(c>='0'&&c<='9')
    {
        sum=((sum*10)+c-'0');
        c=getchar();
    }
    return sum*flag;
}
int n;
int m;
int dis[15][twx][twx];
int asd[twx][twx];
int qwe[twx][twx];
int ans=0;
void init()
{
    n=read();
    m=read();
    memset(dis,0x3f,sizeof dis);
    memset(asd,0x3f,sizeof asd);
    f1(1,m,i)
    {
        int x=read();
        int y=read();
        int z=read();
        dis[0][x][y]=z;
    }
}
void work()
{
    f1(1,n,i)
    {
        dis[0][i][i]=0;//自环
    }
    f1(0,8,z)
    {
        f1(1,n,k)
        {
            f1(1,n,i)
            {
                f1(1,n,j)
                {
                    dis[z+1][i][j]=min(dis[z+1][i][j],dis[z][i][k]+dis[z][k][j]);
                }
            }
        }  
    }
    f1(1,n,i)
    {
        asd[i][i]=0;
    }
    f2(8,0,z)
    {
        f1(1,n,i)
        {
            f1(1,n,j)
            {
                qwe[i][j]=asd[i][j];
            }
        }
        f1(1,n,k)
        {
            f1(1,n,i)
            {
                f1(1,n,j)
                {
                    qwe[i][j]=min(qwe[i][j],min(dis[z][i][k]+asd[k][j],dis[z][k][j]+asd[i][k]));
                }
            }
        }   
        int flag=0;
        f1(1,n,i)
        {
            if(qwe[i][i]<0)
            {
                flag=1;
                break ;
            }
        }
        if(!flag)
        {
            ans|=1<<z;
            f1(1,n,i)
            {
                f1(1,n,j)
                {
                    asd[i][j]=qwe[i][j];
                }
            }
        }
    }
}
void print()
{
    printf("%d\n",ans>=n?0:ans+1);   
}
int main()
{
    //freopen(".in","r",stdin);
    //freopen(".out","w",stdout);
    init();
    work();
    print();
	return 0;
}