内部收益率(二分法)

在金融中,我们有时会用内部收益率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;
}

你可能感兴趣的:(ACM_数学,湖南省赛,ACM_二分,三分)