BZOJ3570: DZY Loves Physics I

3570: DZY Loves Physics I

Time Limit: 10 Sec  Memory Limit: 256 MB
Submit: 101  Solved: 64
[Submit][Status]

Description

背景

众所周知，DZY是个大学霸，精通数理化。有天，吉丽拿着一道物理题目去问DZY，DZY很快就秒了这题，但是懒得算了，就让你来解决它。

题目描述

现在水平面上有一条无限长的光滑轨道，上面有n个小球(体积忽略)，每个小球的质量都相同。现在这些小球在某些神奇的力量下开始了运动。 给定一个常数C, 每个小球在任意时刻的速度v和加速度a均满足a*v=C. 

吉丽还规定，在任意时刻，每个小球的加速度a和速度v方向相同。

因为有多个小球的存在，在运动中，小球之间可能发生碰撞，这里假设小球间发生的碰撞均为完全弹性碰撞。 

吉丽还会给你常数C，每个小球最开始的位置，运动方向，以及初速度Vi.

吉丽的询问是这样的： 给你整数T,K, 让你求当前系统中，让小球开始运动T时刻后，速率第K小的小球的速率为多少。

注意，每次询问之后，小球回到原位。

当然，为了考察你的随机应变能力，吉丽会在问问题的时候不时地增加新的小球。

Input

第一行包含两个正整数n,C. 

接下来n行，每行包含三个整数Vi,Xi,Di

Vi代表第i个小球的速率

Xi代表第i个小球的初始位置(你可以想象轨道是一个数轴)

Di代表第i个小球的初始运动方向 其中1表示正方向，-1表示反方向。

接下来一行Q 代表操作的个数

接下来Q行，第一个数为操作的种类tp

0 Vj Xj Dj 代表吉丽增加了一个初速度Vj, 位置Xj， 运动方向Dj的小球

1 T K 代表吉丽的一次询问

 

Output

对于每一个询问，输出一行包含一个实数，代表答案。 保留3位小数。

数据范围

n,m<=10^5 Vi<=10^5 |Xi|<=10^9 0<C<=10^9 0<=T<=10^5 

保证数据合法

 

Sample Input

3 10
1 2 1
3 4 1
6 5 1
3
1 5 1
0 8 -5 -1
1 8 2

Sample Output

10.050
13.000

HINT

n,Q<=10^5 Vi<=10^5 |Xi|<=10^9 0<C<=10^9 0<=T<=10^5 

Source

By AcrossTheSky

题解：

经过……推导，然后我们发现只需要求v的第k小就行了。。。

主席树，平衡树？大材小用。

然后用权值线段树是很容易在lgn的时间内求出第k小的

但我们发现树状数组居然也可以。。。

原理自行脑补吧。。。

代码：

#include<cstdio>

#include<cstdlib>

#include<cmath>

#include<cstring>

#include<algorithm>

#include<iostream>

#include<vector>

#include<map>

#include<set>

#include<queue>

#include<string>

#define inf 1000000000

#define maxn (1<<17)

#define maxm 20000000+5

#define eps 1e-10

#define ll long long

#define pa pair<int,int>

#define for0(i,n) for(int i=0;i<=(n);i++)

#define for1(i,n) for(int i=1;i<=(n);i++)

#define for2(i,x,y) for(int i=(x);i<=(y);i++)

#define for3(i,x,y) for(int i=(x);i>=(y);i--)

#define mod 1000000007

using namespace std;

inline int read()

{

    int x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}

    return x*f;

}
int n,q,c,s[maxn+1];
inline void add(int x){for(;x<=maxn;x+=x&(-x))s[x]++;}

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);
    n=read();c=read();
    while(n--)add(read()),read(),read();
    q=read();
    while(q--)
    {
        int ch=read();
        if(!ch)add(read()),read(),read();
        else
        {
            int l=0,r=maxn,t=read(),k=read(),sum=0;
            while(l<r)
            {
                int mid=(l+r)>>1;sum+=s[mid];
                //cout<<l<<' '<<mid<<' '<<r<<' '<<sum<<' '<<k<<endl;
                if(sum>=k)r=mid,sum-=s[mid];else l=mid+1;
            }
            printf("%.3f
",(double)sqrt((ll)2*c*t+(ll)l*l));
        }
    }

    return 0;

}