pku 1980 DFS

/*
题目链接http://poj.org/problem?id=1980
题意描述：输入p,q,a,n 给定一个分数p/q. 求满足 p/q = 1/a1 + 1/a2 + 1/a3 +...+ 1/ai 所有的组合情况。并且有2个限制条件，i<=n,a1*a2*...*ai <= a
解题思路：由于n<=7.并且a1*a2*...*ai<12000。 可以用暴力搜索，搜索过程有一个很关键的剪枝，就是判断剩余的分数是否大于当前分母所能达到的最大值，若大于，则可以跳出。具体见

代码
*/

/* 题目链接http://poj.org/problem?id=1980 题意描述：输入p,q,a,n 给定一个分数p/q. 求满足 p/q = 1/a1 + 1/a2 + 1/a3 +...+ 1/ai 所有的组合情况。并且有2个限制条件，i<=n,a1*a2*...*ai <= a 解题思路：由于n<=7.并且a1*a2*...*ai<12000。 可以用暴力搜索，搜索过程有一个很关键的剪枝，就是判断剩余的分数是否大于当前分母所能达到的最大值，若大于，则可以跳出。具体见 代码 */ #include <iostream> using namespace std; int p,q,a,n; int ans; int gcd(int a, int b) //最大公约数 { return b ? gcd(b,a%b) : a; } int lcm(int a, int b) //最小公倍数 { return a*b/gcd(a,b); } //上一层的分母k， 分数和的分子son，分数和的分母mother， 分母积product， 由多少个分数组成cnt void dfs(int k, int son, int mother, int product, int cnt) { //超出n个数或者分数大于给定的分数值 if(cnt == n+1 || son*q > p*mother) { return ; } //满足条件 if(son == p && mother == q) { ans++; return ; } for(int i=k; i<=a/product; i++) { //关键剪枝，预判断。若剩余的分数大于当前分母所能达到的最大值 if((i*son+(n-cnt)*mother)*q < p*i*mother) // son/mother + (n-cnt) / i < p/q { return ; } int flag = lcm(mother,i); int tempSon,tempMother,temp; //分数求和... tempMother = flag; tempSon = flag*son/mother + flag / i; temp = gcd(tempMother,tempSon); tempMother = tempMother /temp; tempSon = tempSon / temp; //往下搜 dfs(i,tempSon,tempMother,product*i,cnt+1); } } int main() { while(cin>>p>>q>>a>>n,p+q+a+n) { ans = 0; int flag = gcd(p,q); //化简给定的分数 p = p / flag; q = q / flag; dfs(1,0,1,1,0); cout<<ans<<endl; } return 0; }