1112 KGold

给出N个人在0时刻的财富值M[i](所有人在0时刻的财富互不相等),以及财富增长速度S[i],随着时间的推移,某些人的财富值会超越另外一些人。如果时间足够长,对于财富增长最快的人来说,他的财富将超越所有其他对手。
求发生的前10000次超越,分别是谁超过了谁?如果总的超越次数不足10000,则输出所有超越,如果1次超越都不会发生,则输出No Solution。
输出按照超越发生的时间排序,同一时刻发生的超越,按照超越者的编号排序,如果编号也相同,则按照被超越者的编号排序。所有排序均为递增序。
Input
第1行:N,表示人的数量。(1 <= N <= 10000)
第2 - N + 1行:每行2个数,分别是初始的财富值M[i],和财富增长速度S[i]。(0 <= M[i] <= 10^5, 1 <= S[i] <= 100)
Output
输出前10000次超越,超越者和被超越者的编号。如果总的超越次数不足10000,则输出所有。如果1次超越都不会发生,则输出No Solution。
输出按照超越发生的时间排序,同一时刻发生的超越,按照超越者的编号排序,如果超越者编号也相同,则按照被超越者的编号排序。所有排序均为递增序。
Input示例
4
1 100
2 100
3 1
4 50
Output示例
2 3
1 3
2 4
1 4

题解:先用二分找到10000次超越所在的时间up,然后枚举两条直线判断相交点的时间是否up以内,若是,放入优先队列中维护,输出优先队列前10000个元素即可。

代码:

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <cstring>
#include <string>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <vector>
#include <set>
#include <stack>
#include <map>
#include <cmath>
#include <ctime>
using namespace std;
typedef long long ll;
typedef pair<ll, int> P;
const int INF = 0x3f3f3f3f;
const ll LINF = 0x3f3f3f3f3f3f3f3f;
const double PI = acos(-1.0);
const double eps = 1e-10;
const int maxn = 1e4+7;
const int mod = 1e9+7;

int n;
struct G{
    int m,s,id;
    G(int m,int s,int id):m(m),s(s),id(id){}
    G(int m,int id):m(m),id(id){}
    G(){}
    bool operator < (const G &b)const{
        return m<b.m;
    }
} q[maxn],p[maxn];

struct U{
    double tim;
    int x,y;
    U(double tim,int x,int y):tim(tim),x(x),y(y){}
//    bool operator <(const U &b)const
//    {
//        return fabs(tim-b.tim)<eps?(x==b.x?y<b.y:x<b.x):tim<b.tim;
//    }
    bool operator <(const U &b)const{
        return fabs(tim-b.tim)<eps?(x==b.x?y>b.y:x>b.x):tim>b.tim;
    }
    U(){}
};
priority_queue<U> que;
int check(double a){
    int sum = 0;
    for(int i = 1;i <= n;i++){
        p[i] = G(q[i].m+q[i].s*a,i);
    }
    sort(p+1,p+n+1);
    for(int i = 1;i <= n;i++){
        if(i>p[i].id){
            sum+=i-p[i].id;
        }
    }
    return sum;
}

int main()
{
    scanf("%d",&n);
    for(int i = 1;i <= n;i++){
        int a,b;
        scanf("%d %d",&a,&b);
        q[i] = G(a,b,i);
    }
    sort(q+1,q+n+1);
    double l = 0,r = 1e6+7;
    while(r-l>eps){
        double mid = (l+r)/2;
        if(check(mid)>10000){
            r = mid;
        }
        else l = mid;
    }
    double up = l;
    double nx = 0;
    for(int i = 2;i <= n;i++){
        for(int j = 1;j < i;j++){
            nx = 1.0*(q[i].m-q[j].m)/(q[j].s-q[i].s);
            if(nx>eps&&nx<=up){
                que.push(U(nx,q[j].id,q[i].id));
                if(que.size()>=20000) break;
            }
        }
        if(que.size()>=20000) break;
    }
    if(!que.size()) return puts("No Solution"),0;
    int cnt = 0;
    while(que.size()&&cnt<10000){
        U now = que.top();
        que.pop();
        printf("%d %d\n",now.x,now.y);
        cnt++;
    }
    return 0;
}

总结:一开始以为是n次方的增长的,想来是曲线,变没什么想法,后来知道是直线,思路也就很简单了。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值