P1763 埃及分数

最新推荐文章于 2024-08-23 15:54:50 发布

NNNNwang

最新推荐文章于 2024-08-23 15:54:50 发布

阅读量165

点赞数

分类专栏：搜索文章标签：深度优先 java 算法

本文链接：https://blog.csdn.net/qq_62525850/article/details/128488981

版权

搜索专栏收录该内容

2 篇文章 0 订阅

订阅专栏

该博客文章展示了一个基于Java的程序，用于解决特定的数学问题。代码通过递归和动态规划寻找一组分数，使得分母尽可能小。算法涉及高斯欧几里得算法计算最大公约数，并进行深度优先搜索。程序最终输出一系列分数及其对应的最优解。

摘要由CSDN通过智能技术生成

package newPro;
import java.util.*;
public class pro11 {
	static final int N=100000;
	static int a,b;
	static long ans[]=new long[N];
	static long v[]=new long[N];
	static int maxd;
	
	public static void main(String args[])
	{ 
		Scanner in=new Scanner(System.in);
		a=in.nextInt();
		b=in.nextInt();
		for( maxd=1;;maxd++)
		{
			Arrays.fill(ans, -1);
			if(dfs(0,get_first(a,b),a,b)) break; 
			//System.out.println(maxd);
		}
			
		for(int i=0;i<=maxd;i++)
		{
			System.out.print(ans[i]+" ");
		}
		
	}
	
	public static long get_first(long aa,long bb)
	{
		return bb/aa+1;//
	}
	public static boolean better(int d)
	{
		for(int i=d;i>=0;i--)
		{
			if(v[i]!=ans[i]) {
				return ans[i]==-1||v[i]<ans[i];//最小的分数要越大，即分母要越小好
			}
		}
		return false;
	}
	
	public static long gcd(long x,long y)
	{
		return y==0? x:gcd(y,x%y);
	}
	public static boolean dfs(int d,long from,long aa,long bb)
	{
		if(d==maxd)
		{
			//判断是否为最约束
			if(aa!=1) return false;
			v[d]=bb;//加入最后的分数，保存的为加入分数的倒数，即原来的分母
			if(better(d)) {ans=Arrays.copyOf(v,d+1);}
			return true;
		}
		boolean ok=false;
		from=Math.max(from, get_first(aa,bb));//上一个循环直接记录的是i+1,因此在当前循环要再次判断from
		for(long i=from;;i++)
		{
			if(bb*(maxd-d+1)<=aa*i) break;
			v[d]=i;
			long a2,b2;
			a2=aa*i-bb;
			b2=bb*i;
			long g=gcd(a2,b2);
			if(dfs(d+1,i+1,a2/g,b2/g)) ok=true;//使用的是i+1
		}
		return ok;
	}
}