洛谷—— P2196 挖地雷

https://www.luogu.org/problem/show?pid=2196

题目背景

NOIp1996提高组第三题

题目描述

在一个地图上有N个地窖（N<=20），每个地窖中埋有一定数量的地雷。同时，给出地窖之间的连接路径。当地窖及其连接的数据给出之后，某人可以从任一处开始挖地雷，然后可以沿着指出的连接往下挖（仅能选择一条路径），当无连接时挖地雷工作结束。设计一个挖地雷的方案，使某人能挖到最多的地雷。

输入输出格式

输入格式：

输入文件mine.in有若干行。

第1行只有一个数字，表示地窖的个数N。

第2行有N个数，分别表示每个地窖中的地雷个数。

第3行至第N+1行表示地窖之间的连接情况：

第3行有n-1个数（0或1），表示第一个地窖至第2个、第3个、…、第n个地窖有否路径连接。如第3行为1 1 0 0 0 … 0，则表示第1个地窖至第2个地窖有路径，至第3个地窖有路径，至第4个地窖、第5个、…、第n个地窖没有路径。

第4行有n-2个数，表示第二个地窖至第3个、第4个、…、第n个地窖有否路径连接。

… …

第n+1行有1个数，表示第n-1个地窖至第n个地窖有否路径连接。（为0表示没有路径，为1表示有路径）。

输出格式：

输出文件wdl.out有两行数据。

第一行表示挖得最多地雷时的挖地雷的顺序，各地窖序号间以一个空格分隔，不得有多余的空格。

第二行只有一个数，表示能挖到的最多地雷数。

输入输出样例

输入样例#1：

5
10 8 4 7 6
1 1 1 0
0 0 0
1 1
1

输出样例#1：

1 3 4 5
27

建立超级源点s，从s跑一边最短路，记录每个在最短路上的点的前驱、处理一下输出便可、

 1 #include <cstdio>
 2 #include <queue>
 3
 4 #define max(a,b) (a>b?a:b)
 5 inline void read(int &x)
 6 {
 7     x=0; register char ch=getchar();
 8     for(; ch>‘9‘||ch<‘0‘; ) ch=getchar();
 9     for(; ch>=‘0‘&&ch<=‘9‘; ch=getchar()) x=x*10+ch-‘0‘;
10 }
11 const int INF(0x3f3f3f3f);
12 const int N(26);
13 int val[N],head[N],sumedge;
14 struct Edge {
15     int v,next,w;
16     Edge(int v=0,int next=0,int w=0):v(v),next(next),w(w){}
17 }edge[N*N];
18 inline void ins(int u,int v,int w)
19 {
20     edge[++sumedge]=Edge(v,head[u],w); head[u]=sumedge;
21 }
22
23 bool inq[N];
24 int dis[N],pre[N],ans[N];
25 std::queue<int>que;
26 inline void SPFA(int s)
27 {
28     que.push(s);
29     for(int u,v; !que.empty(); )
30     {
31         u=que.front(); que.pop(); inq[u]=0;
32         for(int i=head[u]; i; i=edge[i].next)
33         {
34             v=edge[i].v;
35             if(dis[v]<dis[u]+edge[i].w)
36             {
37                 pre[v]=u;
38                 dis[v]=dis[u]+edge[i].w;
39                 if(!inq[v]) inq[v]=1,que.push(v);
40             }
41         }
42     }
43 }
44
45 int Presist()
46 {
47     int n; read(n);
48     for(int i=1; i<=n; ++i)
49         read(val[i]),ins(0,i,val[i]);
50     for(int i=1,x; i<n; ++i)
51       for(int j=i+1; j<=n; ++j)
52         { read(x); if(x) ins(i,j,val[j]); }
53     SPFA(0); int maxdis=0,t,cnt=0;
54     for(int i=1; i<=n; ++i)
55         if(maxdis<dis[i]) maxdis=dis[t=i];
56     for(; t; t=pre[t]) ans[++cnt]=t;
57     for(int i=cnt; i; --i) printf("%d ",ans[i]);
58     printf("\n%d",maxdis);
59     return 0;
60 }
61
62 int Aptal=Presist();
63 int main(int argc,char*argv[]){;}