#include #define rep(i,n) for(int i=0;i<(n);i++) using namespace std; using lint=long long; template class mcf_graph{ struct edge{ int to,rev; capa_t capa,flow; cost_t cost; edge(int to,int rev,const capa_t& capa,const cost_t& cost,const capa_t& flow):to(to),rev(rev),capa(capa),cost(cost),flow(flow){} }; vector> G; static constexpr capa_t CAPA_INF=numeric_limits::max(); static constexpr cost_t COST_INF=numeric_limits::max(); public: mcf_graph(){} mcf_graph(int n):G(n){} void add_directed_edge(int u,int v,const capa_t& capa,const cost_t& cost){ G[u].emplace_back(v,G[v].size() ,capa, cost,0); G[v].emplace_back(u,G[u].size()-1, 0,-cost,0); } pair minimum_cost_flow(int s,int t){ int n=G.size(); vector pre(n); vector d(n),pot(n); auto augment=[&]()->pair{ rep(u,n) d[u]=(u==s?0:COST_INF); // Dijkstra bool ok=false; priority_queue> Q; Q.emplace(0,s); while(!Q.empty()){ int u=Q.top().second; cost_t cost=-Q.top().first; Q.pop(); if(cost0) { cost_t cost2=cost+e.cost+pot[u]-pot[e.to]; if(cost2 G(n+1); rep(i,m){ int u,v; lint c1,c2; scanf("%d%d%lld%lld",&u,&v,&c1,&c2); u--; v--; G.add_directed_edge(u,v,1,c1); G.add_directed_edge(u,v,1,c2); G.add_directed_edge(v,u,1,c1); G.add_directed_edge(v,u,1,c2); } G.add_directed_edge(n-1,n,2,0); printf("%lld\n",G.minimum_cost_flow(0,n).second); return 0; }