ダイクストラ+チェーンフォワードスター+優先キューの最適化

招待状POJ-1511

説明
n-1人の個人がポイント1から残りのn-1プロパガンダポイントまで開始し、次にポイント1に戻って結果を報告し、すべての往復パスの最小合計を見つけます

入力
入力はTケースで構成されます。入力の最初の行には、正の整数Tのみが含まれています。
次はNとM、1 <= N、M <= 1000000です。これは、N個のポイントとN個のポイントを接続するM個のエッジを意味します。
次に、M行があり、各行には3つの値U、V、Wが含まれています。これは、ポイントUからポイントVまでの距離にWが必要であることを示しています。グラフが接続されていると想定できます。
これは一方向のチャネルであることに注意してください。

出力
それぞれの場合について、合計パスの最小値を示す行を印刷します。

サンプル入力
2
2 2
1 2 13
2 1 33
4 6
1 2 10
2 1 60
1 3 20
3 4 10
2 4 5
4 1 50

サンプル出力
46
210

**考え:**この質問では、最短の往復パスの最小値を見つけることができます。ダイクストラを使用して、順方向と逆方向の両方を見つけ、それらを合計することができます。
**注:**質問のデータ範囲が大きすぎます。隣接行列を使用してグラフを保存すると、間違いなくTになるため、チェーンフォワードスターストレージマップを使用することにしました(チェーンフォワードスターがそうであることがわかっただけです)。激しい)しかし、私たちが予期していなかったのは---、ダイクストラは優先キュー最適化またはTを使用せず、それはずさんで、後で優先キュー+チェーンフォワードスターパスを使用しました--AC
コードは次のとおりです:

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cctype>
#include<iomanip>
#include<map>
#include<vector>
#include<list>
#include<deque>
#include<stack>
#include<queue>
#include<set>
#include<cctype>
#include<string>
#include<stdexcept>
#include<fstream>
#include<sstream>
#define mem(a,b) memset(a,b,sizeof(a))
#define mod 10000007
#define debug() puts("what the fuck!")
#define dedebug() puts("what the fuck!!!")
#define ll long long
#define ull unsigned long long
#define speed {
    
    ios::sync_with_stdio(false); cin.tie(0); cout.tie(0); };
using namespace std;
const double PI = acos(-1.0);
const int maxn = 1e6 + 50;
const int N = 35;
const int INF = 0x3f3f3f3f;
const int inf = 0xfffffff;//比INF小,防止累加爆int
const double esp_0 = 1e-6;
const double gold = (1 + sqrt(5)) / 2;
int gcd(int x, int y) {
    
    
	return y ? gcd(y, x % y) : x;
}
int n, m;
int flag;
struct node {
    
    
	int next, to, w;
}edge[maxn],redge[maxn];
int head[maxn], rhead[maxn], cnt, rcnt;
int dis[maxn], rdis[maxn], vis[maxn];
void add(int u, int v, int w) {
    
    
	edge[cnt].to = v;
	edge[cnt].next = head[u];
	edge[cnt].w = w;
	head[u] = cnt++;
}
void radd(int u, int v, int w) {
    
    
	redge[rcnt].to = v;
	redge[rcnt].next = rhead[u];
	redge[rcnt].w = w;
	rhead[u] = rcnt++;
}
void init() {
    
    
	mem(dis, INF);
	mem(rdis, INF);
	mem(head, -1);
	mem(rhead, -1);
}
struct que {
    
    
	int id, dist;
	que() {
    
    };
	que(int id, int dist) :id(id), dist(dist) {
    
    };
	bool operator<(const que& a)const{
    
    
		return dist > a.dist;
	}

};
void dijkstra(int dis[],int head[],node edge[]) {
    
    
	mem(vis, 0);
	dis[1] = 0;
	priority_queue<que>q;
	q.push(que(1, dis[1]));
	while (!q.empty()) {
    
    
		que step = q.top();
		q.pop();
		if (vis[step.id])continue;
		vis[step.id] = 1;
		for (int i = head[step.id]; i != -1; i = edge[i].next) {
    
    
			int value = edge[i].to;
			if (!vis[value] && dis[value] > dis[step.id] + edge[i].w) {
    
    
				dis[value] = dis[step.id] + edge[i].w;
				q.push(que(value, dis[value]));
			}
		}
	}
}
int main() {
    
    
	int t;
	scanf("%d", &t);
	while (t--) {
    
    
		init();
		scanf("%d%d", &n, &m);
		for (int i = 1; i <= m; ++i) {
    
    
			int u, v, w;
			scanf("%d%d%d", &u, &v, &w);
			add(u, v, w);
			radd(v, u, w);
		}
		dijkstra(dis,head,edge);
		dijkstra(rdis, rhead,redge);
		ll ans = 0;
		for (int i = 1; i <= n; ++i) {
    
    
//			printf("%d %d\n", dis[i], rdis[i]);
			ans += dis[i] + rdis[i];
		}
		printf("%lld\n", ans);
	}
	return 0;
}

優先キューで最適化されていないエラーコード

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cctype>
#include<iomanip>
#include<map>
#include<vector>
#include<list>
#include<deque>
#include<stack>
#include<queue>
#include<set>
#include<cctype>
#include<string>
#include<stdexcept>
#include<fstream>
#include<sstream>
#define mem(a,b) memset(a,b,sizeof(a))
#define mod 10000007
#define debug() puts("what the fuck!")
#define dedebug() puts("what the fuck!!!")
#define ll long long
#define ull unsigned long long
#define speed {
    
    ios::sync_with_stdio(false); cin.tie(0); cout.tie(0); };
using namespace std;
const double PI = acos(-1.0);
const int maxn = 1e6 + 50;
const int N = 35;
const int INF = 0x3f3f3f3f;
const int inf = 0xfffffff;//比INF小,防止累加爆int
const double esp_0 = 1e-6;
const double gold = (1 + sqrt(5)) / 2;
int gcd(int x, int y) {
    
    
	return y ? gcd(y, x % y) : x;
}
int n, m;
struct node {
    
    
	int next, to, w;
}edge[maxn],redge[maxn];
int head[maxn], rhead[maxn], cnt, rcnt;
int dis[maxn], rdis[maxn], vis[maxn];
void add(int u, int v, int w) {
    
    
	edge[cnt].to = v;
	edge[cnt].next = head[u];
	edge[cnt].w = w;
	head[u] = cnt++;
}
void radd(int u, int v, int w) {
    
    
	redge[rcnt].to = v;
	redge[rcnt].next = rhead[u];
	redge[rcnt].w = w;
	rhead[u] = rcnt++;
}
void init() {
    
    
	mem(dis, INF);
	mem(rdis, INF);
	mem(head, -1);
	mem(rhead, -1);
}
void dijkstra(int dis[],int head[],node edge[]) {
    
    
	mem(vis, 0);
	dis[1] = 0;
	for (int i = 1; i <= n; ++i) {
    
    
		int min = INF;
		int step = -1;
		for (int j = 1; j <= n; ++j) {
    
    
			if (!vis[j] && min > dis[j]) {
    
    
				min = dis[j];
				step = j;
			}
		}
		if (step == -1)break;
		vis[step] = 1;
		for (int i = head[step]; i != -1; i = edge[i].next) {
    
    
			int now = edge[i].to;
			if (dis[now] > dis[step] + edge[i].w) {
    
    
				dis[now] = dis[step] + edge[i].w;
			}
		}
	}
}
int main() {
    
    
	int t;
	scanf("%d", &t);
	while (t--) {
    
    
		init();
		scanf("%d%d", &n, &m);
		for (int i = 1; i <= m; ++i) {
    
    
			int u, v, w;
			scanf("%d%d%d", &u, &v, &w);
			add(u, v, w);
			radd(v, u, w);
		}
		dijkstra(dis,head,edge);
		dijkstra(rdis, rhead,redge);
		ll ans = 0;
		for (int i = 1; i <= n; ++i) {
    
    
//			printf("%d %d\n", dis[i], rdis[i]);
			ans += dis[i] + rdis[i];
		}
		printf("%lld\n", ans);
	}
	return 0;
}

おすすめ

転載: blog.csdn.net/qq_40924271/article/details/107763277