結果
問題 | No.1413 Dynamic Sushi |
ユーザー | tanimani364 |
提出日時 | 2021-04-04 11:29:59 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 2,899 bytes |
コンパイル時間 | 2,037 ms |
コンパイル使用メモリ | 208,308 KB |
実行使用メモリ | 10,880 KB |
最終ジャッジ日時 | 2024-06-08 02:52:03 |
合計ジャッジ時間 | 16,970 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 3 ms
9,472 KB |
testcase_01 | TLE | - |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | TLE | - |
testcase_05 | TLE | - |
testcase_06 | TLE | - |
testcase_07 | TLE | - |
testcase_08 | TLE | - |
testcase_09 | TLE | - |
testcase_10 | TLE | - |
testcase_11 | TLE | - |
testcase_12 | TLE | - |
testcase_13 | -- | - |
testcase_14 | -- | - |
testcase_15 | -- | - |
testcase_16 | -- | - |
testcase_17 | -- | - |
testcase_18 | -- | - |
testcase_19 | -- | - |
testcase_20 | -- | - |
testcase_21 | -- | - |
testcase_22 | -- | - |
testcase_23 | -- | - |
testcase_24 | -- | - |
ソースコード
#include <bits/stdc++.h> //#include<boost/multiprecision/cpp_int.hpp> //#include<boost/multiprecision/cpp_dec_float.hpp> //#include <atcoder/all> #define rep(i, a) for (int i = (int)0; i < (int)a; ++i) #define rrep(i, a) for (int i = (int)a - 1; i >= 0; --i) #define REP(i, a, b) for (int i = (int)a; i < (int)b; ++i) #define RREP(i, a, b) for (int i = (int)a - 1; i >= b; --i) #define repl(i, a) for (ll i = (ll)0; i < (ll)a; ++i) #define pb push_back #define eb emplace_back #define all(x) x.begin(), x.end() #define rall(x) x.rbegin(), x.rend() #define popcount __builtin_popcount #define fi first #define se second using ll = long long; constexpr ll mod = 1e9 + 7; constexpr ll mod_998244353 = 998244353; constexpr ll INF = 1LL << 60; // #pragma GCC target("avx2") // #pragma GCC optimize("O3") // #pragma GCC optimize("unroll-loops") //using lll=boost::multiprecision::cpp_int; //using Double=boost::multiprecision::number<boost::multiprecision::cpp_dec_float<1024>>;//仮数部が1024桁 template <class T> inline bool chmin(T &a, T b) { if (a > b) { a = b; return true; } return false; } template <class T> inline bool chmax(T &a, T b) { if (a < b) { a = b; return true; } return false; } template <typename T> T mypow(T x, T n, const T &p = -1) { //x^nをmodで割った余り if (p != -1) { x %= p; } T ret = 1; while (n > 0) { if (n & 1) { if (p != -1) ret = (ret * x) % p; else ret *= x; } if (p != -1) x = (x * x) % p; else x *= x; n >>= 1; } return ret; } using namespace std; //using namespace atcoder; void solve() { int n; double w; cin>>n>>w; vector<double>x(n),y(n),r(n),v(n),a(n); rep(i,n)cin>>x[i]>>y[i]>>r[i]>>v[i]>>a[i]; auto dis=[](pair<double,double> p,pair<double,double> q)->double{ return (p.first-q.first)*(p.first-q.first)+(p.second-q.second)*(p.second-q.second); }; auto calc=[](double t,double x,double y,double r,double v,double a)->pair<double,double>{ double pi=acos(-1); return {x+r*cos((v*t+a)*pi/180),y+r*sin((v*t+a)*pi/180)}; }; constexpr double inf=1e9; vector<vector<double>>dp((1<<n)+5,vector<double>(n+1,inf)); dp[0][n]=0;//nは座標(0,0) rep(bit,1<<n){ rep(j,n+1){ rep(k,n){ if(j!=n&&!(bit&1<<j))continue; if(bit&1<<k)continue; double pre_t=dp[bit][j]; pair<double,double> pre_pos={0,0}; if(j!=n)pre_pos=calc(pre_t,x[j],y[j],r[j],v[j],a[j]); double t=0; double left=pre_t,right=1e3; rep(cnt,500){ double mid=(left+right)/2; auto cur_pos=calc(mid,x[k],y[k],r[k],v[k],a[k]); if(dis(pre_pos,cur_pos)<=w*w*(mid-pre_t)*(mid-pre_t))right=mid; else left=mid; } t=right; chmin(dp[bit|1<<k][k],t); } } } double ans=inf; rep(i,n)chmin(ans,dp[(1<<n)-1][i]); cout<<ans<<"\n"; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(15); solve(); return 0; }