OJ算法实验题 内部收益率

思路:使用二分法求方程的近似解。

#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后答案正确。对自己来说也是学到了，还要弥补的东西太多了。