hdu-1874 SPFA

noip day1 少得100分,看错题了。。。

/*
 * hdu-1874 SPFA
 * mike-w
 * 2011-11-12
 */
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define MAXN 500
#define MAXM 5000
#define QSIZE 5000
#define INF 999999

typedef struct _lnode
{
	int node,weight,next;
}lnode;

int hhead[MAXN],ID;
int que[MAXN],inq[MAXN],head,tail,len;
lnode list[MAXM];
int M,N,S,T;
int d[MAXN];

int push(int e)
{
	que[tail++]=e;
	len++;
	if(tail==QSIZE)
		tail=0;
	return 0;
}

int pop(int *e)
{
	*e=que[head++];
	len--;
	if(head==QSIZE)
		head=0;
	return 0;
}

int insert(int start,int end,int weight)
{
	int pos,nxt;
	if(hhead[start]==-1)
	{
		pos=hhead[start]=ID++;
		list[pos].node=end;
		list[pos].weight=weight;
		list[pos].next=-1;
	}
	else
	{
		pos=hhead[start];
		nxt=list[pos].next;
		pos=list[pos].next=ID++;
		list[pos].node=end;
		list[pos].weight=weight;
		list[pos].next=nxt;
	}
	return 0;
}

int SPFA(void)
{
	int i,x,pos;
	head=tail=len=0;
	for(i=0;i<N;i++)
	{
		inq[i]=head=tail=len=0;
		d[i]=INF;
	}

	push(S);
	inq[S]=1;
	d[S]=0;
	while(len>0)
	{
		pop(&x);
		inq[x]=0;
		for(pos=hhead[x];pos!=-1;pos=list[pos].next)
		{
			if(d[x]+list[pos].weight<d[list[pos].node])
			{
				d[list[pos].node]=d[x]+list[pos].weight;
				if(!inq[list[pos].node])
					push(list[pos].node),inq[list[pos].node]=1;
			}
		}
	}
	if(d[T]==INF)
		return -1;
	else
		return d[T];
}

int main(void)
{
	int i,t1,t2,t3;
#ifndef ONLINE_JUDGE
	freopen("in","r",stdin);
#endif
	while(scanf("%d%d",&N,&M)!=EOF)
	{
		for(i=0;i<=N;i++)
			hhead[i]=-1;
		ID=0;
		for(i=1;i<=M;i++)
		{
			scanf("%d%d%d",&t1,&t2,&t3);
			insert(t1,t2,t3);
			insert(t2,t1,t3);
		}
		scanf("%d%d",&S,&T);
		printf("%d\n",SPFA());
	}
	return 0;
}



 

以下是hdu4310的Java解法: ```java import java.util.*; import java.io.*; public class Main { static int MAXN = 100010; static int MAXM = 200010; static int INF = 0x3f3f3f3f; static int n, m, s, t, cnt; static int[] head = new int[MAXN]; static int[] dis = new int[MAXN]; static boolean[] vis = new boolean[MAXN]; static int[] pre = new int[MAXN]; static int[] cur = new int[MAXN]; static class Edge { int to, next, cap, flow, cost; public Edge(int to, int next, int cap, int flow, int cost) { this.to = to; this.next = next; this.cap = cap; this.flow = flow; this.cost = cost; } } static Edge[] edge = new Edge[MAXM]; static void addEdge(int u, int v, int cap, int flow, int cost) { edge[cnt] = new Edge(v, head[u], cap, flow, cost); head[u] = cnt++; edge[cnt] = new Edge(u, head[v], 0, 0, -cost); head[v] = cnt++; } static boolean spfa() { Arrays.fill(dis, INF); Arrays.fill(vis, false); Queue<Integer> q = new LinkedList<>(); q.offer(s); dis[s] = 0; vis[s] = true; while (!q.isEmpty()) { int u = q.poll(); vis[u] = false; for (int i = head[u]; i != -1; i = edge[i].next) { int v = edge[i].to; if (edge[i].cap > edge[i].flow && dis[v] > dis[u] + edge[i].cost) { dis[v] = dis[u] + edge[i].cost; pre[v] = u; cur[v] = i; if (!vis[v]) { vis[v] = true; q.offer(v); } } } } return dis[t] != INF; } static int[] MCMF(int s, int t) { int flow = 0, cost = 0; while (spfa()) { int f = INF; for (int u = t; u != s; u = pre[u]) { f = Math.min(f, edge[cur[u]].cap - edge[cur[u]].flow); } for (int u = t; u != s; u = pre[u]) { edge[cur[u]].flow += f; edge[cur[u] ^ 1].flow -= f; cost += edge[cur[u]].cost * f; } flow += f; } return new int[]{flow, cost}; } public static void main(String[] args) { Scanner in = new Scanner(new BufferedReader(new InputStreamReader(System.in))); int T = in.nextInt(); for (int cas = 1; cas <= T; cas++) { n = in.nextInt(); m = in.nextInt(); s = 0; t = n + m + 1; cnt = 0; Arrays.fill(head, -1); for (int i = 1; i <= n; i++) { int c = in.nextInt(); addEdge(s, i, c, 0, 0); } for (int i = 1; i <= m; i++) { int c = in.nextInt(); addEdge(i + n, t, c, 0, 0); } for (int i = 1; i <= n; i++) { for (int j = 1; j <= m; j++) { int c = in.nextInt(); addEdge(i, j + n, INF, 0, c); } } int[] ans = MCMF(s, t); System.out.printf("Case #%d: %d\n", cas, ans[1]); } } } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值