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));
}