在金融中,我们有时会用内部收益率IRR来评价项目的投资财务效益,它等于使得投资净现值NPV等于0的贴现率。换句话说,给定项目的期数T、初始现金流CF0和项目各期的现金流CF1, CF2, …,CFT,IRR是下面方程的解:
为了简单起见,本题假定:除了项目启动时有一笔投入(即初始现金流CF0 < 0)之外,其余各期均能赚钱(即对于所有i=1,2,…,T,CFi > 0)。根据定义,IRR可以是负数,但不能大于-1。
Input
输入文件最多包含25组测试数据,每个数据占两行,第一行包含一个正整数T(1<=T<=10),表示项目的期数。第二行包含T+1个整数:CF0, CF1, CF2, …, CFT,其中CF0 < 0, 0 < CFi < 10000 (i=1,2,…,T)。T=0表示输入结束,你的程序不应当处理这一行。
Output
对于每组数据,输出仅一行,即项目的IRR,四舍五入保留小数点后两位。如果IRR不存在,输出”No”,如果有多个不同IRR满足条件,输出”Too many”(均不含引号)
Sample Input
1
-1 2
2
-8 6 9
0
Sample Output
1.00
0.50
分析:给了一个IRR为未知数的方程,求方程的解,这是一个一元高次方程,它的解一定存在,而题目中有一处错误,, 应该是错误吧, IRR的范围应该是不小于-1吧? 看样例就不是大于。。
然后这里IRR的具体范围不知道,就设定一个最小值l=-1.0,最大值r尽可能的大一些,我用的1e8,小了的话可能WA吧。 然后二分查找
#include
#include
#include
#include
#include
using namespace std;
#define mem(a,n) memset(a,n,sizeof(a))
int r1,r2,r3,c1,c2,c3;
int dxy[8][2]= {-1,0,0,-1,1,0,0,1,-1,-1,-1,1,1,-1,1,1};
const int N=15;
const int eps=1e-6;
int main()
{
int T;
double a[N];
while(~scanf("%d",&T)&&T)
{
for(int i=0; i<=T; i++)
scanf("%lf",&a[i]);
double l=-1.0,r=1e8;
double mid,sum=0,t;
for(int i=0;i<50;i++)
{
sum=0,t=1.0;
mid=(l+r)/2;
for(int i=1; i<=T; i++)
{
t/=(1.0+mid);
sum+=t*a[i];
}
if(sum+a[0]>0)
l=mid;
else r=mid;
}
printf("%.2lf\n",mid);
}
return 0;
}