Problem:
微生物実験(Bug Party)

あなたはJust Odd Inventions 社を知っているだろうか? この会社の業務は「ただ奇妙な発明(just odd inventions)」をすることである．ここでは略してJOI 社と呼ぶ．

JOI 社では多くの微生物を1 つのシャーレに生きたまま閉じ込める研究をしている．調査対象の微生物はN 匹存在し，1, 2, ... , N という番号がついている．各微生物はシャーレに閉じ込められると，foo (fatally odd object) と呼ばれる有害物質を一瞬のうちにシャーレ内に放出する．各微生物が放出するfoo の量は知られている．シャーレに閉じ込められた全微生物が放出したfoo はシャーレ内の各微生物が均等に摂取する．各微生物のfoo 許容量は知られており，この量を超えてfoo を摂取すると微生物は死んでしまう．

微生物i のfoo 放出量はai ミリグラム，foo 許容量はbi ミリグラムである．すなわち，シャーレに微生物i1, i2, ... , ik を閉じ込めたとき，シャーレ内の各微生物はそれぞれ(ai1 + ai2 + ... + aik)/k ミリグラムのfoo を摂取し，シャーレ内の微生物i は，この摂取量がbi より大きいと死んでしまう．

JOI 社からの委託により，あなたは出来るだけ多くの微生物をシャーレに生きたまま閉じ込めなければならない．ただし，微生物の死骸はシャーレ内の環境に悪影響を及ぼすため，シャーレ内のどの微生物もfoo の摂取で死んではならない．

なお，JOI 社が「ただ奇妙な発明」をすることでどうやって利益を得ているかは，今もって謎であり，JOI社内でも社長以外の誰も知らない．

課題

調査対象の微生物数と，各微生物のfoo 放出量とfoo 許容量が与えられたとき，1 つのシャーレに閉じ込めることができる微生物数の最大値を求めるプログラムを作成せよ．

制限

1 ≤ N ≤ 300000 = 3 × 105    調査対象の微生物の数

1 ≤ ai ≤ 100000 = 105    微生物i のfoo 放出量(ミリグラム)

1 ≤ bi ≤ 100000 = 105    微生物i のfoo 許容量(ミリグラム)

入力

標準入力から以下の入力を読み込め．

1 行目には整数N が書かれており，調査対象の微生物がN 匹存在することを表す．

続くN 行には各微生物の情報が書かれている．i + 1 行目(1 ≤ i ≤ N) には，正整数 ai, bi が空白を区切りとして書かれており，微生物i のfoo 放出量がai ミリグラムでfoo 許容量がbi ミリグラムであることを表す.

出力

標準出力に，1 つのシャーレに閉じ込めることができる微生物数の最大値を1 行で出力せよ．

注意

この問題では，扱う整数の範囲が32 ビットに収まらない可能性があることに注意せよ．

採点基準

採点用データのうち，配点の30% については，N ≤ 1000 を満たす．

入出力の例

入力例

6
12 8
5 9
2 4
10 12
6 7
13 9

出力例

3

この例では，微生物2, 4, 5 をシャーレに入れれば，放出されるfoo の合計量は5 + 10 + 6 = 21 ミリグラムとなり，それぞれの微生物が摂取するfoo の量は21/3 = 7 ミリグラムとなる．微生物2, 4, 5 のfoo 許容量はそれぞれ9, 12, 7 ミリグラムなので，シャーレ内のどの微生物も死なない．また，4 匹以上の微生物をどの微生物も死なないようにシャーレに入れることはできない．

問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C++ (judge verdict: Accepted)

#include <iostream>
#include <iomanip>
#include <vector>
#include <set>
#include <map>
#include <list>
#include <iterator>
#include <queue>
#include <stack>
#include <deque>
#include <sstream>
#include <algorithm>
#include <numeric>
#include <string>
#include <cmath>
#include <complex>
#include <functional>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cassert>
#include <cstdlib>
#include <cctype>

#define REP(i,m,n) for(int i=m;i<n;i++)
#define rep(i,n) for(int i=0;i<n;i++)
#define rp(i,c) rep(i,(int)c.size())
#define fr(i,c) for(__typeof((c).begin()) i=(c).begin();i!=(c).end();i++)
#define mp make_pair
#define pb push_back
#define all(c) (c).begin(),(c).end()
#define rall(u) (u).rbegin(),(u).rend()
#define dbg(x) cerr<<#x<<" = "<<(x)<<endl
#define v_delete(a,b) (a).erase(remove((a).begin(), (a).end(), b), (a).end())
#define v_unique(a) (a).erase(unique((a).begin(), (a).end()), (a).end())
#define VV(T) vector<vector< T > >

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef pair<int,int> pi;
const int inf=1<<28;
const double EPS=1e-9;
const long long INF = 100001000010000100ll;

template<class T>int SIZE(T a){return a.size();}
template<class T>string IntToString(T num){string res;stringstream ss;ss<<num;return ss.str();}
template<class T>T StringToInt(string str){T res=0;for(int i=0;i<SIZE(str);i++)res=(res*10+str[i]-'0');return res;}
template<class T>T gcd(T a,T b){if(b==0)return a;return gcd(b,a%b);}
template<class T>T lcm(T a,T b){return a/gcd(a,b)*b;}

template <class T> void input(T& a, int n) { for (int i = 0; i < n; ++i) cin >> a[i]; }
template <class T> void input(T* a, int n) { for (int i = 0; i < n; ++i) cin >> a[i]; }
typedef pair<int,int> pii;


bool compare(const pii& lhs, const pii& rhs){
    if((lhs.first - lhs.second) != (rhs.first - rhs.second))
        return (lhs.first - lhs.second) > (rhs.first - rhs.second);
    else
        return lhs.second < rhs.first;
}

int N;
vector<pii>d;
vi capal;

int maxh;

class BIT{
private:
    vector<ll> bit;
    int n;
public:
    BIT(){}
    BIT(int _n){
        n = _n;
        bit.assign(n+3, 0);
    }
    void init(int _n){
        n = _n;
        bit.assign(n+3, 0);
    }
    void add(int i, int v){
        for(;i <= n; i += (i & -i))bit[i]+=v;
    }
    ll sum(int i){
        ll s=0;
        for(;i > 0;i -= (i & -i))s+=bit[i];
        return s;
    }
};

BIT kosuu;
BIT hoshutu;

bool chk(int mid, int cap){
    return kosuu.sum(mid) * cap >= hoshutu.sum(mid);
}

int main(){
    cin >> N;
    rep(i, N){
        int a,b;
        cin >> a >> b;
        maxh=max(maxh, a);
        d.pb(mp(-b, a));
        //capal.pb(b);
    }
    kosuu.init(maxh);
    hoshutu.init(maxh);

    sort(all(d));
    //sort(all(capal), greater<int>());
    //capal.erase(unique(all(capal)), capal.end());
    int cur=0;
    ll ans=0;
    rp(i, d){
        //printf("info:%d %d\n", d[i].first, d[i].second);
        kosuu.add(d[i].second, 1);
        hoshutu.add(d[i].second, d[i].second);

        int ma=maxh+1;
        int mn=0;
        while(ma - mn > 1){
            int mid=(ma+mn)/2;
            if(chk(mid, -d[i].first)){
                mn=mid;
            }else{
                ma=mid;
            }
        }

        ll k1 = kosuu.sum(mn);
        ll k2 = kosuu.sum(mn+1);
        ll s = hoshutu.sum(mn);
        for(;k1 < k2;){
        	if((k1 + 1) * (-d[i].first) < s + (mn+1)){
        		break;
        	}else{
        		//printf("%d :: %d %d\n",k1, (k1 + 1) * -d[i].first , s + (mn+1));
        		k1++;
        		s += (mn+1);
        	}
        }
        if(ans < k1){
            ans = k1;
        }
        //puts("");
        //printf("%d %d\n", kosuu.sum(mn), hoshutu.sum(mn));

    }
    cout << ans << endl;


    return 0;

}