SPFA算法,简单回溯记录路径。
翻译来自scy
Part Acquisition部分收购
题意:
奶牛们被送往太空去执行一个任务:为它们的农场收购一台挤奶机器。它们正通过
一个含有N颗行星的星系。(1<=N<=50 000),每颗行星都是一个交易站。
奶牛们决定了使用 K种(1<=K<=1000)种商品在这些星球中进行交易。这些星球没
有发展货币,所以星球交易都是在货易系统下工作的:所有的交易都是用一种货物交换
另外一种货物。假设互相交易的两种货物是不同的。
奶牛们从地球出发,带着一罐高品质的干草(货物1号),它们想最终用罐干草换回
一台挤奶机器(货物K号)。你们经常喝奶牛们的牛奶,奶牛们请你帮忙:确定如何与
该星系所有星球交易,最后得到得到货物K号,奶牛们要你构造最好的交易序列。如果无
论如何都不可能最后换得货物K号,那么输出“-1”。
Input 输入
* 第一行:N和K
* 第 2..N+1行: 第 i+1 行 包含两个整数 ai 和 bi, 代表星球i用货物ai去换取货物bi.
Output 输出
* 第一行: 一个整数T,即为用货物1号换得货物k号的最少交换序列。(或者 -1。如果无法
换得货物K号 ).
* 第2..T+1行: 列出交换序列中每个货物的编号.
Sample Input
6 5
1 3
3 2
2 3
3 1
2 5
5 4
Sample Output
4
1
3
2
5
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<queue>
#include<vector>
#include<iostream>
using namespace std;
const int N=51000,K=1100,INF=(int)1e6;
int n,k,len;
bool vis[N];
int first[N],d[N],pre[N];
vector<int> ans;
struct node{
int x,y,d,next;
}a[N];
void ins(int x,int y)
{
len++;
a[len].x=x;a[len].y=y;
a[len].next=first[x];first[x]=len;
}
void print(int x)
{
ans.push_back(x);
if(pre[x]) print(pre[x]);
}
int main()
{
int i,j;
scanf("%d%d",&n,&k);
len=0;
memset(first,0,sizeof(first));
for(i=1;i<=n;i++)
{
int x,y;
scanf("%d%d",&x,&y);
ins(x,y);
}
d[1]=0;
for(i=2;i<=n;i++) d[i]=INF;
memset(vis,0,sizeof(vis));
memset(pre,0,sizeof(pre));
queue<int> q;
q.push(1);
while(!q.empty())
{
int x=q.front();q.pop();
vis[x]=1;
for(i=first[x];i;i=a[i].next)
{
int y=a[i].y;
if(d[x]+a[i].d<d[y])
{
d[y]=d[x]+a[i].d;
pre[y]=x;
if(!vis[y])
{
vis[y]=1;
q.push(y);
}
}
}
vis[x]=0;
}
if(d[k]>=INF) printf("-1\n");
else
{
print(k);
printf("%d\n",ans.size());
for(i=ans.size()-1;i>=0;i--) printf("%d\n",ans[i]);
}
}