2018 Multi-University Training Contest 6 L Pinball【受力分析】

最新推荐文章于 2019-08-09 11:05:27 发布

原创最新推荐文章于 2019-08-09 11:05:27 发布 · 481 阅读

1 ·

CC 4.0 BY-SA版权

Mathematics 专栏收录该内容

66 篇文章

订阅专栏

本文介绍了一种基于物理原理的运动模拟算法，该算法通过分解重力加速度来实现弹跳和匀加速直线运动的计算，并给出了具体的实现代码。通过计算不同方向上的运动周期，可以得到物体的运动轨迹。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

http://acm.hdu.edu.cn/contests/contest_showproblem.php?cid=807&pid=1012

将物体的重力加速度分解。
一部分用来做弹起运动，一部分用来做匀加速直线运动。
可以计算得出两个方向上的运动周期，两者之商为答案。

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
using namespace std;
typedef long long int LL;
const int N = 1e5 + 10;
const int mod = 998244353;
const double g = 9.8;
const double pi = acos(-1.0);
double t, vx, vy, v;
int T;
double a, b, x, y;
int main()
{
    scanf("%d",&T);
    while(T--){
        scanf("%lf%lf%lf%lf",&a,&b,&x,&y);
        t = sqrt(2.0*(y-b/a*(-x))/g);
        vx = 0; v = vy = g*t;
        int ans = 0;
        double da = atan(b/a), db = pi/2-da;
        vx = v * cos(fabs(db-da));
        vy = v * sin(fabs(db-da));
        double temp = -x;
        while(temp > 0){
            ans ++;
            if(db >= da) t = vy*2/g;
            else t = 0;
            double B = 2*(vy-b*vx/a), C = 2*b*vx*t/a;
            double tt = (-B + sqrt(B*B + 4*g*C))/2/g;
            t = t + tt;
            vy = vy + g*tt;
            temp -= vx * t;
            if(vx == 0) db = pi/2;
            else db = atan(vy/vx) - da;
            v = sqrt(vx*vx + vy*vy);
            vx = v * cos(fabs(db-da));
            vy = v * sin(fabs(db-da));
        }
        printf("%d\n", ans);
    }
    return 0;
}