思路:迭代加深搜索。首先先枚举由多少个数构成。然后把当前剩下的分子和分母等传进去。如果发现最优解,更新值。
代码
#include<bits/stdc++.h>
using namespace std;
#define int long long
int a,b,dm,ans[100001],f[100001];
bool flag;
int gcd(int a,int b){
int t;
while(t=a%b){
a=b;
b=t;
}
return b;
}
void dfs(int a,int b,int x,int last){
if(x==dm){
if(a!=1)
return;
ans[x]=b;
if(ans[x]<=ans[x-1])
return;
if(!flag||f[x]>ans[x]){
for(int i=1;i<=x;i++)
f[i]=ans[i];
}
flag=1;
return;
}
for(int i=max(last+1,b/a);i<b*(dm-x+1)/a;++i){
ans[x]=i;
int g=gcd(a*i-b,b*i);
int aa=(a*i-b)/g;
int bb=b*i/g;
dfs(aa,bb,x+1,i);
}
}
signed main(){
cin>>a>>b;
ans[0]=1;
flag=0;
for(dm=1;;dm++){
dfs(a,b,1,b/a);
if(flag)break;
}
for(int i=1;i<=dm;i++)
cout<<f[i]<<' ';
return 0;
}