题目描述
数轴上有n个闭区间:D1,...,Dn。
其中区间Di用一对整数[ai, bi]来描述,满足ai < bi。
已知这些区间的长度之和至少有10000。
所以,通过适当的移动这些区间,你总可以使得他们的“并”覆盖[0, 10000]——也就是说[0, 10000]这个区间内的每一个点都落于至少一个区间内。
你希望找一个移动方法,使得位移差最大的那个区间的位移量最小。
具体来说,假设你将Di移动到[ai+ci, bi+ci]这个位置。你希望使得maxi{|ci|} 最小。输入
输入的第一行包含一个整数n,表示区间的数量。
接下来有n行,每行2个整数ai, bi,以一个空格分开,表示区间[ai, bi]。
保证区间的长度之和至少是10000。
输出
输出一个数字,表示答案。如果答案是整数,只输出整数部分。如果答案不是整数,输出时四舍五入保留一位小数。
样例输入
2 10 5010 4980 9980样例输出
20
题目思路:很显然,我们要二分枚举结果,但是如何判断当前的结果行或不行是关键。
分析样例可以知道,答案可能存在0.5,1.5,2.5等。那么我们将所有数值扩大两倍,最后再除2就可以了。
接下来我们讨论应该采用什么贪心策略能够正确的判断当前枚举的数值行与不行。
- 区间乱序:假设有区间1~10,2-4,来覆盖1~12的区间。遍历的每个区间要保证覆盖左侧未覆盖的点,第一个区间移动距离为0,第二个区间的移动距离为8。分析可知,正确的结果应为2,说明这种策略是不合理的。
- 按照区间左端点对区间从小到大进行排序:假设有区间1~10,2-4,来覆盖1~12的区间。排序后仍为1~10,2~4;分析可知这个也是不合理的。
- 按照区间右端点对区间从小到大进行排序:假设有区间1~10,2-4,来覆盖1~12的区间。排序后为2~4,1~10。这种策略可以保证当前样例的正确性。从0开始更新覆盖的范围,每次使用一个可能的有边界尽量小的区间来更新覆盖围。以使右边界更大的区间剩余下来,尽量往右移动来覆盖更大的区域。注意:采用这种策略能够在一定程度上尽量保证结果的正确性,但是还会存在一定的特殊样例结果不对。因为这个题目的数据量比较小,所以我们可以采用多次遍历区间来解决,具体看代码。
采用第三种策略的实现代码:
判断函数:
bool judge(int x){//x为枚举的结果 int k=0;//未覆盖的最左端点 vector<node> tmp(vec); while(1){ int ok=0; for(int i=0;i<tmp.size() ;i++){//寻找可行区间,找到即break if(tmp[i].l-x<=k&&tmp[i].r+x>=k){//每个区间可以覆盖的范围 ok=1; if(k<=tmp[i].l+x)k=k+tmp[i].r-tmp[i].l;//k在区间左端点的可移动范围内 else k=tmp[i].r+x;//k不在 tmp.erase(tmp.begin() +i);//删除已经用过的区间 break;//每找到一个可行区间,结束 } } if(ok==0||k>=N)break;//未找到可行区间或所有区间已经覆盖 } return k>=N; }
AC代码:
//仔细观察样例我们发现,答案存在0.5, //那么我们将所有东西扩大两倍,最后再除2就可以了。 #include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <cmath> #include <map> #include <set> #include <stack> #include <queue> #include <vector> #include <string> #define cla(a, sum) memset(a, sum, sizeof(a)) #define rap(i, m, n) for(int i=m; i<=n; i++) #define rep(i, m, n) for(int i=m; i>=n; i--) #define bug printf("???\n") using namespace std; typedef long long ll; typedef unsigned long long ull; typedef pair<int, int> P; typedef pair<ll, ll> Pl; const int Inf = 0x3f3f3f3f; const double eps = 1e-8; const int mod=998244353; const int maxn = 1e4+5; const int N=2e4; template <typename T> void read(T &x){ x = 0; int f = 1; char ch = getchar(); while (!isdigit(ch)) {if (ch == '-') f = -1; ch = getchar();} while (isdigit(ch)) {x = x * 10 + ch - '0'; ch = getchar();} x *= f; } int n; struct node{ int l; int r; }; vector<node>vec; bool cmp(node u,node v){ return u.r <v.r ; } bool judge(int x){ int k=0; vector<node> tmp(vec); while(1){ int ok=0; for(int i=0;i<tmp.size() ;i++){ if(tmp[i].l-x<=k&&tmp[i].r+x>=k){ ok=1; if(k<=tmp[i].l+x)k=k+tmp[i].r-tmp[i].l; else k=tmp[i].r+x; tmp.erase(tmp.begin() +i); break; } } if(ok==0||k>=N)break; } return k>=N; } int main() { //ios::sync_with_stdio(false),cin.tie(0); scanf("%d",&n); int l,r; rap(i,1,n){ scanf("%d%d",&l,&r); l=2*l;r=2*r; vec.push_back({l,r}); } sort(vec.begin() ,vec.end() ,cmp); l=0,r=N; double ans=0; while(l<=r){ int mid=(l+r)>>1; if(judge(mid)){ r=mid-1; ans=mid; } else l=mid+1; } ans=ans/2.0; cout<<ans<<endl; return 0; }