粒子 - 二分

Mys_C_K

于 2018-12-20 20:50:41 发布

阅读量177

点赞数

分类专栏：二分

本文链接：https://blog.csdn.net/Mys_C_K/article/details/85139365

版权

二分专栏收录该内容

5 篇文章 0 订阅

订阅专栏

涨姿势了为啥每次运算强制类型转化(double)会比一开始就全开double慢一个数量级啊

题目大意：
数轴上有两个起点，分别发射正反粒子各N，正反粒子相撞会湮灭，保证任意一个粒子都会碰到另一个相反的粒子，相同的粒子相撞不会有任何问题。问前k对相撞的粒子是？ $n\le4e5,k\le100$
题解：
二分出下一次有粒子相撞的时刻T。
做k次即可。

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define db double
#define mp make_pair
#define fir first
#define sec second
#define gc getchar()
#define debug(x) cerr<<#x<<"="<<x
#define sp <<" "
#define ln <<endl
using namespace std;
typedef pair<db,db> pdd;
inline int inn()
{
    int x,ch;while((ch=gc)<'0'||ch>'9');
    x=ch^'0';while((ch=gc)>='0'&&ch<='9')
        x=(x<<1)+(x<<3)+(ch^'0');return x;
}
const int N=50010;const db eps=1e-3;
pair<pdd,int> A[N],B[N];int dA[N],dB[N],n;db L;
inline int check(db T)
{
    db mxA=0,mxB=0;
    rep(i,1,n)
    {
        if(!dA[i]) mxA=max(mxA,(T-A[i].fir.fir)*A[i].fir.sec);
        if(!dB[i]) mxB=max(mxB,(T-B[i].fir.fir)*B[i].fir.sec);
        if(mxA+mxB>=L) return 1;
    }
    return 0;
}
db low;
inline int solve()
{
    db l=low,r=0,mid=(l+r)*0.5;
    r=max(r,A[n].fir.fir+L/A[n].fir.sec);
    r=max(r,B[n].fir.fir+L/B[n].fir.sec);
    while(l<r-eps)
    {
        if(!check(mid)) l=mid;
        else r=mid;mid=(l+r)*0.5;
    }
    low=l;
    int x=0,y=0;db mxA=-1,mxB=-1,v;
    rep(i,1,n)
    {
        if(!dA[i]&&(v=(l-A[i].fir.fir)*A[i].fir.sec)>mxA) mxA=v,x=i;
        if(!dB[i]&&(v=(l-B[i].fir.fir)*B[i].fir.sec)>mxB) mxB=v,y=i;
    }
    return printf("%d %d\n",A[x].sec,B[y].sec),dA[x]=dB[y]=1;
}
inline bool cmp(const pair<pdd,int> &a,const pair<pdd,int> &b)
{   return a.fir.fir+L/a.fir.sec<b.fir.fir+L/b.fir.sec;  }
int main()
{
    n=inn(),L=inn();int k=inn();
    rep(i,1,n) A[i].fir.fir=inn(),A[i].fir.sec=inn(),A[i].sec=i;
    rep(i,1,n) B[i].fir.fir=inn(),B[i].fir.sec=inn(),B[i].sec=i;
    sort(A+1,A+n+1,cmp),sort(B+1,B+n+1,cmp);
    while(k--) solve();return 0;
}