PAT 1063 计算谱半径

2019-04-14 11:55发布

1063 计算谱半径(20 分)

在数学中,矩阵的“谱半径”是指其特征值的模集合的上确界。换言之,对于给定的 n 个复数空间的特征值 { a1" role="presentation">a1​ +b1i" role="presentation">b1i,⋯,an+bni" role="presentation">an+bni },它们的模为实部与虚部的平方和的开方,而“谱半径”就是最大模。
现在给定一些复数空间的特征值,请你计算并输出这些特征值的谱半径。
输入格式:
输入第一行给出正整数 N(10000" role="presentation">10000)是输入的特征值的个数。随后 N 行,每行给出 1 个特征值的实部和虚部,其间以空格分隔。注意:题目保证实部和虚部均为绝对值不超过 1000 的整数。
输出格式:
在一行中输出谱半径,四舍五入保留小数点后 2 位。
输入样例: 5 0 1 2 0 -1 0 3 3 0 -3 输出样例: 4.24


解析

浮点比较:
这里写图片描述
完整的可以看算法笔记
小数点两位四舍五入的函数 double round(double max){ int temp = max*1000; int temp1 = temp%10; if(temp1<5){ return (double)(temp-temp1)/1000; } else{ return (double)(temp-temp1+10)/1000; } } code: #include using namespace std; const double eps = 1e-8; #define Less(a,b) (((a)-(b))<(-eps)) double round(double max){ int temp = max*1000; int temp1 = temp%10; if(temp1<5){ return (double)(temp-temp1)/1000; } else{ return (double)(temp-temp1+10)/1000; } } int main() { int N,a,b; double max = 0,temp; cin>>N; while(N--){ cin>>a>>b; temp = sqrt(a*a+b*b); if(Less(max,temp)) max=temp; } printf("%.2lf",round(max)); }