TOJ 1883 ZOJ 1298 POJ 1135 Domino Effect / Dijkstra

把第一个推倒 求所有多米诺倒下的最短时间 不是每个点上有一个骨牌 是路上都有 每个点只是关键的骨牌

只有一种最优解 先求出1到各点的最短路 求最短路中的最大值 这是一个关键点的情况

在枚举每一条边 u v w 如果(d[u] + d[v] + w) / 2 > max 就更新一下 因为是要求最短路的最大值  这是2个关键点的情况

Dijkstra

#include <cstring>
#include <cstdio>
#include <vector>
#include <queue>
using namespace std;
const int maxn = 510;
const double INF = 999999999;

struct HeapNode
{
	double d;
	int u;
	bool operator < (const HeapNode& rhs) const
	{
		return d > rhs.d;
	}
};

struct Edge
{
	int from, to;
	double dist;
};
vector <Edge> edges;
vector <int> G[maxn];
bool done[maxn];
double d[maxn];
int n, m, t;

void AddEdge(int from, int to, double dist)
{
	edges.push_back((Edge){from, to, dist});
	edges.push_back((Edge){to, from, dist});
	m = edges.size();
	G[from].push_back(m-2);
	G[to].push_back(m-1);
}

void Dijkstra()
{
	priority_queue <HeapNode> Q;
	for(int i = 1; i <= n; i++)
		d[i] = INF;
	d[1] = 0;
	memset(done, 0, sizeof(done));
	Q.push((HeapNode){0, 1});
	while(!Q.empty())
	{
		HeapNode x = Q.top();Q.pop();
		int u = x.u;
		if(done[u])
			continue;
		done[u] = true;
		//printf("%d\n", u);
		for(int i = 0; i < G[u].size(); i++)
		{
			Edge& e = edges[G[u][i]];
			if(d[e.to] > d[u] + e.dist)
			{
				d[e.to] = d[u] + e.dist;
				Q.push((HeapNode){d[e.to], e.to});
			} 
		}
	}
}
int main()
{
	int cas = 1;
	while(scanf("%d %d", &n, &t) && (n+t))
	{
		edges.clear();
		for(int i = 1;i <= n; i++)
			G[i].clear();
		for(int i = 1; i <= t; i++)
		{
			int u, v;
			double w;
			scanf("%d %d %lf", &u, &v, &w);
			AddEdge(u, v, w);
		}
		
		Dijkstra();
		double ans = d[1];
		int p1 = 1, p2 = 0;
		int flag = 1;
		for(int i = 1; i <= n; i++)
		{
			if(ans < d[i])
			{
				ans = d[i];
				p1 = i;
			}
		}
		for(int i = 1; i <= n; i++)
		{
			for(int j = 0; j < G[i].size(); j++)
			{
				Edge e = edges[G[i][j]];
				double s = (d[i] + d[e.to] + e.dist) / 2;
				if(s > ans)
				{
					flag = 0;
					ans = s;
					if(i < e.to)
					{
						p1 = i;
						p2 = e.to;
					}
					else
					{
						p2 = i;
						p1 = e.to;
					}
				}
			}
		}
		if(flag)
		{
			printf("System #%d\nThe last domino falls after %.1f seconds, at key domino %d.\n\n", cas++, ans, p1);
		}
		else
		{
			printf("System #%d\nThe last domino falls after %.1f seconds, between key dominoes %d and %d.\n\n", cas++, ans, p1, p2);
		}
	}
	return 0;
}

Bellmford

#include <cstring>
#include <cstdio>
#include <vector>
#include <queue>
using namespace std;
const int maxn = 510;
const double INF = 999999999;

struct HeapNode
{
	double d;
	int u;
	bool operator < (const HeapNode& rhs) const
	{
		return d > rhs.d;
	}
};

struct Edge
{
	int from, to;
	double dist;
};
vector <Edge> edges;
vector <int> G[maxn];
bool done[maxn];
double d[maxn];
int cnt[maxn];
bool inq[maxn];
int n, m, t;

void AddEdge(int from, int to, double dist)
{
	edges.push_back((Edge){from, to, dist});
	edges.push_back((Edge){to, from, dist});
	m = edges.size();
	G[from].push_back(m-2);
	G[to].push_back(m-1);
}

bool Bellmanford()
{
	queue <int> Q;
	memset(inq, 0, sizeof(inq));
	memset(cnt, 0, sizeof(cnt));
	for(int i = 1; i <= n; i++)
		d[i] = INF;
	d[1] = 0;
	Q.push(1);
	while(!Q.empty())
	{
		int u = Q.front(); Q.pop();
		inq[u] = false;
		for(int i = 0; i < G[u].size(); i++)
		{
			Edge& e = edges[G[u][i]];
			if(d[e.to] > d[u] + e.dist)
			{
				d[e.to] = d[u] + e.dist;
				
				if(!inq[e.to])
				{
					Q.push(e.to);
					inq[e.to] = true;
					if(++cnt[e.to] > n)
						return false;
				}
			}
		}
	}
}
int main()
{
	int cas = 1;
	while(scanf("%d %d", &n, &t) && (n+t))
	{
		edges.clear();
		for(int i = 1;i <= n; i++)
			G[i].clear();
		for(int i = 1; i <= t; i++)
		{
			int u, v;
			double w;
			scanf("%d %d %lf", &u, &v, &w);
			AddEdge(u, v, w);
		}
		
		Bellmanford();
		double ans = d[1];
		int p1 = 1, p2 = 0;
		int flag = 1;
		for(int i = 1; i <= n; i++)
		{
			if(ans < d[i])
			{
				ans = d[i];
				p1 = i;
			}
		}
		for(int i = 1; i <= n; i++)
		{
			for(int j = 0; j < G[i].size(); j++)
			{
				Edge e = edges[G[i][j]];
				double s = (d[i] + d[e.to] + e.dist) / 2;
				if(s > ans)
				{
					flag = 0;
					ans = s;
					if(i < e.to)
					{
						p1 = i;
						p2 = e.to;
					}
					else
					{
						p2 = i;
						p1 = e.to;
					}
				}
			}
		}
		if(flag)
		{
			printf("System #%d\nThe last domino falls after %.1f seconds, at key domino %d.\n\n", cas++, ans, p1);
		}
		else
		{
			printf("System #%d\nThe last domino falls after %.1f seconds, between key dominoes %d and %d.\n\n", cas++, ans, p1, p2);
		}
	}
	return 0;
}



 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值