結果
問題 | No.2010 Magical Floor |
ユーザー |
![]() |
提出日時 | 2022-07-15 22:40:43 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 15 ms / 2,000 ms |
コード長 | 1,155 bytes |
コンパイル時間 | 1,790 ms |
コンパイル使用メモリ | 169,300 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-27 19:34:34 |
合計ジャッジ時間 | 2,619 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 21 |
ソースコード
#include <bits/stdc++.h>using namespace std;typedef long long int ll;typedef pair<ll,ll> P;typedef vector<ll> VI;typedef vector<VI> VVI;#define REP(i,n) for(ll i=0;i<(n);i++)#define ALL(v) v.begin(),v.end()template<typename T> bool chmax(T &x, const T &y) {return (x<y)?(x=y,true):false;};template<typename T> bool chmin(T &x, const T &y) {return (x>y)?(x=y,true):false;};constexpr ll MOD=1000000007;constexpr ll INF=2e18;double f(ll h, ll a, ll b){double x=(double)b*b*h*h/(a*a-b*b);if(x<0)return INF;elsereturn sqrt(x);}void solve(){ll n, x; cin >> n >> x;VI h(n+1), a(n); h[0]=0;REP(i,n) cin >> h[i+1];REP(i,n) cin >> a[i];double ans=INF;REP(i,n){double d=0, msum=0;for(int j=i-1;j>=0;j--){double res=f(h[j+1]-h[j],a[j],a[i]);msum+=res;d+=sqrt(res*res+(h[j+1]-h[j])*(h[j+1]-h[j]))*a[j];}if(msum*2>x)continue;d*=2;d+=(x-msum*2)*a[i];chmin(ans,d);}cout << setprecision(10) << ans << endl;}int main(){int t; cin >> t;while(t--) solve();}