QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#166824 | #7110. Kuririn MIRACLE | ucup-team1479 | WA | 86ms | 4388kb | C++14 | 1.5kb | 2023-09-06 19:01:28 | 2023-09-06 19:01:28 |
Judging History
answer
//prepare for coding{{{
#include<bits/stdc++.h>
#define RELEASE
#ifndef RELEASE
#define FL
#define DB(...) fprintf(stderr,__VA_ARGS__)
#endif
typedef double DB;
const DB EPS=1e-8;
using namespace std;
const DB PAI=acos(-1);
inline int Read(){
char c=getchar();
int res=0;
bool b=0;
while(c>'9'||c<'0')
b=(c=='-'),c=getchar();
while(c>='0'&&c<='9')
res=(res<<3)+(res<<1)+c-'0',c=getchar();
return b?-res:res;
}
DB v,r,d;
//}}}
//solve{{{
const int B=5000;
inline DB Ftheta(DB theta){
DB sn=sin(theta);
return 2/(sqrt(sn*sn*4+12)-sn*2)/v;
}
inline DB Simpson(DB lim){
DB d=lim/B,res=Ftheta(0)+Ftheta(lim);
for(int i=1;i<B;i+=2)
res+=4*Ftheta(i*d);
for(int i=2;i<B;i+=2)
res+=2*Ftheta(i*d);
return res*d/3;
}
//}}}
//main{{{
inline void Solve(){
scanf("%lf%lf%lf",&v,&r,&d);
DB ll=PAI/2,rr=PAI,ans=PAI;
while(rr-ll>EPS){
DB mid=(ll+rr)/2,t=Simpson(mid)*r*2,alpha=atan2(sin(mid)*2*r,d-2*r-v*t+cos(mid)*2*r);
DB x=-cos(mid)*2*r+2*EPS*cos(alpha)-EPS,y=sin(mid)*2*r-2*EPS*sin(alpha);
if(x*x+y*y>4*r*r)
ans=rr=mid;
else
ans=ll=mid;
}
//printf("%lf\n",ans);
DB anst=Simpson(ans)*r*2;
DB x=cos(ans)*2*r+d-r*2-anst*v,y=sin(ans)*2*r;
//printf("%lf %lf %lf\n",anst,x,y);
printf("%.10lf\n",min(anst+sqrt(x*x+y*y)/v/2,d/v));
}
int main(){
#ifdef FL
freopen(".in","r",stdin);
freopen(".out","w",stdout);
#endif
int T=Read()+1;
while(--T)
Solve();
return 0;
}
//}}}
//《象》曰:泽灭木,大过。君子以独立不惧,遁世无
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 2ms
memory: 4288kb
input:
1 2.00 3 30.0
output:
8.3105799339
result:
ok found '8.3105799', expected '8.3105799', error '0.0000000'
Test #2:
score: -100
Wrong Answer
time: 86ms
memory: 4388kb
input:
100 8.66 6.05 71.59 8.44 8.67 82.55 4.23 9.68 24.74 6.63 9.95 7.45 1.2 7.27 97.55 2.08 2.94 31.61 8.48 5.56 20.16 7.35 5.55 52.48 1.35 8.15 51.23 3.88 9.02 72.37 7.03 1.97 57.03 1.78 4.44 93.44 5.5 6.11 52.94 1.09 9.64 69.01 1.01 3.54 16.78 1.63 4.55 44.06 5.82 5.51 5.92 5.62 2.93 76.01 4.75 4.43 73...
output:
4.4578371227 5.4883608090 5.8486997636 1.1236802413 43.2664800535 8.3054901382 2.3773584906 4.0154859695 37.9481481481 13.6188238200 4.1551236436 27.1789630071 6.1517775153 55.6333419215 16.6138613861 15.0939944612 1.0171821306 6.9496235027 8.1504302992 8.0094898462 1.7697642164 4.5411392405 10.4731...
result:
wrong answer 10th numbers differ - expected: '18.6520619', found: '13.6188238', error = '0.2698489'