「P3385」模板负环（spfa

Posted 2021-01-17 qwerta

题目描述

暴力枚举/SPFA/Bellman-ford/奇怪的贪心/超神搜索

输入输出格式

输入格式：

 

第一行一个正整数T表示数据组数，对于每组数据：

第一行两个正整数N M，表示图有N个顶点，M条边

接下来M行，每行三个整数a b w，表示a->b有一条权值为w的边（若w<0则为单向，否则双向）

 

输出格式：

 

共T行。对于每组数据，存在负环则输出一行"YE5"(不含引号)，否则输出一行"N0"(不含引号)。

 

输入输出样例

输入样例#1： 复制

2
3 4
1 2 2
1 3 4
2 3 1
3 1 -3
3 3
1 2 3
2 3 4
3 1 -8

输出样例#1： 复制

N0
YE5

说明

nleq 2000n≤2000mleq 3000m≤3000-10000leq wleq 10000−10000≤w≤10000Tleq 10T≤10建议复制输出格式中的字符串。 本题数据感谢@negiizhao的精心构造，请不要使用玄学算法 本题数据有更新

题解

枚举每个点，如果该点未被遍历过，则以其为起点跑spfa。（防止图不连通

跑的过程中记一下每个点被经历的次数，如果被经历（学术说法松弛？）了超过n次，说明这里出现了负环。

算是spfa的一个小技巧？

————————————

然后是年度悲情大戏之——我的STLqueue不可能这么慢↓

我的妈我这么爱用的STLqueue这么慢？！！

老姐你把我的时间直接拖长了一倍不止啊？！！（气哭

还以为自己哪里写的太丑了......

 1 /*
 2 qwerta 
 3 P3385 【模板】负环 Accepted 
 4 100
 5 代码 C++，1.24KB
 6 提交时间 2018-11-03 15:37:36
 7 耗时/内存 1759ms, 15192KB
 8 */
 9 #include<iostream>
10 #include<cstring>
11 #include<cstdio>
12 #include<queue>
13 using namespace std;
14 struct emm{
15     int e,f,l;
16 }a[6003];
17 int h[2003];
18 int tot;
19 void con(int x,int y,int l)
20 {
21     //cout<<"con "<<x<<" "<<y<<" "<<l<<endl;
22     a[++tot].f=h[x];
23     h[x]=tot;
24     a[tot].e=y;
25     a[tot].l=l;
26     return;
27 }
28 const int INF=1e9+7;
29 int d[2003];
30 int c[2003];
31 bool sf[2003];
32 bool vis[2003];
33 int q[4000003];//时间减半的手写队列
34 int main()
35 {
36     //freopen("a.in","r",stdin);
37     int t;
38     scanf("%d",&t);
39     while(t--)
40     {
41         tot=0;
42         memset(h,0,sizeof(h));
43         int n,m;
44         scanf("%d%d",&n,&m);
45         for(int i=1;i<=m;++i)
46         {
47             int x,y,l;
48             scanf("%d%d%d",&x,&y,&l);
49             if(l<0)
50             con(x,y,l);
51             else
52             con(x,y,l),con(y,x,l);
53         }
54         int flag=0;
55         memset(d,127,sizeof(d));
56         memset(c,0,sizeof(c));
57         memset(vis,0,sizeof(vis));
58         for(int s=n;s>=1&&!flag;--s)
59         if(!vis[s])
60         {
61             memset(sf,0,sizeof(sf));
62             d[s]=0;sf[s]=1;q[1]=s;
63             int l=0,r=1;
64             while(l<r&&!flag)
65             {
66                 int x=q[++l];
67                 vis[x]=1;
68                 for(int i=h[x];i;i=a[i].f)
69                 if(d[a[i].e]>d[x]+a[i].l)
70                 {
71                     d[a[i].e]=d[x]+a[i].l;
72                     c[a[i].e]++;
73                     if(c[a[i].e]>=n){flag++;break;}
74                     if(!sf[a[i].e])
75                     {
76                         sf[a[i].e]=1;
77                         q[++r]=a[i].e;
78                     }
79                 }
80                 sf[x]=0;
81             }
82             d[s]=INF;
83         }
84         if(flag)cout<<"YE5"<<endl;
85         else cout<<"N0"<<endl;
86     }
87     return 0;
88 }