【noip模拟题】藏宝图(prim)

好神的一题。。

一开始没想多久就看题解了QAQ。。

首先我们发现,这棵树任意两个点的边一定是最小的(即所有其它这两个点的路径都比这条边大,才有可能出解)

然后生成树后再算距离判断即可。。

注意特判n=1.。。。。。。。。。。。。。。。。。

还有出题人一点都不良心。。。说好的0<=dist[i][j]<=10^9呢。。。搞得我爆了好几次。。。。。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <vector>
#include <map>
using namespace std;
typedef long long ll;
#define pii pair<int, int>
#define mkpii make_pair<int, int>
#define pdi pair<double, int>
#define mkpdi make_pair<double, int>
#define pli pair<ll, int>
#define mkpli make_pair<ll, int>
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define error(x) (!(x)?puts("error"):0)
#define printarr2(a, b, c) for1(_, 1, b) { for1(__, 1, c) cout << a[_][__]; cout << endl; }
#define printarr1(a, b) for1(_, 1, b) cout << a[_] << ‘\t‘; cout << endl
inline const ll getint() { ll r=0, k=1; char c=getchar(); for(; c<‘0‘||c>‘9‘; c=getchar()) if(c==‘-‘) k=-1; for(; c>=‘0‘&&c<=‘9‘; c=getchar()) r=r*10+c-‘0‘; return k*r; }
inline const int max(const int &a, const int &b) { return a>b?a:b; }
inline const int min(const int &a, const int &b) { return a<b?a:b; }

const int N=2505, oo=~0u>>1;
int n, vis[N], from[N], ihead[N], cnt;
ll sum[N], child[N], mp[N][N], mp2[N][N], d[N];
struct dat { int next, to; ll w; }e[N<<1];
priority_queue<pli, vector<pli>, greater<pli> >q;
void add(int u, int v, ll w) {
	e[++cnt].next=ihead[u]; ihead[u]=cnt; e[cnt].to=v; e[cnt].w=w;
	e[++cnt].next=ihead[v]; ihead[v]=cnt; e[cnt].to=u; e[cnt].w=w;
}
void prim() {
	for1(i, 1, n) d[i]=oo;
	d[1]=0; q.push(mkpli(0, 1));
	while(!q.empty()) {
		int u=q.top().second; q.pop(); if(vis[u]) continue; vis[u]=1;
		for1(v, 1, n) if(!vis[v] && mp[u][v] && d[v]>mp[u][v]) {
			d[v]=mp[u][v];
			from[v]=u;
			q.push(mkpli(d[v], v));
		}
	}
	for1(i, 2, n) if(d[i]!=oo) add(from[i], i, d[i]);
}
void dfs(int x, int fa) {
	for1(i, 1, n) if(vis[i]) mp2[i][x]=mp2[x][i]=mp2[i][fa]+mp2[fa][x];
	int y; child[x]=0; sum[x]=0; vis[x]=1;
	for(int i=ihead[x]; i; i=e[i].next) if((y=e[i].to)!=fa) {
		mp2[x][y]=mp2[y][x]=e[i].w;
		dfs(y, x);
		sum[y]+=e[i].w;
		sum[x]+=e[i].w;
		child[y]++;
		child[x]++;
	}
}
void init() {
	for1(i, 1, n) vis[i]=ihead[i]=from[i]=child[i]=sum[i]=0;
	for1(i, 1, n) for1(j, 1, n) mp2[i][j]=mp[i][j]=0;
	cnt=0;
	while(!q.empty()) q.pop();
}
int main() {
	int cs=getint();
	while(cs--) {
		read(n); if(n==1) { puts("Yes"); puts("1"); n=getint(); continue; }
		init();
		int flag=0;

		for1(i, 1, n) for1(j, 1, n) read(mp[i][j]);
		prim();
		for1(i, 1, n) if(d[i]==oo) { puts("No"); flag=1; break; }
		if(flag) continue;

		for1(i, 1, n) vis[i]=0;
		dfs(1, 0);

		//for1(i, 1, n) { for1(j, 1, n) printf("%d ", mp2[i][j]); puts(""); }
		for1(i, 1, n) { for1(j, 1, n) if(mp2[i][j]!=mp[i][j]) { flag=1; break; } if(flag) break; }
		if(!flag) {
			double mx=-1e199; int ans;
			// for1(i, 1, n) dbg(sum[i]);
			// for1(i, 1, n) dbg(child[i]);
			for1(i, 1, n) if((double)sum[i]/(double)child[i]>mx) { mx=sum[i]/child[i]; ans=i; }
			printf("Yes\n%d\n", ans);
		}
		else puts("No");
	}
	return 0;
}

  



背景

Czy爬上黑红树,到达了一个奇怪的地方……

题目描述

Czy发现了一张奇怪的藏宝图。图上有n个点,m条无向边。已经标出了图中两两之间距离dist。但是czy知道,只有当图刚好又是一颗树的时候,这张藏宝图才是真的。如果藏宝图是真的,那么经过点x的边的边权平均数最大的那个x是藏着宝物的地方。请计算这是不是真的藏宝图,如果是真的藏宝之处在哪里。

格式

输入数据第一行一个数T,表示T组数据。

对于每组数据,第一行一个n,表示藏宝图上的点的个数。

接下来n行,每行n个数,表示两两节点之间的距离。

输出一行或两行。第一行”Yes”或”No”,表示这是不是真的藏宝图。

若是真的藏宝图,第二行再输出一个数,表示哪个点是藏宝之处。

样例输入

2

3

0 7 9

7 0 2

9 2 0

3

0 2 7

2 0 9

7 9 0

样例输出

Yes

1

Yes

3

样例解释:第一棵树的形状是1--2--3。1、2之间的边权是7,2、3之间是2。

第二棵树的形状是2--1--3。2、1之间的边权是2,1、3之间是7。

数据范围

对于30%数据,n<=50,1<=树上的边的长度<=10^9。

对于50%数据,n<=600.

对于100%数据,1<=n<=2500,除30%小数据外任意0<=dist[i][j]<=10^9,T<=5

时间: 2024-10-08 01:13:57

【noip模拟题】藏宝图(prim)的相关文章

[BZOJ入门OJ2092][Noip模拟题]舞会

2092: [Noip模拟题]舞会 Time Limit: 20 Sec  Memory Limit: 256 MB Submit: 9  Solved: 5 [Submit][Status][Web Board] Description 学校举行舞会啦,一共有N个人参加,所有人站成一排,从左开始编号,最左边的人编号为1 ,最右边的为N.每个人跳舞的熟练度我们用一个整数表示,第i个人的熟练度为Ai,每次熟 练度最接近的一对相邻男女会出列跳舞,如果有多对那么最左边的那一对会先出列,请你给 出出列跳

8.22 NOIP 模拟题

  8.22 NOIP 模拟题 编译命令 g++ -o * *.cpp gcc -o * *.c fpc *.pas 编译器版本 g++/gcc 4.9.2 fpc 2.6.2 评测环境 64 位 Linux, 3.3GHZ CPU 评测软件 Lemon 评测方式 忽略行末空格和回车 特别注意:c/c++ 选手使用 printf 输出 64 位整数请使用%lld 1 注意事项 A 债务 文件名                            输入文件             输出文件  

noip模拟题题解集

最近做模拟题看到一些好的题及题解. 升格思想: 核电站问题 一个核电站有N个放核物质的坑,坑排列在一条直线上.如果连续M个坑中放入核物质,则会发生爆炸,于是,在某些坑中可能不放核物质. 任务:对于给定的N和M,求不发生爆炸的放置核物质的方案总数 输入:输入文件只一行,两个正整数N,M( 1<N<50,2≤M≤5) 输出:输出文件只有一个正整数S,表示方案总数. 运用升格思想.设N个坑不会发生爆炸的方案数是f[N],那么我们假设N以前的坑的方案 都已知了,那么我们只需要考虑第N个坑如何放即可(顺

9.22 NOIP模拟题

吉林省信息学奥赛 2017 冬令营                                                                                                                                                                                                      day2 嗯 今天题目好水 嗯 今天爆零了好伤心...... 值得纪念

6.19 noip模拟题(题目及解析转自 hzwer 2014-3-15 NOIP模拟赛)

Problem 1 高级打字机(type.cpp/c/pas) [题目描述] 早苗入手了最新的高级打字机.最新款自然有着与以往不同的功能,那就是它具备撤销功能,厉害吧. 请为这种高级打字机设计一个程序,支持如下3种操作: 1.T x:在文章末尾打下一个小写字母x.(type操作) 2.U x:撤销最后的x次修改操作.(Undo操作) (注意Query操作并不算修改操作) 3.Q x:询问当前文章中第x个字母并输出.(Query操作) 文章一开始可以视为空串. [输入格式] 第1行:一个整数n,表

【noip模拟题】天神下凡(贪心)

vijos某次模拟赛原题... 处理出每个圆的一级祖先就行了... 其实没有那么麻烦,贪心即可出解. 我们将每个圆转换成线段后按左端点小右端点大的方法排序 然后维护一个栈: 对于每一个圆i 如果栈顶右端点比圆i的右端点小,则出栈,直到栈空 否则i的一级祖先就是栈顶,并且加入i到栈. 证明: 因为左端点排序,所以问题转换为找一个最小的右端点能够包含此线段 假如栈顶的右端点比当前右端点小,显然对于所有将来的线段,不可能包含将来的线段(或者说,尽管能,也不是最优解,因为最优解就是当前i) 然后如果有n

【noip模拟题】挖掘机(模拟题+精度)

这题直接模拟. 可是我挂在了最后两个点上QAQ.唯一注意的是注意精度啊...用来double后边转成整数就忘记用longlong...sad #include <cstdio> #include <cstring> #include <cmath> #include <string> #include <iostream> #include <algorithm> #include <queue> #include <

【noip模拟题】迎接仪式(dp+特殊的技巧)

好神的一题... 这是一道DP题,本题的难点在于状态的确定,由于调整是任意的,很难划分状态,我们略微修改一下调整的形式:把一次’j’和’z’交换看做两次变换:’j’->’z’;’z’->’j’ (zz交换和jj交换是没有意义的,不作考虑);于是最多’j’->'z' 'z'->'j'各K次. F[i,j,k]:=F[i-2,j-x,k-y]+1;If A[i-1]=j then x=0 else x=1If A[i]=z then y=0 else y=1 然后最后要看是否有那么多个

XJOI NOIP模拟题1

第一题 分析: 开始想的是贪心,取每列均值最大一段. 应该是01分数规划,具体看代码 代码: program gold; var a:array[0..100000]of int64; n,i,m,j,x:longint; function max(x,y:real):real; begin if x>y then max:=x else max:=y; end; function cheak(x:real):boolean; var i,j:longint; s,ans:real; begin