Time Limit: 2000MS | Memory Limit: 262144KB | 64bit IO Format: %I64d & %I64u |
Description
The administration of the Tomsk Region firmly believes that it‘s time to become a megacity (that is, get population of one million). Instead of improving the demographic situation, they decided to achieve its goal by expanding the boundaries of the city.
The city of Tomsk can be represented as point on the plane with coordinates (0; 0). The city is surrounded with n other locations, the i-th one has coordinates (xi, yi) with the population of ki people. You can widen the city boundaries to a circle of radius r. In such case all locations inside the circle and on its border are included into the city.
Your goal is to write a program that will determine the minimum radius r, to which is necessary to expand the boundaries of Tomsk, so that it becomes a megacity.
Input
The first line of the input contains two integers n and s (1 ≤ n ≤ 103; 1 ≤ s < 106) — the number of locatons around Tomsk city and the population of the city. Then n lines follow. The i-th line contains three integers — the xi and yi coordinate values of the i-th location and the number ki of people in it (1 ≤ ki < 106). Each coordinate is an integer and doesn‘t exceed 104 in its absolute value.
It is guaranteed that no two locations are at the same point and no location is at point (0; 0).
Output
In the output, print "-1" (without the quotes), if Tomsk won‘t be able to become a megacity. Otherwise, in the first line print a single real number — the minimum radius of the circle that the city needs to expand to in order to become a megacity.
The answer is considered correct if the absolute or relative error don‘t exceed 10 - 6.
Sample Input
Input
4 9999981 1 12 2 13 3 12 -2 1
Output
2.8284271
Input
4 9999981 1 22 2 13 3 12 -2 1
Output
1.4142136
Input
2 11 1 9999972 2 1
Output
-1
Source
Codeforces Round #242 (Div. 2)
#include <iostream> #include <stdio.h> #include <algorithm> #include <math.h> using namespace std; struct sa { int k; double r; }data[10050]; bool cmp(sa x,sa y) { return x.r<y.r; } int main() { int i,n,s,a,b; cin>>n>>s; for(i=0;i<n;i++) { cin>>a>>b>>data[i].k; data[i].r=sqrt(a*a+b*b); } sort(data,data+n,cmp); for(i=0;i<n;i++) { s+=data[i].k; if(s>=1000000) { printf("%.7lf",data[i].r); return 0; } } cout<<-1<<endl; return 0; }