题意:
给出一个无向图,类似三角形的样子,然后给出边的权值,问找一条从第一个点到最后一个点的路径,要求每一条边只能走一次,并且权值和最大,点可以重复走。
思路:
首先观察这个图可以发现,所有的点的度数都是偶数。然后由每条边只能走一次知道,这个是和欧拉路相关的,是欧拉道路,不是欧拉回路,因为题目要求是从一个点到另一个点。但是图的所有点的度数都是偶数,那么想办法让图中的第一个点和最后一个点度数变为奇数,其他点的度数都是偶数。这个就比较巧妙,去掉从第一个点到最后一个点的一条无重复点的路径,除了起点和终点度数减1,其它点的度数都减2,目的就达到了。由于题目要求最后走的边的权值和最大,所以去掉的边的权值尽量小,那么从起点到终点求一个最短路即可。
求路径的方法是首先标记最短路上的边,然后从起点或终点开始dfs,走过的每条边标记(注意这是无向图,所以反向路径也要标记),当一个点再无边可走的时候,就把它放入路径中,这样可以保证求出的一定是一个欧拉道路。
代码:
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <vector>
#include <queue>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int N = 305;
const int inf = 0x3f3f3f3f;
struct edge
{
int u,v,cost;
edge(int u,int v,int cost):u(u),v(v),cost(cost){}
edge(){}
};
int a[N][N],b[N][N],c[N][N];
int mp[N][N];
vector<edge> es;
vector<int> G[N*N];
vector<pii> anc;
void adde(int u,int v,int cost)
{
es.push_back(edge(u,v,cost));
es.push_back(edge(v,u,cost));
int sz = es.size();
G[u].push_back(sz-2);
G[v].push_back(sz-1);
}
ll dis[N*N];
int rev[N*N];
bool vis[N*N];
bool used[N*N*8];
pii rid[N*N];
int cnt;
void spfa()
{
for (int i = 0;i <= cnt;i++) dis[i] = 1e18;
memset(vis,0,sizeof(vis));
memset(rev,0,sizeof(rev));
vis[1] = 1;
dis[1] = 0;
queue<int> q;
q.push(1);
while (!q.empty())
{
int u = q.front();
q.pop();
vis[u] = 0;
for (int i = 0;i < G[u].size();i++)
{
edge &e = es[G[u][i]];
int v = e.v;
if (dis[v] > dis[u] + e.cost)
{
dis[v] = dis[u] + e.cost;
rev[v] = G[u][i];
if (!vis[v])
{
vis[v] = 1;
q.push(v);
}
}
}
}
}
void dfs(int u)
{
for (int i = 0;i < G[u].size();i++)
{
int id = G[u][i];
edge &e = es[id];
if (!used[id])
{
used[id] = used[id^1] = 1;
int v = e.v;
dfs(v);
}
}
anc.push_back(rid[u]);
}
void init()
{
cnt = 0;
anc.clear();
es.clear();
for (int i = 0;i < N * N;i++) G[i].clear();
}
int main()
{
int t;
scanf("%d",&t);
while (t--)
{
int n;
scanf("%d",&n);
init();
ll ans = 0;
for (int i = 1;i < n;i++)
{
for (int j = 1;j <= i;j++)
{
scanf("%d",&a[i][j]);
ans += a[i][j];
}
}
for (int i = 1;i < n;i++)
{
for (int j = 1;j <= i;j++)
{
scanf("%d",&b[i][j]);
ans += b[i][j];
}
}
for (int i = 1;i < n;i++)
{
for (int j = 1;j <= i;j++)
{
scanf("%d",&c[i][j]);
ans += c[i][j];
}
}
for (int i = 1;i <= n;i++)
{
for (int j = 1;j <= i;j++)
{
mp[i][j] = ++cnt;
rid[cnt] = pii(i,j);
}
}
for (int i = 1;i < n;i++)
{
for (int j = 1;j <= i;j++)
{
int x = mp[i][j];
int y = mp[i+1][j];
adde(x,y,a[i][j]);
y = mp[i+1][j+1];
adde(x,y,b[i][j]);
x = mp[i+1][j];
y = mp[i+1][j+1];
adde(x,y,c[i][j]);
}
}
spfa();
ans -= dis[cnt];
memset(used,0,sizeof(used));
for (int i = cnt;i != 1;i = es[rev[i]].u)
{
used[rev[i]] = used[rev[i]^1] = 1;
}
dfs(cnt);
printf("%lld\n",ans);
printf("%d\n",(int)anc.size());
for (int i = 0;i < anc.size();i++)
{
printf("%d %d%c",anc[i].first,anc[i].second,i == anc.size() - 1 ? '\n' : ' ');
}
}
return 0;
}