20180520模拟赛T1——math

【问题描述】

小美有 n 个点 m 条边。

让你给每个点一个正整数编号。

每条边有两个属性，相连的两个点的编号的 GCD 和 LCM。

题目保证整张图连通。

让你构造出一个编号。

【输入格式】

从文件 math.in 中读入数据。

第一行两个正整数 n 和 m。

接下去m行每行 4 个正整数 xi,yi,gcdi,lcmi。

【输出格式】

输出到文件 math.out 中。

如果是有解：

第一行一行 YES。

第二行 n 个数字表示编号。

否则输出一行NO。

【样例】

【样例输入】

1 0

【样例输出】

YES 1

【样例输入】

2 1
1 2 1 3

【样例输出】

YES 1 3

【样例输入】

3 2
3 2 1 2
3 1 1 10

【样例输出】

YES 5 1 2

【样例输入】

2 1
1 2 3 7

【样例输出】

NO

【数据规模】

对于\(100\%\)的数据\(2\le n \le 100,n-1 \le m \le n*(n-1)/2,1\le gcdi, lcmi \le 10^6\)。

题解

这题其实就是一道暴力题，是根据\([a, b] \times (a, b) = a\times b\)。已知\(a\times b\)，暴力枚举\(a\)，然后把整张图dfs一遍判断其正确性即可。

不得不说：大力出奇迹!

另外，记得开long long。

题目加强版：CF 60C。

代码

#include <cctype>
#include <cstdio>
#include <cstring>

typedef long long LL;
#define int long long

#define dd c = getchar()
inline void read(int& x)
{
	x = 0;
	char dd;
	bool f = false;
	for(; !isdigit(c); dd)
		if(c == '-')
			x = -x;
	for(; isdigit(c); dd)
		x = (x<<1) + (x<<3) + (c^48);
	if(f) x = -x;
}
#undef dd

inline int gcd(int __n, int __m)
{
	while (__n)
	{
		int __t = __m % __n;
		__m = __n;
		__n = __t;
	}
	return __m;
}

const int maxn = 105;

int n, m;

int ans[maxn];

struct edge
{
	LL cheng;
	int t;
	int ne;
	int gcdd;
} e[maxn*maxn];

int first[maxn];
bool vis[maxn];

int mm;
inline void add_edge(int f, int t, LL cheng, int gcdd) 
{
	e[++mm].ne = first[f];
	e[mm].t = t;
	e[mm].cheng = cheng;
	e[mm].gcdd = gcdd;
	first[f] = mm;

	e[++mm].ne = first[t];
	e[mm].t = f;
	e[mm].cheng = cheng;
	e[mm].gcdd = gcdd;
	first[t] = mm;
}

inline bool dfs(int n, int last)//dfs检验答案正确性 
{
	vis[n] = true;
	for(int i = first[n]; i; i = e[i].ne)
	{
		int to = e[i].t;
		if(to == last) continue;
		if(vis[to])
		{
			if((ans[to]*ans[n] != e[i].cheng) || (gcd(ans[to], ans[n]) != e[i].gcdd))
				return false;
		}
		else
		{
			ans[to] = e[i].cheng / ans[n];
			if(gcd(ans[to], ans[n]) != e[i].gcdd)
				return false;
			if(!dfs(to, n))
				return false;
		}
	}
	return true;
}

inline void print()
{
	puts("YES");
	for(int i = 1; i <= n; ++i)
		printf("%lld ", ans[i]);
}

inline void search()
{
	vis[1] = true;
	if(!first[1])
	{
		puts("NO");
		return;
	}
	LL k = e[first[1]].cheng;
	for(int i = 1; i*i <= k; ++i)//暴力枚举所有可能情况 
		if(!(k%i))
		{
			ans[1] = i;
			memset(vis, 0, sizeof(vis));
			if(dfs(1, 0))
			{
				print();
				return;
			}
			ans[1] = k/i;
			memset(vis, 0, sizeof(vis));
			if(dfs(1, 0))
			{
				print();
				return;
			}
		}
	puts("NO");
}

signed main()
{
	freopen("math.in", "r", stdin);
	freopen("math.out", "w", stdout);
	scanf("%lld%lld", &n, &m);
	if(n == 1)
	{
		puts("YES\n1");
		return 0;
	}
	for(int i = 1; i <= m; ++i)
	{
		int f, t, gcdd, lcmm;
		read(f), read(t), read(gcdd), read(lcmm);
		add_edge(f, t, (LL)(gcdd*lcmm), gcdd);
	}
	search();
	return 0;
}