内部收益率（二分法）_二分法求到期收益率

在金融中，我们有时会用内部收益率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<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
#include<algorithm>
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;
}
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38