[SDOI 2010]スターウォーズ

昨日は2つのネットワークストリームを連続して作りました・・・これが1つ目です
これは道図論題が肯定的で,図は全部君にあげた.
問題はモデリング方法です
問題は各ポイントへのアクセスがちょうど1回であることを要求している(私は最初この条件を見ていなかった......)
総時間が最も短いことを要求して、問題をいくつかの経典の図論の問題に転化して例えば最も短いことを試みます
残念ながらだめなので、自然とネットの流れを思い浮かべます(グループの中には「すべてネットの流れができる」という戯言があります.例えばA+B......)
さらに分析すると単純なネットワークストリームではだめであり,費用フローが必要であることが分かった.
各ポイントにアクセスするのはちょうど1回で、パスオーバーライドに似ていますが、実は少し似ています.
各星を2つの点に分解し、uとu'
私たちは各テーマに与えられたエッジ(u,v)に対して,ネットワークストリームにエッジ(u,v')を1つ追加し,トラフィックは1であり,費用は時間である.
その後、最初に任意の点に移動することができ、stからv'に1つのエッジを接続し、流量は1であり、費用は位置決め時間である.
各v′からedに1つのエッジを接続し、容量は1、費用は0であり、各ポイントの入度が1であり、1回のみアクセスできることを示す.
stから各uに1つのエッジを接続し、容量は1であり、費用は0であり、uが(u,v')を通ってv'に達するようにする.
では、この図の最小費用フローが答えです.
今までEK版しか書けなかった費用の流れ・・・スピード超遅い・・バカ×一つは...

//Lib
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>

#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<queue>
using namespace std;
//Macro
#define	rep(i,a,b)	for(int i=a,tt=b;i<=tt;++i)
#define	drep(i,a,b)	for(int i=a,tt=b;i>=tt;--i)
#define	erep(i,e,x)	for(int i=x;i;i=e[i].next)
#define	irep(i,x)	for(__typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define	read()	(strtol(ipos,&ipos,10))
#define	sqr(x)	((x)*(x))
#define	pb	push_back
#define	PS	system("pause");
typedef	long long	ll;
typedef	pair<int,int>	pii;
const int oo=~0U>>1;
const double inf=1e100;
const double eps=1e-6;
string name="", in=".in", out=".out";
//Var
struct E
{
	int next,node,cap,v;
}e[40008];
queue<int> q;
int n,m,st,ed,ans,tot=1;
int h[1608],dis[1608],pre[1608],epre[1608];
bool vis[1608];
void add(int a,int b,int c,int d)
{
	e[++tot].next=h[a];e[tot].node=b;e[tot].cap=c;e[tot].v=d;h[a]=tot;
	e[++tot].next=h[b];e[tot].node=a;e[tot].cap=0;e[tot].v=-d;h[b]=tot;
}
void Init()
{
	scanf("%d%d",&n,&m);
	st=n*2+1;ed=n*2+2;
	int a,b,c;
	rep(i,1,n)
	{
		scanf("%d",&c);
//		add(st,i,1,0);
		add(st,n+i,1,c);
		add(n+i,ed,1,0);
	}
	rep(i,1,m)
	{
		scanf("%d%d%d",&a,&b,&c);
		if(a>b)swap(a,b);
		add(a,b+n,1,c);
	}
}
bool SPFA()
{
	memset(dis,55,sizeof dis);
	q.push(st);dis[st]=0;
	int u,v;bool flag=false;
	while(!q.empty())
	{
		u=q.front();q.pop();vis[u]=false;
		erep(i,e,h[u])
			if(e[i].cap&&dis[v=e[i].node]>dis[u]+e[i].v)
			{
				dis[v]=dis[u]+e[i].v;
				pre[v]=u;epre[v]=i;
				if(!vis[v])q.push(v),vis[v]=true;
				if(v==ed)flag=true;
			}
	}
	return flag;
}
int Flow()
{
	int tmp,ret=0,flow=oo;
	tmp=ed;
	while(tmp!=st)
	{
		flow=min(flow,e[epre[tmp]].cap);
		ret+=e[epre[tmp]].v;
		tmp=pre[tmp];
	}
	tmp=ed;
	while(tmp!=st)
	{
		e[epre[tmp]].cap-=flow;
		e[epre[tmp]^1].cap+=flow;
		tmp=pre[tmp];
	}
	return ret*flow;
}
void Work()
{
	while(SPFA())
		ans+=Flow();
	printf("%d
",ans);
}
int main()
{
//	freopen((name+in).c_str(),"r",stdin);
//	freopen((name+out).c_str(),"w",stdout);
	Init();
	Work();
	return 0;
}