思路:使用二分法求方程的近似解。
#include<iostream>
#include<cmath>
#include<iomanip>
using namespace std;
int main(){
double IRR,NPV,l,r;
int T;
while(cin>>T&&T!=0){
double *p=new double[T+1];
for(int i=0;i<T+1;i++){
cin>>p[i];//输入系数
}
NPV=0.0;l=1E6;r=-1+1e-6;
for(int i=0;i<1000;i++){
IRR=(l+r)/2.0;
NPV=0.0;
for(int j=0;j<T+1;j++){
NPV+=1.0*p[j]/pow((1+IRR),j);
}
if(fabs(NPV)<1e-9) break;//达成要求
if(NPV>0) r=IRR;
else l=IRR;
}
cout<<std::fixed<<setprecision(2)<<IRR<<endl;
}
return 0;
}
一开始只有例题答案正确,提交上去总是错误百分之50,后来查找资料对比别人的代码,发现自己的上界设的太小,改成1e6后答案正确。对自己来说也是学到了,还要弥补的东西太多了。