Problem:
Problem E: Immortal Jewels

English text is not available in this practice contest.

ある時ある貴族が，ある貧乏な国のおてんばで勇敢なお姫様に惚れ込み結婚を申し込んだ．お姫様は貴族にある条件を出した．その条件とは「不死の宝石」と呼ばれている宝石を大量に持ってくることであった．不死の宝石は，ある山の特定の場所でしか取ることができない非常に希少な宝石である．しかも大変壊れやすいため，採取するには特別な方法が必要だった．

不死の宝石は円形をしており，二次元空間上に複数存在している．この宝石を取るには，ある特別な金属の棒でそれらを吸着する必要がある．金属の棒は無限の長さを持つ直線であり，太さは無視することができる．宝石は一つ一つ異なる強さの磁力を持っており，金属がその磁力に反応するほど十分に近ければ宝石は吸着される．具体的には，金属と宝石の表面との距離を d ，宝石の磁力の強さを m としたとき，

0 ≤ d ≤ m

であれば宝石は金属に吸着される．逆に，金属の棒と宝石がその磁力よりも離れている場合は吸着できない．また，棒が少しでも宝石を貫通してしまった場合にも，その宝石が壊れてしまうため吸着できない．

例を見てみよう．下の図は二次元空間上に置かれた宝石の例である．宝石1から宝石6まであり，磁力はそれぞれ1, 0, 1, 1, 1, 2であるとする．

図 E-1: 宝石の配置例

上の図に加えて金属の棒を配置した例が下の図である．宝石を吸着した結果も表に示してある．この例の場合，宝石3は磁力の届く範囲より離れており，また宝石4は棒に貫通されているため吸着できないが，残りの4つは全て吸着することができる．

図 E-2: 金属の棒の配置例

宝石名

磁力

金属との距離

吸着できるか

宝石1

1

約0.21

できる

宝石2

0

0

できる

宝石3

1

約5.37

できない

宝石4

1

貫通している

できない

宝石5

1

約0.97

できる

宝石6

2

約0.53

できる

表 E-3: 吸着結果

貴族は全財産を注ぎ込み，特別な金属の棒を必死に探し求めた．しかしながら，この金属も非常に貴重であったため，結局たった一本しか入手することができなかった．したがって吸着のチャンスはただ一回しか無い．

あなたはある貴族に仕えるプログラマーである．あなたの仕事は，与えられた二次元上の宝石の配置に対して，金属の棒をうまく配置したときに，最大で何個の宝石を吸着することができるかを求めるプログラムを書くことである．

Input

入力は複数のデータセットから成り， 1つのデータセットは以下の形式で与えられる．

N

x1 y1 r1 m1

x2 y2 r2 m2

...

xN yN rN mN

データセットの最初の行は宝石の数 N (1 ≤ N ≤ 50) を表している．続く N 行の各行には4つの整数 xi, yi, ri, mi (-1000 ≤ xi, yi ≤ 1000, 1 ≤ ri ≤ 100, 0 ≤ mi ≤ 100) が記述されており，宝石の位置，大きさ，および磁力を表す．すなわち，宝石 i は中心を(xi, yi)とし半径が ri であるような円形をしており，その磁力は mi である．宝石は互いに重なり合うことは無い．

入力の終わりは，0のみからなる行で表される．

Output

各データセットについて， 一度に吸着することができる宝石の最大数を1行に出力せよ．

Sample Input

6
-2 -2 1 1
2 2 2 0
5 7 1 1
8 0 3 1
13 4 1 1
16 1 1 2
3
0 0 2 1
10 0 2 1
0 10 2 1
3
0 0 2 1
10 0 2 1
0 6 2 1
3
0 0 2 1
10 0 2 1
0 4 2 1
1
0 0 1 1
0

Output for the Sample Input

4
2
3
3
1

Solution in C++11 (judge verdict: Accepted)

#include <iostream>
#include <vector>
#include <algorithm>
#include <cassert>
constexpr double eps = 1e-7;
namespace boost {template <class T>struct integer_iterator {T a;bool operator != (integer_iterator const & it) const { return a != it.a; }T operator * () const { return a; }integer_iterator & operator ++ () { ++ a; return *this; }integer_iterator operator ++ (int) { return { a ++ }; }};template <class T>struct integer_range {T l, r;typedef integer_iterator<T> iterator;iterator begin() const { return { l }; }iterator end  () const { return { r }; }};template <class T>integer_range<T> irange(T l, T r) { assert (l <= r); return { l, r }; }template <class T>struct integer_iterator_with_step {T a, d, i;bool operator != (integer_iterator_with_step const & it) const { return a != it.a or d != it.d or i != it.i; }T operator * () const { return a+d*i; }integer_iterator_with_step & operator ++ () { ++ i; return *this; }integer_iterator_with_step operator ++ (int) { return { a, d, i ++ }; }};template <class T>struct strided_integer_range {T l, r, s;typedef integer_iterator_with_step<T> iterator;iterator begin() const { return { l, s, 0 }; }iterator end  () const { return { l, s, (r - l) / s }; }};template <class T>strided_integer_range<T> irange(T l, T r, T s) { assert (s > 0 ? l <= r : s < 0 ? l >= r : false); return { l, r, s }; }}
#include <utility>
#include <cmath>
#define ttt template <typename T>
#define pca point<T> const &
ttt struct point { T x, y; };ttt point<T> operator + (pca a, pca b) { return { a.x+b.x, a.y+b.y }; }ttt point<T> operator - (pca a, pca b) { return { a.x-b.x, a.y-b.y }; }ttt point<T> operator - (pca a) { return { -a.x, -a.y }; }ttt point<T> operator * (T a, pca b) { return { a*b.x, a*b.y }; }ttt std::pair<T,T> to_pair(pca a) { return { a.x, a.y }; }ttt bool operator == (pca a, pca b) { return to_pair(a) == to_pair(b); }ttt bool operator != (pca a, pca b) { return to_pair(a) != to_pair(b); }ttt bool operator <  (pca a, pca b) { return to_pair(a) <  to_pair(b); }ttt bool operator <= (pca a, pca b) { return to_pair(a) <= to_pair(b); }ttt bool operator >= (pca a, pca b) { return to_pair(a) >= to_pair(b); }ttt bool operator >  (pca a, pca b) { return to_pair(a) >  to_pair(b); }ttt T length_squared(pca p) { return p.x*p.x + p.y*p.y; }ttt double length(pca p) { return sqrt(length_squared(p)); }ttt point<T> normalized(pca a) { return (1 / length(a)) * a; }ttt T   dot(pca p, pca q) { return p.x * q.x + p.y * q.y; }ttt T cross(pca p, pca q) { return p.x * q.y - p.y * q.x; }ttt int ccw(pca a, pca b, pca c) { T x = cross(b - a, c - a); return x > 0 ? 1 : x < 0 ? -1 : 0; }
#undef ttt
#undef pca
template <typename T> std::istream & operator >> (std::istream & input, point<T> & p) { return input >> p.x >> p.y; }template <typename T> std::ostream & operator << (std::ostream & output, const point<T> & p) { return output << p.x << ' ' << p.y; }
template <typename T> struct circle { point<T> p; T r; };
template <typename T> struct line { point<T> a, b; };template <typename T> struct segment { point<T> a, b; };template <typename T> T length_squared(const segment<T> & p) { return length_squared(p.a - p.b); }template <typename T> double length(const segment<T> & p) { return length(p.a - p.b); }
template <typename T>T sq(T x) { return x * x; }
template <typename T>point<T> rotate_90(point<T> const & a) {return { - a.y, a.x };}template <typename T>point<T> rotate_270(point<T> const & a) {return { a.y, - a.x };}template <typename T>point<T> rotate_180(point<T> const & a) {return { - a.x, - a.y };}point<double> rotate(point<double> const & a, double th) {return {a.x * cos(th) - a.y * sin(th),a.x * sin(th) + a.y * cos(th) };}
bool does_include(circle<double> const & a, point<double> const & b) {return length(a.p - b) < a.r + eps;}bool is_on_circumference(point<double> const & a, circle<double> const & b) {return std::abs(length(b.p - a) - b.r) < eps;}// commutable
bool does_include(circle<double> const & a, circle<double> const & b) {return length(a.p - b.p) < std::abs(a.r - b.r) + eps;}// commutable
bool does_intersect(circle<double> const & a, circle<double> const & b) {return length(a.p - b.p) < a.r + b.r + eps;}// commutable
bool does_circumscribe(circle<double> const & a, circle<double> const & b) {return std::abs(length(a.p - b.p) - (a.r + b.r)) < eps;}// commutable
bool does_inscribe(circle<double> const & a, circle<double> const & b) {return std::abs(length(a.p - b.p) - std::abs(a.r - b.r)) < eps;}std::vector<point<double> > intersection_points(circle<double> const & a, circle<double> const & b) {std::vector<point<double> > result;if (not does_intersect(a,b)) {// nop
} else if (does_include(a,b) and not does_inscribe(a,b)) {// truly included
// nop
} else if (does_inscribe(a,b) or does_circumscribe(a,b)) {result.push_back(a.r * normalized(b.p - a.p));} else {//       c
//      /|        // a.r / |  b.r
//    /th|          //   a-------b
//       l
double l = length(b.p - a.p);double ac = ( sq(a.r) - sq(b.r) + sq(l) ) / ( 2 * l ); // a * cos(th)
double as = sqrt( sq(a.r) - sq(ac) );point<double> e = normalized(b.p - a.p);result.push_back(a.p + ac*e + as*rotate_90(e));result.push_back(a.p + ac*e + as*rotate_270(e));}return result;}std::vector<line<double> > common_tangents(circle<double> const & a, point<double> const & b) {std::vector<line<double> > result;if (does_include(a,b) and not is_on_circumference(b,a)) {// truly included
// nop
} else if (is_on_circumference(b,a)) {point<double> e = normalized(a.p - b);result.push_back({b + rotate_90(e),b + rotate_270(e) });} else {double theta = asin(a.r / length(b - a.p));result.push_back({ b, b + rotate(a.p - b, + theta) });result.push_back({ b, b + rotate(a.p - b, - theta) });}return result;}std::vector<line<double> > common_inner_tangents(circle<double> const & a, circle<double> const & b) {point<double> p = a.p + (a.r / (a.r + b.r)) * (b.p - a.p);return common_tangents(a,p);}std::vector<line<double> > common_outer_tangents(circle<double> a, circle<double> b) {if (b.r < a.r) std::swap(a,b);std::vector<line<double> > result;if (does_include(a,b) and not does_inscribe(b,a)) {// truly included
// nop
} else if (does_inscribe(b,a)) {point<double> e = normalized(a.p - b.p);result.push_back({b.p + a.r*e + rotate_90(e),b.p + a.r*e + rotate_270(e) });} else {point<double> d = b.p - a.p;double theta = asin((b.r - a.r) / length(d));point<double> dp = rotate(d, + theta);point<double> dm = rotate(d, - theta);point<double> ep = rotate_90(normalized(dp));point<double> em = rotate_270(normalized(dm));result.push_back({ a.p + a.r*ep, a.p + a.r*ep + dp });result.push_back({ a.p + a.r*em, a.p + a.r*em + dm });}return result;}std::vector<line<double> > common_tangents(circle<double> const & a, circle<double> const & b) {std::vector<line<double> > result;for (auto l : common_inner_tangents(a,b)) result.push_back(l);for (auto l : common_outer_tangents(a,b)) result.push_back(l);return result;}double distance(line<double> const & a, point<double> const & b) {point<double> p = b   - a.a;point<double> q = a.b - a.a;return std::abs( cross(p, q) / length(q) );}
using namespace std;
using namespace boost;
int main() {
    int age = 0;
    while (true) {
        int n; cin >> n;
        if (n == 0) break;
        vector<point<double> > p(n); vector<double> r(n), m(n);
        for (int i : irange(0,n)) cin >> p[i] >> r[i] >> m[i];
        if (n <= 2) {
            cout << n << endl;
            continue;
        }
        vector<circle<double> > cs(2*n);
        for (int i : irange(0,n)) {
            cs[2*i].p = p[i];
            cs[2*i].r = r[i];
            cs[2*i+1].p = p[i];
            cs[2*i+1].r = r[i]+m[i];
        }
        vector<line<double> > ls;
        for (int i : irange(0,2*n)) {
            for (int j : irange(i+1,2*n)) {
                for (auto l : common_tangents(cs[i], cs[j])) {
                    ls.push_back(l);
                }
            }
        }
        int result = 0;
        for (auto l : ls) {
            int jewels = 0;
            for (int i : irange(0,n)) {
                double d = distance(l, p[i]);
                if (r[i] < d + eps and d < r[i]+m[i] + eps) {
                    jewels += 1;
                }
            }
            result = max(result, jewels);
        }
        cout << result << endl;
    }
    return 0;
}