递归算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”。
样例输入
5 5 5 1
1 1 4 2
0 0 0 0
样例输出
YES

NO

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int DIS=1e-6;
double a[4];
int iszero(double x)
{
    if(fabs(x)<=DIS)
    return 1;
    else return 0;
}
int cal(double a[],int n)
{
    if(n==1)
    {
        if(iszero(a[0]-24))
            return 1;
        else
            return 0;
    }
    double b[5];
    memset(b,0,sizeof(b));
    for(int i=0;i<n-1;i++)
    {
        for(int j=i+1;j<n;j++)//任取两个数
    {
        int m=0;
        for(int k=0;k<n;k++)
        {
            if(k!=j&&k!=i)
            {
                b[m++]=a[k];
            }
        }
        b[m]=a[i]+a[j];
        if(cal(b,m+1))return 1;
        b[m]=a[i]-a[j];
        if(cal(b,m+1))return 1;
        b[m]=a[j]-a[i];
        if(cal(b,m+1))return 1;
        b[m]=a[i]*a[j];
        if(cal(b,m+1))return 1;
        if(!iszero(a[i]))
        {
        b[m]=a[j]/a[i];
        if(cal(b,m+1))return 1;
        }
        if(!iszero(a[j]))
        {
        b[m]=a[i]/a[j];
        if(cal(b,m+1))return 1;
        }
    }
    }
    return 0;
}
int main()
{
    while(1)
    {
   // memset(a,0,sizeof(a));
    for(int i=0;i<4;i++)
        cin>>a[i];
    if(a[0]==0&&a[1]==0&&a[2]==0&&a[3]==0)break;
    if(cal(a,4))cout<<"YES"<<endl;
    else cout<<"NO"<<endl;
    }
    return 0;
}