BIT1060 Is It A Tree?

题意

树的数据结构定义为

根节点没有一条边指向它 

除根节点外每个节点必有刚好一条边指向它 

从根节点刚好只有一条路径可以到达每个点

样例有多组,每组样例以0 0结尾

如果样例是-1 -1,程序终止

每组样例的输入是n对数,一对数a和b代表a到b有一条有向边

解法:

原来想怎么解来着忘了,反正弄了一堆很麻烦的玩意还RE了

后来用

入度为0的点只有一个,其余都是入度为1的点这个性质过了

#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;
int counter[1010];//每个点的入度
bool exist[1010];//在此组样例中出现的点
int main()
{
	int a,b;
	int kk=1;
	while(scanf("%d %d",&a,&b))
	{
		if(a==0&&b==0)
		{
			printf("Case %d is a tree.\n",kk++);
		}
		if(a==-1&&b==-1)
		{
			break;
		}
		memset(counter,0,sizeof(counter));
		memset(exist,false,sizeof(exist));
		counter[b]++;
		exist[a]=true;
		exist[b]=true;
		while(scanf("%d %d",&a,&b),a||b)
		{
			counter[b]++;
			exist[a]=true;
			exist[b]=true;
		}
		bool ans=true;
		int num=0;//出现过的且入度为0的点
		for(int i=1;i<=1000;i++)//所有点的入度应为0或1
		{
			if(counter[i]>1)
			{
				ans=false;
			}
			if(counter[i]==0&&exist[i])
			{
				num++;
			}
		}
		if(num!=1)
		{
			ans=false;
		}
		if(ans)
		{
			printf("Case %d is a tree.\n",kk++);
		}
		else
		{
			printf("Case %d is not a tree.\n",kk++);
		}
	}
	return 0;
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
R R version 4.2.2 (2022-10-31) -- "Innocent and Trusting" Copyright (C) 2022 The R Foundation for Statistical Computing Platform: x86_64-conda-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors.Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. library(ape) setwd("/ifs1/User/dengwei/NTF_data/7.14/rooted_species_tree") species_tree <- read.tree("species_tree.treefile")> compare_trees <- function(gene_tree_file, species_tree) { gene_tree <- read.tree(gene_tree_file) diff_count <- comparePhylo(gene_tree, species_tree, force.rooted = TRUE) return(diff_count) } batch_compare_trees <- function(gene_tree_folder, species_tree) { gene_tree_files <- list.files(path = gene_tree_folder, pattern = ".treefile", full.names = TRUE) diff_counts <- data.frame(Gene_Tree_File = gene_tree_files, Diff_Count = numeric(length(gene_tree_files)), stringsAsFactors = FALSE) for (i in seq_along(gene_tree_files)) { gene_tree_file <- gene_tree_files[i] diff_counts$Diff_Count[i] <- compare_trees(gene_tree_file, species_tree) } return(diff_counts) } gene_tree_folder <- "/ifs1/User/dengwei/NTF_data/7.14/rooted_gene_tree" diff_counts <- batch_compare_trees(gene_tree_folder, species_tree) Error in if (n1 == n2) paste("Both trees have the same number of tips:", : the condition has length > 1
07-15
这个错误可能是由于某个基因树文件与物种树的物种数量不匹配引起的。请确保所有的基因树文件都具有相同的物种数量,并且与物种树的物种数量一致。你可以检查一下基因树文件和物种树文件中的物种数量是否匹配。以下是修改后的代码,包括错误处理,以帮助你找出问题所在: ```R library(ape) # 设置文件路径 species_tree_file <- "/ifs1/User/dengwei/NTF_data/7.14/rooted_species_tree/species_tree.treefile" gene_tree_folder <- "/ifs1/User/dengwei/NTF_data/7.14/rooted_gene_tree" # 读取物种树 species_tree <- read.tree(species_tree_file) # 定义函数:比较基因树和物种树的差异数量 compare_trees <- function(gene_tree_file, species_tree) { gene_tree <- read.tree(gene_tree_file) diff_count <- comparePhylo(gene_tree, species_tree, force.rooted = TRUE) return(diff_count) } # 定义函数:批量比较基因树和物种树 batch_compare_trees <- function(gene_tree_folder, species_tree) { gene_tree_files <- list.files(path = gene_tree_folder, pattern = ".treefile", full.names = TRUE) diff_counts <- data.frame(Gene_Tree_File = gene_tree_files, Diff_Count = numeric(length(gene_tree_files)), stringsAsFactors = FALSE) for (i in seq_along(gene_tree_files)) { gene_tree_file <- gene_tree_files[i] diff_counts$Diff_Count[i] <- tryCatch({ compare_trees(gene_tree_file, species_tree) }, error = function(e) { NA }) } return(diff_counts) } # 执行批量比较基因树和物种树 diff_counts <- batch_compare_trees(gene_tree_folder, species_tree) ``` 在这个修改后的代码中,我添加了一个`tryCatch`块来处理可能出现的错误。如果在比较基因树和物种树时出现错误,将会返回NA作为差异数量。这可以帮助你找出具体哪个基因树文件导致了错误。请确保在运行代码之前,已经正确设置了`species_tree_file`和`gene_tree_folder`变量的值。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值