HDU-4803-Poor Warehouse Keeper（贪心+二分）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4803

题目大意：给出两个数x，y。有两种操作，

1.将x+1，同时y+y/x;

2.x不变，y+1；

问最少操作多少次，能够将1，1，变为x，y。

如果不能够变成x，y形式，输出-1。

思路：显然，如果x>y的话，必定不可能。然后我们只有考虑1，1如何变成x，y即可。

经过简单的样例推导，可以发现，如果想要尽可能少的操作。那么使得每个a对应的b尽可能的变大。

那么我们的工作就变成了对于从1~x中的每一个数，从1开始，都找到他们对应的最大的b即可。

一开始我是每次b+1，很明显超时，比赛的时候心态有点炸，想到了用二分，但是总是写挂（还是不熟悉）。

赛后静下心来多看了看，就轻易的得出了结论。还是心态问题啊，其实还剩半个小时的，完全有时间重新写一个二分出来。

看起来挺简单，其实坑还不少。注意每次+1是整数，double和int之间的转换。

ACCode：

//#pragma comment(linker, "/STACK:1024000000,1024000000")
  
#include<stdio.h>
#include<string.h> 
#include<math.h> 
   
#include<map>  
#include<set>
#include<deque> 
#include<queue> 
#include<stack> 
#include<bitset>
#include<string> 
#include<fstream>
#include<iostream> 
#include<algorithm> 
using namespace std; 
  
#define ll long long 
#define Pair pair<ll,ll>
//#define max(a,b) (a)>(b)?(a):(b)
//#define min(a,b) (a)<(b)?(a):(b)
#define clean(a,b) memset(a,b,sizeof(a))// ??
//std::ios::sync_with_stdio(false);
//  register
const int MAXN=2500+10;
const int INF32=0x3f3f3f3f;
const ll INF64=0x3f3f3f3f3f3f3f3f;
const ll mod=1e9+7;
const double PI=acos(-1.0);
const double EPS=1.0e-8;


double x,y;
int cnt;

void DFS(double a,double b){//计算a 时 b的最大值 
//	printf("a=%lf b=%lf cnt=%d\n",a,b,cnt);
	double la=x-a,nb=b,price=b/a;
	int l=0,r=(int)y-(int)b,mid;
	while(l<=r){//最多 加多少 
		mid=(l+r)>>1;
		price=(b+mid*1.0)/a;nb=b+mid*1.0;
		if(la*price+nb<y+1.0) l=mid+1;
		else r=mid-1;
	}nb=r*1.0+b;price=nb/a;
	cnt+=r;//得出b的最大值 
	if((x-a)<=EPS&&(a-x)>=-EPS) return ;
	++cnt;//加成，cnt+1 
	DFS(a+1,price*(a+1));
}
int main(){
	while(~scanf("%lf%lf",&x,&y)){
		if(x-y>EPS){
			printf("-1\n");continue;
		}cnt=0;
		DFS(1,1);
		printf("%d\n",cnt);
	}
}


/*

Sample Input
1 1
3 8
9 31
Sample Output
0
5
11

*/