http://poj.org/problem?id=2076
#include<iostream>
#include"algorithm"
using namespace std;
int dist[101][101],n,m;
void init()
{
int i,a,b,d,j,k;
cin>>n>>m;
for(i=0;i<n;i++)
for(j=0;j<n;j++)
dist[i][j]=9999;
for(i=0;i<m;i++)
{
cin>>a>>b>>d;
a--;b--;
if(d<dist[a][b])
{
dist[a][b]=dist[b][a]=d;
}
}
for(i=0;i<n;i++)
dist[i][i]=0;
for(k=0;k<n;k++)
for(i=0;i<n;i++)
for(j=0;j<n;j++)
if(dist[i][k]+dist[k][j]<dist[i][j])
dist[i][j]=dist[i][k]+dist[k][j];
return ;
}
int answer[101];
int c1,c2;
int cmp(int a,int b)
{
return dist[a][c2]-dist[a][c1]>dist[b][c2]-dist[b][c1]||
(dist[a][c2]-dist[a][c1]==dist[b][c2]-dist[b][c1]&&a<b);
}
void doit()
{
int index[101];long best=99999999,now;
int bestk;
int i,j;
for(c2=0;c2<n;c2++)
for(c1=c2+1;c1<n;c1++)
{
bestk=-1;
for(i=0,j=0;i<n;i++)
if(i!=c1&&i!=c2)
index[j++]=i;
std::sort(&index[0],&index[j],cmp);
now=0;
for(i=0;i<n;i++)
if(i!=c1&&i!=c2)now+=(n-1)*dist[i][c2];
now+=(n-1)*dist[c1][c2];
if(now<best)
{
bestk=0;
best=now;
}
for(i=1;i<n-1;i++)
{
now-=i*(n-i)*dist[c1][c2];
now+=(i+1)*(n-i-1)*dist[c1][c2];
now+=(dist[index[i-1]][c1]-dist[index[i-1]][c2])*(n-1);
if(now<best)
{
bestk=i;
best=now;
}
}
if(bestk>=0)
{
for(i=0;i<bestk;i++)
answer[index[i]]=c1+1;
for(;i<n-2;i++)
answer[index[i]]=c2+1;
answer[c1]=0;
answer[c2]=0;
}
}
return;
}
int main()
{
int t,i;
cin>>t;
while(t--)
{
init();
doit();
for(i=0;i<n;i++)
{
if(i)cout<<' ';
cout<<answer[i];
}
cout<<endl;
}
return 0;
}