poj 1980 Unit Fraction Partition…

这个用了double来判断相等与否，精度不太好把握，而且需要注意的东西也多一些

      搜索顺下，t从小到大

剪枝  tmp>sum  tmp为接下来所能选择的最大的分母，sum为剩下的值

      (1/t)*n<sum 不再搜索 剩下的值都选最大仍不能满足

       当乘积大于a时 退出

 

#include <iostream>
#include<stdio.h>
#include<string.h>
#include<math.h>
using namespace std;
int ans;
double e=0.000000000001;
int dfs(double sum,int t,int n,double a)
{
//    printf("%f %d %d %f\n",sum,t,n,a);

    if(fabs(sum)<e)
    {
        if(a>=(1-e)) ans++;
        return(0);
    }

    double tmp=(double)1/a;
    if(tmp-e>sum) return(0);

    for(int i=t;;i++)
    {
        double tmp=(double)1/(double)i;
        if(tmp-e>sum) continue;
        if(tmp*n+e<sum) break;
        if(i-e>a) break;
        dfs(sum-tmp,i,n-1,a/i);
    }
}

int main()
{
    int p,q,a,n;
    while(1)
    {
        scanf("%d %d %d %d",&p,&q,&a,&n);
        if(p==0&&q==0&&a==0&&n==0) break;
        double sum=(double)p/(double)q;
        ans=0;
        dfs(sum,1,n,(double)a);
        printf("%d\n",ans);
    }
    return 0;
}