[sgu]185. Two shortest 网络流

我要疯啦~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

网上看见的解法都是spfa+最大流，我偏去搞最小费用最大流……

写完了觉得C++压代码就是爽……于是交，MLE，才看见空间4096K……

不死心瞎改，修改过程中反复出现了WA CE RE TLE MLE……

Edge结构体里全改short发现空间不超了，又开始RE……

然后发现有一个成员变量必须int,算了下3200+K 应该到3900K左右可以搞……交，MLE

发现因为有结构体好像也在吃空间，全改成纯数组，AC了……一共交28遍，历时9小时

现在想想貌似用临接矩阵就可以水过去了？总之题思路真的算不上难T_T只能说我太弱了，颓一个月的后果……

/*
  两遍spfa显然过不了，反例想下也很好举……
  很容易能想到网络流建模，这方面还是比较水的……
  两种寻路方法
  一种是先spfa再枚举找出来dist[a]+len(a,b)==dist[b]的边搞出来个“最短路图”
  两条路径就是两道流，最大流很容易搞……
  另一种是最小费用最大流，建图也没什么可说的……不过网上不通行说明这么做效率有问题？
  改进下应该其实算是两遍spfa的正确姿势吧，思路和朴素最大流相似，既然会走错路就加上回边
  实现起来和mcmf还是有区别的，只进行两次1-n网络增广判断下两次增量是不是都一样
  输出的话两种方法一样，直接dfs输出就行了
*/
#include <cstdio>
#include <cstring>
#include <cstdlib>
const int N=401,M=320002,INF=999999999;

struct Caf{int flow,cost;};
int g[N],q[N],prep[N],pree[N],dist[N],flow[N],tot=1,ps,pt,n,m,next[M];
short p[M],c[M],k[M];
bool v[N];
int min(int a,int b) {return a<b?a:b;}
void add(int a,int b,int cc,int d) {
  tot++; p[tot]=b; c[tot]=cc; k[tot]=d; next[tot]=g[a]; g[a]=tot;
  tot++; p[tot]=a; c[tot]=0; k[tot]=-d; next[tot]=g[b]; g[b]=tot;
}
Caf spfa() {
  int head=0,tail=1;
  memset(v,false,sizeof(v)); prep[ps]=pree[ps]=0;
  for (int i=1; i<=N; i++) dist[i]=INF;
  dist[ps]=0; q[1]=ps; flow[ps]=INF;
  while (head!=tail) {
    v[q[head=head%N+1]]=false;
    for (int tmp=g[q[head]]; tmp; tmp=next[tmp])
      if ((dist[p[tmp]]>dist[q[head]]+(int)k[tmp])&&c[tmp]) {
         dist[p[tmp]]=dist[q[head]]+(int)k[tmp];
         prep[p[tmp]]=q[head]; pree[p[tmp]]=tmp;
         flow[p[tmp]]=min(flow[q[head]],c[tmp]);
         if (!v[p[tmp]]) v[q[tail=tail%N+1]=p[tmp]]=true;
      }
  }
  if (dist[pt]==INF) return (Caf){0,0};
  for (int tmp=pt; tmp!=ps; tmp=prep[tmp])
    c[pree[tmp]]-=flow[pt], c[pree[tmp]^1]+=flow[pt];
  return (Caf){flow[pt],flow[pt]*dist[pt]};
}
Caf mcmf() {return spfa();}
void print(int pt) {
  if (pt==n) {printf("%d\n",pt); return;}
  for (int tmp=g[pt]; tmp; tmp=next[tmp])
    if ((c[tmp]==0)&&(tmp%2==0)) {
      printf("%d ",pt); c[tmp]=99;
      print(p[tmp]); return;
    }
}
int main() {
  memset(g,0,sizeof(g));
  //===============
  scanf("%d%d",&n,&m);
  while (m--) {
    int a,b,c;
    scanf("%d%d%d",&a,&b,&c);
    add(a,b,1,c); add(b,a,1,c);
  }
  ps=1; pt=n; Caf t1=spfa(),t2=spfa();
  if ((t1.flow==t2.flow==1)&&(t1.cost==t2.cost)) print(1),print(1);
  else puts("No solution");
  //system("pause");
  return 0;
}