标签：
|
description |
|
对于每一个非负有理数，我们知道它一定能划归成某些特殊真分数之和，特殊真分数要满足它们的分子为1，但是我们知道，对于无穷级数1/2+1/3+1/4…。虽然，它是发散的，但是改级数增长得极为缓慢，例如到了数百万之后，和也在18~19左右。 若干年来，不断有人宣称发现了该级数的特殊性质，这些都对这个问题的研究起到了深远的影响。 你的任务来了，要求给你个真分数，你需要将其化简为最少的若干特殊真分数之和，你要输出这个序列（序列按递增序）。 如果有不同的方案，则分数个数相同的情况下使最大的分母最小。若还相同，则使次大的分母最大……以此类推。 如：2/3=1/2+1/6，但不允许2/3=1/3+1/3，因为加数中有相同的。 对于一个分数a/b，表示方法有很多种，但是哪种最好呢？ 首先，加数少的比加数多的好，其次，加数个数相同的，最小的分数越大越好。如： 19/45=1/3 + 1/12 + 1/180 19/45=1/3 + 1/15 + 1/45 19/45=1/3 + 1/18 + 1/30 19/45=1/4 + 1/6 + 1/180 19/45=1/5 + 1/6 + 1/18 最好的是最后一种，因为18 比180, 45, 30,都小。 |
对于此类搜索问题，搜索深度没有明显的上界，而且加数的选择在理论上也是无限的，也就是说宽度搜索连一层都扩展不完。
利用迭代加深搜索（IDA*）可以解决上面的问题。一方面我们要解决，利用BFS从小到大枚举深度上限maxd，虽然理论上深度是无限的，但是只要保证有解，深度值必然在有限的时间内能枚举到。
另一方面，我们还要解决每次枚举的值无限的问题，可以借助maxd来“剪枝”，以此题为例，每一次枚举的的分母个数必然有一个结束值，否则就会出现无限下去的尴尬了。那么当扩展到第i层时，第i层分数值为1/e，那么接下来由于分母是以递增顺序进行的，那么分数值都不会大于1/e，如果c/d（前i个分数之和）+1/e*（maxd-i）<a/b，可以直接break;因为起码深度不够，要再加深度。
基本框架：
for（maxd=1；；maxd++） { if(dfs(0,,,)) { ok=1; break; } }
代码如下：
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #define MAXN 20000 //一定要用64位的整数，由于里面存在化归为同分母的操作，需要相乘 using namespace std; int64_t maxd; int64_t v[MAXN],ans[MAXN];//v是暂时存放的满足题意的分母的数组，ans是记录满足题意的最优分母值的数组。 int64_t gcd(int64_t a,int64_t b)//求最大公约数 { int64_t m; while(a!=0) { m=b%a; b=a; a=m; } return b; } int64_t getfirst(int64_t a,int64_t b)//取比a/b小的最大分数，分子必须为1 { int64_t i,j; for(i=2;;i++) { if(b<a*i) { break; } } return i; } bool better(int64_t d)//必要的判断，这组分数之和是否满足最优解 { int64_t i; bool flag=true; if(ans[d]==-1)//由于初始化ans是-1，如果是第一次出现的满足题意的解，返回true return true; for(i=d;i>=0;i--) { if(v[i]==ans[i])//从高位进行判断大小，题意要求 continue; else if(v[i]>ans[i]) { //return false; flag=false; break; } else { break; } } if(flag==true) return true; else return false; } bool dfs(int64_t a,int64_t b,int64_t from,int64_t d)//深度为d { int64_t aa,bb,g,i; bool ok; // cout<<from<<endl; if(d==maxd) { if(a!=1) return false; for(i=0;i<=d-1;i++) { if(v[i]==b) { return false; } } v[d]=b; sort(v,v+d); if(better(d)) memcpy(ans,v,sizeof(int64_t)*(d+1)); return true; } ok=false; //重要！！！ from=max(from,getfirst(a,b));//枚举起点，去上一次加一的分母值和比a/b小的最大分数的分母中更大的。 for(i=from;;i++) { //剪枝，如果c/d（前i个分数之和）+1/e*（maxd-i）<a/b，可以直接break; if(b*(maxd+1-d)<=a*i) break; v[d]=i; aa=a*i-b; bb=b*i; g=gcd(aa,bb); aa=aa/g; bb=bb/g;//约分 if(dfs(aa,bb,i+1,d+1)) ok=true; } return ok; } int main() { int64_t a,b,aa,bb,i,g; while(cin>>a>>b) { if(a==0) { cout<<a<<"/"<<b<<"=0"<<endl; continue; } memset(ans,-1,sizeof(ans)); g=gcd(a,b); aa=a/g; bb=b/g; if(aa==1) { printf("%d/%d=%d/%d\n",a,b,aa,bb); } else { for(maxd=1;;maxd++) { if(dfs(aa,bb,getfirst(aa,bb),0)) break; } cout<<a<<"/"<<b<<"="; for(i=0;i<=maxd-1;i++) cout<<"1/"<<ans[i]<<"+"; cout<<"1/"<<ans[i]; cout<<endl; } } return 0; }
标签：
原文地址：http://www.cnblogs.com/fancy-itlife/p/4297269.html
