彈性碰撞·POJ-3684-Physics Experiment

題目大意：

n個小球在一圓筒內每隔一秒垂直下落，問t秒後n個小球的位置。

解題思路：

如果只有一個小球，那麼t時間小球的位置很好確定。

但是如果多個，像以前Ants那道題一樣就算碰觸也不算他反彈，直接視作繼續前進，但是這裏有個問題就是小球是有半徑的，就算是一起釋放上方小球也要多上2*r，所以除卻第一個小球，其餘小球都比下一個多2*r，但是因爲我們是要按照編號來輸出，我們按照ants的做法得不到對應編號的小球位置，但是我們可以知道一件事情 ：下面的小球的位置一定比上面的小球低，因此最後排序就行了。

但是我在做的時候卡到了兩件事情上：

1、沒考慮時間小於n的情況，這個時候t時間內未釋放的小球的高度就是H。

2、我竟然忘記乘上i了；

 H[i] + 2 * r * i/ 100.0

AC代碼：

#include <map>
#include <set>
#include <stack>
#include <cmath>
#include <queue>
#include <bitset>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <fstream>
#include <cstdlib>
#include <sstream>
#include <cstring>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;
#define   maxn          1010
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   ms(x,y)      memset(x,y,sizeof(x))
#define   rep(i,n)      for(int i=0;i<(n);i++)
#define   repf(i,a,b)   for(int i=(a);i<=(b);i++)
#define   PI           pair<int,int>
//#define   mp            make_pair
#define   FI            first
#define   SE            second
#define   IT            iterator
#define   PB            push_back
#define   Times         10
typedef   long long     ll;
typedef   unsigned long long ull;
typedef   long double   ld;
typedef   pair<int ,int > P;
//#define N 100
const double eps = 1e-10;
const double  pi = acos(-1.0);
const  ll    mod = 1e9+7;
const  int   inf = 0x3f3f3f3f;
const  ll    INF = (ll)1e18+300;
const int   maxd = 1000000 + 10;

#define g 10.0
int n;
double h, r, t;
double H[maxd];

double solve(double time) {
    if(time < 0) {
        return h;
    }
    double t = sqrt(2 * h / g);
    int k = (int)(time / t);
    if(k&1) {
        double temp = (k + 1) * t - time;
        return h - g / 2 * temp * temp;
    }
    else {
        double temp = time - k * t;
        return h - g / 2 * temp * temp;
    }
}
int main() {
    int T;
    cin >> T;
    while(T --) {
        cin >> n >> h >> r >> t;
        for (int i = 0; i < n; i++) {
            H[i] = solve(t - i);
        }
        sort(H, H + n);
        for (int i = 0; i < n; i++) {
            if(i != 0) {
                cout << " ";
            }
            if(i == 0) {
                printf("%.2f", H[i]);
                //cout << H[i] ;
            }
            else {
                printf("%.2f", H[i] + 2 * r * i/ 100.0);
                //cout << H[i] + 2 * r / 100.0 ;
            }
        }
        cout << endl;
    }
}