結果

問題 No.1776 Love Triangle 2 (Hard)
ユーザー 👑 ygussanyygussany
提出日時 2021-11-19 18:49:53
言語 C
(gcc 12.3.0)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 3,813 bytes
コンパイル時間 1,130 ms
コンパイル使用メモリ 32,504 KB
実行使用メモリ 8,760 KB
最終ジャッジ日時 2023-09-15 07:37:45
合計ジャッジ時間 13,113 ms
ジャッジサーバーID
(参考情報)
judge12 / judge13
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 0 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 0 ms
4,380 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 TLE -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
testcase_43 -- -
testcase_44 -- -
testcase_45 -- -
testcase_46 -- -
testcase_47 -- -
testcase_48 -- -
testcase_49 -- -
testcase_50 -- -
testcase_51 -- -
testcase_52 -- -
testcase_53 -- -
testcase_54 -- -
testcase_55 -- -
testcase_56 -- -
testcase_57 -- -
testcase_58 -- -
testcase_59 -- -
testcase_60 -- -
testcase_61 -- -
testcase_62 -- -
testcase_63 -- -
testcase_64 -- -
testcase_65 -- -
testcase_66 -- -
testcase_67 -- -
testcase_68 -- -
testcase_69 -- -
testcase_70 -- -
testcase_71 -- -
testcase_72 -- -
testcase_73 -- -
testcase_74 -- -
testcase_75 -- -
testcase_76 -- -
testcase_77 -- -
testcase_78 -- -
testcase_79 -- -
testcase_80 -- -
testcase_81 -- -
testcase_82 -- -
testcase_83 -- -
testcase_84 -- -
testcase_85 -- -
testcase_86 -- -
testcase_87 -- -
testcase_88 -- -
testcase_89 -- -
testcase_90 -- -
testcase_91 -- -
testcase_92 -- -
testcase_93 -- -
testcase_94 -- -
testcase_95 -- -
testcase_96 -- -
testcase_97 -- -
testcase_98 -- -
testcase_99 -- -
testcase_100 -- -
testcase_101 -- -
testcase_102 -- -
testcase_103 -- -
testcase_104 -- -
testcase_105 -- -
testcase_106 -- -
testcase_107 -- -
testcase_108 -- -
testcase_109 -- -
testcase_110 -- -
testcase_111 -- -
testcase_112 -- -
testcase_113 -- -
testcase_114 -- -
testcase_115 -- -
testcase_116 -- -
testcase_117 -- -
testcase_118 -- -
testcase_119 -- -
testcase_120 -- -
testcase_121 -- -
testcase_122 -- -
testcase_123 -- -
testcase_124 -- -
testcase_125 -- -
testcase_126 -- -
testcase_127 -- -
testcase_128 -- -
testcase_129 -- -
testcase_130 -- -
testcase_131 -- -
testcase_132 -- -
testcase_133 -- -
testcase_134 -- -
testcase_135 -- -
testcase_136 -- -
testcase_137 -- -
testcase_138 -- -
testcase_139 -- -
testcase_140 -- -
testcase_141 -- -
testcase_142 -- -
testcase_143 -- -
testcase_144 -- -
testcase_145 -- -
testcase_146 -- -
testcase_147 -- -
testcase_148 -- -
testcase_149 -- -
testcase_150 -- -
testcase_151 -- -
testcase_152 -- -
testcase_153 -- -
testcase_154 -- -
testcase_155 -- -
testcase_156 -- -
testcase_157 -- -
testcase_158 -- -
testcase_159 -- -
testcase_160 -- -
testcase_161 -- -
testcase_162 -- -
testcase_163 -- -
testcase_164 -- -
testcase_165 -- -
testcase_166 -- -
testcase_167 -- -
testcase_168 -- -
testcase_169 -- -
testcase_170 -- -
testcase_171 -- -
testcase_172 -- -
testcase_173 -- -
testcase_174 -- -
testcase_175 -- -
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <stdio.h>

#define N_MAX 150
#define M_MAX 12000



typedef struct Edge {
	struct Edge *next;
	int v;
	unsigned int label;
} edge;

void chmin(int* a, int b)
{
	if (*a > b) *a = b;
}

int lex_smaller(int a[], int b[])
{
	int i;
	for (i = 0; i <= a[0]; i++) {
		if (a[i] < b[i]) return 1;
		else if (a[i] > b[i]) return -1;
	}
	return 0;
}

void chlexmin(int a[], int b[])
{
	int i;
	if (lex_smaller(a, b) < 0) for (i = 0; i <= b[0]; i++) a[i] = b[i];
}

void print_ans(int N, int ans[])
{
	int i;
	if (ans[0] > N) {
		printf("-1\n");
		return;
	} else printf("%d\n", ans[0]);
	for (i = 1; i <= ans[0]; i++) printf("%d ", ans[i]);
	printf("%d\n", ans[1]);
}



void DFS_naive_lexmin(edge* adj[], int X, int Y, int Z, int u, int flag[], int tmp[], int ans[])
{
	int i;
	if (flag[u] != 0) {
		if (u == X && flag[Y] != 0 && flag[Z] != 0) for (i = 0; i <= tmp[0]; i++) ans[i] = tmp[i];
		return;
	} else flag[u] = 1;
	
	edge *p;
	tmp[++tmp[0]] = u;
	if (lex_smaller(tmp, ans) > 0) for (p = adj[u]; p != NULL; p = p->next) DFS_naive_lexmin(adj, X, Y, Z, p->v, flag, tmp, ans);
	tmp[0]--;
	flag[u] = 0;
}

// Naive solution for finding the lexmin solution
int naive_lexmin(int N, int M, int X, int Y, int Z, int A[], int B[], int ans[])
{
	static int i, u, w, adj_mat[N_MAX + 1][N_MAX + 1];
	static edge *adj[N_MAX + 1], e[M_MAX * 2], *p;
	for (u = 1; u <= N; u++) for (w = u + 1; w <= N; w++) adj_mat[u][w] = 0;
	for (i = 1; i <= M; i++) {
		u = A[i];
		w = B[i];
		adj_mat[u][w] = 1;
	}
	for (u = 1; u <= N; u++) adj[u] = NULL;
	for (u = 1, i = 0; u <= N; u++) {
		for (w = N; w > u; w--) {
			if (adj_mat[u][w] != 0) continue;
			e[i].v = w;
			e[i].next = adj[u];
			adj[u] = &(e[i++]);
		}
		for (w = u - 1; w >= 1; w--) {
			if (adj_mat[w][u] != 0) continue;
			e[i].v = w;
			e[i].next = adj[u];
			adj[u] = &(e[i++]);
		}
	}
	
	static int tmp[N_MAX + 1], flag[N_MAX + 1] = {}, q[N_MAX + 1], head, tail;
	
	// Y--Z path?
	flag[X] = 1;
	flag[Y] = 1;
	q[0] = X;
	q[1] = Y;
	for (head = 1, tail = 2; head < tail; head++) {
		u = q[head];
		for (p = adj[u]; p != NULL; p = p->next) {
			w = p->v;
			if (flag[w] == 0) {
				flag[w] = 1;
				q[tail++] = w;
			}
		}
	}
	if (flag[Z] == 0) {
		for (head = 0; head < tail; head++) flag[q[head]] = 0;
		ans[0] = N + 1;
		for (i = 1; i <= N + 1; i++) ans[i] = 0;
		return -1;
	} else for (head = 0; head < tail; head++) flag[q[head]] = 0;
	
	// Z--X path?
	flag[Y] = 1;
	flag[Z] = 1;
	q[0] = Y;
	q[1] = Z;
	for (head = 1, tail = 2; head < tail; head++) {
		u = q[head];
		for (p = adj[u]; p != NULL; p = p->next) {
			w = p->v;
			if (flag[w] == 0) {
				flag[w] = 1;
				q[tail++] = w;
			}
		}
	}
	if (flag[X] == 0) {
		for (head = 0; head < tail; head++) flag[q[head]] = 0;
		ans[0] = N + 1;
		for (i = 1; i <= N + 1; i++) ans[i] = 0;
		return -1;
	} else for (head = 0; head < tail; head++) flag[q[head]] = 0;
	
	// X--Y path?
	flag[Z] = 1;
	flag[X] = 1;
	q[0] = Z;
	q[1] = X;
	for (head = 1, tail = 2; head < tail; head++) {
		u = q[head];
		for (p = adj[u]; p != NULL; p = p->next) {
			w = p->v;
			if (flag[w] == 0) {
				flag[w] = 1;
				q[tail++] = w;
			}
		}
	}
	if (flag[Y] == 0) {
		for (head = 0; head < tail; head++) flag[q[head]] = 0;
		ans[0] = N + 1;
		for (i = 1; i <= N + 1; i++) ans[i] = 0;
		return -1;
	} else for (head = 0; head < tail; head++) flag[q[head]] = 0;
	
	// DFS
	tmp[0] = 0;
	ans[0] = N + 1;
	for (i = 1; i <= N + 1; i++) ans[i] = 0;
	DFS_naive_lexmin(adj, X, Y, Z, X, flag, tmp, ans);
	if (ans[0] > N) return -1;
	else return ans[0];
}



int main()
{
	int i, N, M, X, Y, Z, A[M_MAX + 1], B[M_MAX + 1], ans[N_MAX + 2];
	scanf("%d %d", &N, &M);
	scanf("%d %d %d", &X, &Y, &Z);
	for (i = 1; i <= M; i++) scanf("%d %d", &(A[i]), &(B[i]));
	naive_lexmin(N, M, X, Y, Z, A, B, ans);
	print_ans(N, ans);
	fflush(stdout);
	return 0;
}
0