UVA1616 Caravan Robbers（二分答案+小数化分数）

传送门

---

1.题目大意：输入n条线段，把每条线段变成原来的一个子线段，使得改变之后所有线段等长且不相交（端点可以重合）。求最大的相等长度

2.看到这里的最大，应该想到是二分答案，因为考虑贪心的话无从下手。但是这里需要注意的是我们如何枚举线段，这里需要用到前面讲的处理区间相交问题的贪心策略，给右端点排序，右端点相同的左端降序，这样之后设一个变量是当前所有线段的最右端，每次更新该值并当前线段在变成某长度是否合法，如下是浮点数二分的基本框架：

const double eps=1e-8;
double l=?,r=?,ans;
while(r-l>eps){
	double mid=(l+r)/2;
	if(check(...)){
		ans=mid;
		l=mid;  //r=mid;
	}else r=mid;  //l=mid;
}

3.关于小数化分数：

• 因为该分数表示为p/q=d，那么p=d*q，又因为p,q都是整数，因此枚举分母q，根据p=round(d*q)得到最精确的分子p，接着比较上次记录的整数对和这次哪个误差更小
• n的求法，因为1/x已经是以x为分母的最小值，考虑最坏的情况所有m条线段均在一个长度为1的区间内，那么最小的长度为1/m，因此枚举m个分母即可

const double eps=1e-8;
int n;
int ansp=0,ansq=1;  //分母初始不能为0
for(int p,q=1;q<=n;q++){  	//n是最大的分母
	p=round(ans*q);
	if(fabs((double)p/q-ans)<fabs((double)ansp/ansq-ans)){
		ansp=p,ansq=q;
	}
}

4.需要注意此题的eps设为1e-8会WA，此题要求的精度较高

代码：

#include <set>
#include <map>
#include <stack>
#include <queue>
#include <math.h>
#include <cstdio>
#include <string>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <unordered_map>
using namespace std;
#define lowbit(x) (x&(-x))
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> P;
const double eps=1e-9;
const double pi=acos(-1.0);
const int inf=0x3f3f3f3f;
const ll INF=1e18;
const int Mod=1e9+7;
const int maxn=2e5+10;

struct node{
    int l,r;
}a[maxn];
int n;
double ans;

bool cmp(node &p,node &q){
    return p.r==q.r?p.l>q.l:p.r<q.r;
}

bool check(double x){
    double cur=0.0;
    for(int i=0;i<n;i++){
        if(cur<a[i].l) cur=(double)a[i].l;
        if(cur+x>a[i].r){
            return false;
        }
        cur+=x;
    }
    return true;
}

int main(){
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    //ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);

    while(scanf("%d",&n)!=EOF){
        for(int i=0;i<n;i++)
            scanf("%d %d",&a[i].l,&a[i].r);
        sort(a,a+n,cmp);
        double l=0,r=1e6;
        while((r-l)>eps){
            double mid=(double)(l+r)/2;
            if(check(mid)){
                ans=mid;
                l=mid;
            }else r=mid;
        }
        //printf("%f",ans);
        int ansp=0,ansq=1;
        for(int p,q=1;q<=n;q++){
            p=round(ans*q);
            if(fabs((double)p/q-ans)<fabs((double)ansp/ansq-ans)){
                ansp=p,ansq=q;
            }
        }
        printf("%d/%d\n",ansp,ansq);
    }
    return 0;
}