算24

描述

给出4个小于10个正整数，你可以使用加减乘除4种运算以及括号把这4个数连接起来得到一个表达式。现在的问题是，是否存在一种方式使得得到的表达式的结果等于24。

这里加减乘除以及括号的运算结果和运算的优先级跟我们平常的定义一致（这里的除法定义是实数除法）。

比如，对于5，5，5，1，我们知道5 * (5 – 1 / 5) = 24，因此可以得到24。又比如，对于1，1，4，2，我们怎么都不能得到24。

输入

输入数据包括多行，每行给出一组测试数据，包括4个小于10个正整数。最后一组测试数据中包括4个0，表示输入的结束，这组数据不用处理。

输出

对于每一组测试数据，输出一行，如果可以得到24，输出“YES”；否则，输出“NO”。

代码：

#include <iostream>  
#include <cmath>  
using namespace std;  
const double ZERO = 1e-6; 
const int NUM = 4;   
bool flag = false;  
double f[NUM];  


bool Function(int step)
{
if(step == 1)  
    {  
        if(fabs(f[0] - 24) < ZERO)   //注释*
        {  
            flag = true;  
            return true;  
        }  
        else  
            return false;  
    }  
    for(int i = 0; i < step; ++i)  
    {  
        for(int j = i + 1; j < step; ++j)  
        {  
            double tmp_a = f[i], tmp_b = f[j];  
            f[j] = f[step - 1];  
            f[i] = tmp_a + tmp_b;  
            Function(step - 1);  
            f[i] = tmp_a - tmp_b;  
            Function(step - 1);  
            f[i] = tmp_b - tmp_a;  
            Function(step - 1);  
            f[i] = tmp_a * tmp_b;  
            Function(step - 1);  
            if(fabs(tmp_b - 0.0) > ZERO)  
            {  
                f[i] = tmp_a / tmp_b;  
                Function(step - 1);  
            }  
            if(fabs(tmp_a - 0.0) > ZERO)  
            {  
                f[i] = tmp_b / tmp_a;  
                Function(step - 1);  
            }  
            f[i] = tmp_a;  
            f[j] = tmp_b;  
        }  
    }  
    return true;  
}


int main()  
{  
    while(cin >> f[0]>>f[1]>>f[2]>>f[3])  
    {  
        if(f[0] == 0)  
            break;   
        flag = false;  
        Function(NUM);  
        if(flag)  
            cout << "YES" << endl;  
        else  
            cout << "NO" << endl;  
    }  
    return 0;  
}

* float，double分别遵循R32-24,R64-53的标准，

他们的位数分别是23,52；即误差在2^-23,2^-52；

所以float的精度误差在1e-6；double精度误差在1e-15 
所以要判断一个单精度浮点数：则是if( fabs(f) <= 1e-6)；
要判断一个双精度浮点数：则是if( fabs(f) <= 1e-15 )；