week10
Einstein学画画
题目描述
Einstein 学起了画画。
此人比较懒~~,他希望用最少的笔画画出一张画……
给定一个无向图,包含 n n n 个顶点(编号 1 ∼ n 1 \sim n 1∼n), m m m 条边,求最少用多少笔可以画出图中所有的边。
输入格式
第一行两个整数 n , m n, m n,m。
接下来 m m m 行,每行两个数 a , b a, b a,b( a ≠ b a \ne b a=b),表示 a , b a, b a,b 两点之间有一条边相连。
一条边不会被描述多次。
输出格式
一个数,即问题的答案。
样例 #1
样例输入 #1
5 5
2 3
2 4
2 5
3 4
4 5
样例输出 #1
1
提示
对于 50 % 50 \% 50% 的数据, n ≤ 50 n \le 50 n≤50, m ≤ 100 m \le 100 m≤100。
对于 100 % 100\% 100% 的数据, 1 ≤ n ≤ 1000 1 \le n \le 1000 1≤n≤1000, 1 ≤ m ≤ 10 5 1 \le m \le {10}^5 1≤m≤105。
解题方法
此题的解题方法为欧拉路,则可知因为一个连通图只可能有偶数个奇点,那么我们可以从欧拉路的介绍,发现只有2个奇点才能1笔画完,而奇点是指跟这个点相连的边数目有奇数个的点
定理1:存在欧拉路的条件:图是连通的,有且只有2个奇点。
定理2:存在欧拉回路的条件:图是连通的,有0个奇点。
所以我们直接判断某个点的度,是否为奇数就能解决。每多出两个奇点,那么画的次数就会+1,又因为连通图只可能有偶数个奇点所以,答案为奇点数除2
解题代码
#include<iostream>
#include<cstring>
using namespace std;
int main()
{
int n,m;
int num[1010],ans=0;
memset(num,0,sizeof(num));
cin>>n>>m;
for(int i=1;i<=m;i++)
{
int a,b;
cin>>a>>b;
num[a]++;//点的度+1
num[b]++;
}
for(int i=1;i<=n;i++)
{
if(num[i]%2==1) ans++;//计算奇点个数
}
if(!ans) cout<<1;//判断特殊条件
else cout<<ans/2;//每有两个奇点就是一个连通图
return 0;
}
[蓝桥杯 2017 国 C] 合根植物
题目描述
w 星球的一个种植园,被分成 m × n m \times n m×n 个小格子(东西方向 m m m 行,南北方向 n n n 列)。每个格子里种了一株合根植物。
这种植物有个特点,它的根可能会沿着南北或东西方向伸展,从而与另一个格子的植物合成为一体。
如果我们告诉你哪些小格子间出现了连根现象,你能说出这个园中一共有多少株合根植物吗?
输入格式
第一行,两个整数 m m m, n n n,用空格分开,表示格子的行数、列数( 1 < m , n < 1000 1<m,n<1000 1<m,n<1000)。
接下来一行,一个整数 k k k,表示下面还有 k k k 行数据 ( 0 < k < 1 0 5 ) (0<k<10^5) (0<k<105)。
接下来 k k k 行,第行两个整数 a a a, b b b,表示编号为 a a a 的小格子和编号为 b b b 的小格子合根了。
格子的编号一行一行,从上到下,从左到右编号。
比如: 5 × 4 5 \times 4 5×4 的小格子,编号:
1 2 3 4
5 6 7 8
9 10 11 12
13 14 15 16
17 18 19 20
输出格式
一行一个整数,表示答案
样例 #1
样例输入 #1
5 4
16
2 3
1 5
5 9
4 8
7 8
9 10
10 11
11 12
10 14
12 16
14 18
17 18
15 19
19 20
9 13
13 17
样例输出 #1
5
提示
样例解释
时限 1 秒, 256M。蓝桥杯 2017 年第八届国赛
解题方法
比较典型的并查集的问题,需要记住基本的并查集的写法,就能轻松解决。
解题代码
#include<iostream>
#include<cstring>
using namespace std;
int root[1000005],ans=0,num[1000005];
int find(int x)//标准的并查集查找“父亲”的函数
{
if(root[x]==x) return x;//如果“父亲”就是自己,则返回自己
return root[x]=find(root[x]);//寻找“父亲”的“父亲”
}
void unite(int a,int b)
{
int root1=find(a);//查找两个点的“父亲”
int root2=find(b);
if(root1!=root2) root[root2]=root1;//如果“父亲”不一样,则归为同一个“父亲”
}
int main()
{
int n,m,k;
cin>>n>>m;
cin>>k;
memset(num,0,sizeof(num));
for(int i=1;i<=n*m;i++)//初始化
{
root[i]=i;
}
for(int i=1;i<=k;i++)
{
int a,b;
cin>>a>>b;
unite(a,b);//合并为同一个“父亲”
}
for(int i=1;i<=n*m;i++)
{
num[find(i)]=1;//标记每个连通图的“父亲”,便于后面的计算
}
for(int i=1;i<=n*m;i++)
{
if(num[i]) ans++;
}
cout<<ans<<endl;
return 0;
}
[NOIP2017 提高组] 奶酪
题目背景
NOIP2017 提高组 D2T1
题目描述
现有一块大奶酪,它的高度为 h h h,它的长度和宽度我们可以认为是无限大的,奶酪中间有许多半径相同的球形空洞。我们可以在这块奶酪中建立空间坐标系,在坐标系中,奶酪的下表面为 z = 0 z = 0 z=0,奶酪的上表面为 z = h z = h z=h。
现在,奶酪的下表面有一只小老鼠 Jerry,它知道奶酪中所有空洞的球心所在的坐标。如果两个空洞相切或是相交,则 Jerry 可以从其中一个空洞跑到另一个空洞,特别地,如果一个空洞与下表面相切或是相交,Jerry 则可以从奶酪下表面跑进空洞;如果一个空洞与上表面相切或是相交,Jerry 则可以从空洞跑到奶酪上表面。
位于奶酪下表面的 Jerry 想知道,在不破坏奶酪的情况下,能否利用已有的空洞跑 到奶酪的上表面去?
空间内两点 P 1 ( x 1 , y 1 , z 1 ) P_1(x_1,y_1,z_1) P1(x1,y1,z1)、 P 2 ( x 2 , y 2 , z 2 ) P2(x_2,y_2,z_2) P2(x2,y2,z2) 的距离公式如下:
d i s t ( P 1 , P 2 ) = ( x 1 − x 2 ) 2 + ( y 1 − y 2 ) 2 + ( z 1 − z 2 ) 2 \mathrm{dist}(P_1,P_2)=\sqrt{(x_1-x_2)^2+(y_1-y_2)^2+(z_1-z_2)^2} dist(P1,P2)=(x1−x2)2+(y1−y2)2+(z1−z2)2
输入格式
每个输入文件包含多组数据。
第一行,包含一个正整数 T T T,代表该输入文件中所含的数据组数。
接下来是 T T T 组数据,每组数据的格式如下: 第一行包含三个正整数 n , h , r n,h,r n,h,r,两个数之间以一个空格分开,分别代表奶酪中空洞的数量,奶酪的高度和空洞的半径。
接下来的 n n n 行,每行包含三个整数 x , y , z x,y,z x,y,z,两个数之间以一个空格分开,表示空洞球心坐标为 ( x , y , z ) (x,y,z) (x,y,z)。
输出格式
T
T
T 行,分别对应
T
T
T 组数据的答案,如果在第
i
i
i 组数据中,Jerry 能从下表面跑到上表面,则输出 Yes
,如果不能,则输出 No
。
样例 #1
样例输入 #1
3
2 4 1
0 0 1
0 0 3
2 5 1
0 0 1
0 0 4
2 5 2
0 0 2
2 0 4
样例输出 #1
Yes
No
Yes
提示
【输入输出样例 1 1 1 说明】
第一组数据,由奶酪的剖面图可见:
第一个空洞在 ( 0 , 0 , 0 ) (0,0,0) (0,0,0) 与下表面相切;
第二个空洞在 ( 0 , 0 , 4 ) (0,0,4) (0,0,4) 与上表面相切;
两个空洞在 ( 0 , 0 , 2 ) (0,0,2) (0,0,2) 相切。
输出 Yes
。
第二组数据,由奶酪的剖面图可见:
两个空洞既不相交也不相切。
输出 No
。
第三组数据,由奶酪的剖面图可见:
两个空洞相交,且与上下表面相切或相交。
输出 Yes
。
【数据规模与约定】
对于 20 % 20\% 20% 的数据, n = 1 n = 1 n=1, 1 ≤ h 1 \le h 1≤h, r ≤ 1 0 4 r \le 10^4 r≤104,坐标的绝对值不超过 1 0 4 10^4 104。
对于 40 % 40\% 40% 的数据, 1 ≤ n ≤ 8 1 \le n \le 8 1≤n≤8, 1 ≤ h 1 \le h 1≤h, r ≤ 1 0 4 r \le 10^4 r≤104,坐标的绝对值不超过 1 0 4 10^4 104。
对于 80 % 80\% 80% 的数据, 1 ≤ n ≤ 1 0 3 1 \le n \le 10^3 1≤n≤103, 1 ≤ h , r ≤ 1 0 4 1 \le h , r \le 10^4 1≤h,r≤104,坐标的绝对值不超过 1 0 4 10^4 104。
对于 100 % 100\% 100% 的数据, 1 ≤ n ≤ 1 × 1 0 3 1 \le n \le 1\times 10^3 1≤n≤1×103, 1 ≤ h , r ≤ 1 0 9 1 \le h , r \le 10^9 1≤h,r≤109, T ≤ 20 T \le 20 T≤20,坐标的绝对值不超过 1 0 9 10^9 109。
解题方法
此题方法可以使用上题的并查集的方法,具体的做法为记录每个与上表面和下表面相交或相切的球,然后把每个球进行并查集的归类。
最后把每个与上表面和下表面相交或相切的球进行并查集查找“父亲”的判断,如果是同一个“父亲”,则答案为可行。
解题代码
#include<iostream>
#include<cstring>
using namespace std;
long long x[100005],y[100005],z[100005];
int root[1005],cnt1=0,cnt2=0,root1[100005],root2[100005];
int find(int x)//并查集的查找“父亲”
{
if(root[x]==x) return x;
return root[x]=find(root[x]);
}
bool judge(long long x,long long y,long long z,long long x1,long long y1,long long z1,long long r)//判断是否相切或相交
{
if((x-x1)*(x-x1)+(y-y1)*(y-y1)+(z-z1)*(z-z1)<=4*r*r) return true;
else return false;
}
int main()
{
int t;
cin>>t;
for(int i=1;i<=t;i++)
{
long long n,h,r;
int cnt1=0,cnt2=0;
cin>>n>>h>>r;
for(int j=1;j<=n;j++)//并查集初始化
{
root[j]=j;
}
for(int j=1;j<=n;j++)
{
cin>>x[j]>>y[j]>>z[j];
if(z[j]+r>=h)//判断与上表面是否相交或相切
{
cnt1++;root1[cnt1]=j;//进行标记
}
if(z[j]-r<=0)//判断与下表面是否相交或相切
{
cnt2++;root2[cnt2]=j;//进行标记
}
for(int k=1;k<=j;k++)
{
if((x[j]-x[k])*(x[j]-x[k])+(y[j]-y[k])*(y[j]-y[k])>4*r*r) continue;
if(judge(x[j],y[j],z[j],x[k],y[k],z[k],r))//判断并且进行合并
{
int a1=find(j);
int a2=find(k);
if(a1!=a2) root[a2]=a1;
}
}
}
int flag=0;
for(int j=1;j<=cnt1;j++)//列举上表面
{
for(int k=1;k<=cnt2;k++)//列举下表面
{
if(find(root1[j])==find(root2[k]))//判断“父亲”是否为同一个,即是否上表面连通下表面
{
flag=1;
break;
}
}
if(flag) break;
}
if(flag) cout<<"Yes"<<endl;
else cout<<"No"<<endl;
}
return 0;
}
灾后重建
题目背景
B 地区在地震过后,所有村庄都造成了一定的损毁,而这场地震却没对公路造成什么影响。但是在村庄重建好之前,所有与未重建完成的村庄的公路均无法通车。换句话说,只有连接着两个重建完成的村庄的公路才能通车,只能到达重建完成的村庄。
题目描述
给出 B 地区的村庄数
N
N
N,村庄编号从
0
0
0 到
N
−
1
N-1
N−1,和所有
M
M
M 条公路的长度,公路是双向的。并给出第
i
i
i 个村庄重建完成的时间
t
i
t_i
ti,你可以认为是同时开始重建并在第
t
i
t_i
ti 天重建完成,并且在当天即可通车。若
t
i
t_i
ti 为
0
0
0 则说明地震未对此地区造成损坏,一开始就可以通车。之后有
Q
Q
Q 个询问
(
x
,
y
,
t
)
(x,y,t)
(x,y,t),对于每个询问你要回答在第
t
t
t 天,从村庄
x
x
x 到村庄
y
y
y 的最短路径长度为多少。如果无法找到从
x
x
x 村庄到
y
y
y 村庄的路径,经过若干个已重建完成的村庄,或者村庄
x
x
x 或村庄
y
y
y 在第
t
t
t 天仍未重建完成,则需要返回 -1
。
输入格式
第一行包含两个正整数 N , M N,M N,M,表示了村庄的数目与公路的数量。
第二行包含 N N N个非负整数 t 0 , t 1 , … , t N − 1 t_0, t_1,…, t_{N-1} t0,t1,…,tN−1,表示了每个村庄重建完成的时间,数据保证了 t 0 ≤ t 1 ≤ … ≤ t N − 1 t_0 ≤ t_1 ≤ … ≤ t_{N-1} t0≤t1≤…≤tN−1。
接下来 M M M行,每行 3 3 3个非负整数 i , j , w i, j, w i,j,w, w w w为不超过 10000 10000 10000的正整数,表示了有一条连接村庄 i i i与村庄 j j j的道路,长度为 w w w,保证 i ≠ j i≠j i=j,且对于任意一对村庄只会存在一条道路。
接下来一行也就是 M + 3 M+3 M+3行包含一个正整数 Q Q Q,表示 Q Q Q个询问。
接下来 Q Q Q行,每行 3 3 3个非负整数 x , y , t x, y, t x,y,t,询问在第 t t t天,从村庄 x x x到村庄 y y y的最短路径长度为多少,数据保证了 t t t是不下降的。
输出格式
共 Q Q Q行,对每一个询问 ( x , y , t ) (x, y, t) (x,y,t)输出对应的答案,即在第 t t t天,从村庄 x x x到村庄 y y y的最短路径长度为多少。如果在第t天无法找到从 x x x村庄到 y y y村庄的路径,经过若干个已重建完成的村庄,或者村庄x或村庄 y y y在第 t t t天仍未修复完成,则输出 − 1 -1 −1。
样例 #1
样例输入 #1
4 5
1 2 3 4
0 2 1
2 3 1
3 1 2
2 1 4
0 3 5
4
2 0 2
0 1 2
0 1 3
0 1 4
样例输出 #1
-1
-1
5
4
提示
对于 30 % 30\% 30%的数据,有 N ≤ 50 N≤50 N≤50;
对于 30 % 30\% 30%的数据,有 t i = 0 t_i= 0 ti=0,其中有 20 % 20\% 20%的数据有 t i = 0 t_i = 0 ti=0且 N > 50 N>50 N>50;
对于 50 % 50\% 50%的数据,有 Q ≤ 100 Q≤100 Q≤100;
对于 100 % 100\% 100%的数据,有 N ≤ 200 N≤200 N≤200, M ≤ N × ( N − 1 ) / 2 M≤N \times (N-1)/2 M≤N×(N−1)/2, Q ≤ 50000 Q≤50000 Q≤50000,所有输入数据涉及整数均不超过 100000 100000 100000。
解题方法
我们观察本题,是查找最短路径,所以我们会想到使用floyd。
在本题要求下,只需要在每次规定的时间前把修好的路进行更新,即进行floyd求最短路。
解题代码
#include<iostream>
#include<cstring>
using namespace std;
int n,m,q,now=0;
int tt[205],f[205][205];
void floyd(int k)
{
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
{
f[i][j]=min(f[i][j],f[i][k]+f[k][j]);//Floyd模板
}
}
return;
}
int main()
{
cin>>n>>m;
for(int i=0;i<n;i++)
{
cin>>tt[i];
}
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
{
if(i==j) f[i][i]=0;//每个村庄原地不动初始为0
else f[i][j]=0x3ffffff;//先把每个村庄之间的路径设为最大值
}
}
int x,y,d;
for(int i=0;i<m;i++)
{
cin>>x>>y>>d;
f[x][y]=f[y][x]=d;
}
cin>>q;
int a,b,t;
for(int i=0;i<q;i++)
{
cin>>a>>b>>t;
while(tt[now]<=t&&now<n)//在时间要求之间更新路径
{
floyd(now);
++now;
}
if(tt[a]>t||tt[b]>t) cout<<-1<<endl;//如果要求的时间内无法到达,则输出-1
else
{
if(f[a][b]==0x3ffffff) cout<<-1<<endl;//无法通路也输出-1
else cout<<f[a][b]<<endl;
}
}
return 0;
}
[HAOI2006]聪明的猴子
题目描述
在一个热带雨林中生存着一群猴子,它们以树上的果子为生。昨天下了一场大雨,现在雨过天晴,但整个雨林的地表还是被大水淹没着,部分植物的树冠露在水面上。猴子不会游泳,但跳跃能力比较强,它们仍然可以在露出水面的不同树冠上来回穿梭,以找到喜欢吃的果实。
现在,在这个地区露出水面的有N棵树,假设每棵树本身的直径都很小,可以忽略不计。我们在这块区域上建立直角坐标系,则每一棵树的位置由其所对应的坐标表示(任意两棵树的坐标都不相同)。
在这个地区住着的猴子有M个,下雨时,它们都躲到了茂密高大的树冠中,没有被大水冲走。由于各个猴子的年龄不同、身体素质不同,它们跳跃的能力不同。有的猴子跳跃的距离比较远(当然也可以跳到较近的树上),而有些猴子跳跃的距离就比较近。这些猴子非常聪明,它们通过目测就可以准确地判断出自己能否跳到对面的树上。
【问题】现已知猴子的数量及每一个猴子的最大跳跃距离,还知道露出水面的每一棵树的坐标,你的任务是统计有多少个猴子可以在这个地区露出水面的所有树冠上觅食。
输入格式
输入文件monkey.in包括:
第1行为一个整数,表示猴子的个数M(2<=M<=500);
第2行为M个整数,依次表示猴子的最大跳跃距离(每个整数值在1–1000之间);
第3行为一个整数表示树的总棵数N(2<=N<=1000);
第4行至第N+3行为N棵树的坐标(横纵坐标均为整数,范围为:-1000–1000)。
(同一行的整数间用空格分开)
输出格式
输出文件monkey.out包括一个整数,表示可以在这个地区的所有树冠上觅食的猴子数。
样例 #1
样例输入 #1
4
1 2 3 4
6
0 0
1 0
1 2
-1 -1
-2 0
2 2
样例输出 #1
3
提示
【数据规模】
对于40%的数据,保证有2<=N <=100,1<=M<=100
对于全部的数据,保证有2<=N <= 1000,1<=M=500
感谢@charlie003 修正数据
解题方法
此题的解题方法为生成树,相当于并查集和最短路合成得到的“提高版”。
我们只需要把每条边进行Kruskal算法,然后加入到同一个连通图,每次更新最大边。
Kruskal算法就是每次选所有边的“最短”的边,加入同一个连通图,之后同理寻找,但是此题不用在意是否形成回路,只需要将每个树加入同一个连通图并更新最大边,所以类似最小生成树
解题代码
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
int n,m,ans=0;
double num=-1.0;
int a[1000005][3],ty[1000005],root[1000005];
struct tre{
int t1,t2;//标记树的边端点
double dis;//记录边大小
}tree[1000005];
int cmp(tre a,tre b)
{
return a.dis<b.dis;
}
int find(int x)//并查集查找
{
if(root[x]==x) return x;
return root[x]=find(root[x]);
}
int main()
{
cin>>m;
for(int i=1;i<=m;i++)
{
cin>>ty[i];
}
cin>>n;
for(int i=1;i<=n;i++)
{
cin>>a[i][0]>>a[i][1];
}
int k=0;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
if(i!=j)
{
k++;
tree[k].t1=i;//记录树的两端点
tree[k].t2=j;
tree[k].dis=sqrt((a[i][0]-a[j][0])*(a[i][0]-a[j][0])+(a[i][1]-a[j][1])*(a[i][1]-a[j][1]));//计算边的大小
}
}
}
int pos=n;
sort(tree+1,tree+k+1,cmp);//进行Kruskal算法要排序
for(int i=1;i<=n;i++)
{
root[i]=i;//并查集
}
for(int i=1;i<=k;i++)//并查集
{
if(pos==1) break;//判断每棵数是否加入同一个连通图
int s1=find(tree[i].t1);
int s2=find(tree[i].t2);
if(s1!=s2)
{
root[s1]=s2;
pos--;
num=tree[i].dis;
}
}
for(int i=1;i<=m;i++)
{
if(num<=ty[i]) ans++;
}
cout<<ans<<endl;
return 0;
}