QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#102519 | #5433. Absolute Difference | lgvc# | WA | 2ms | 14020kb | C++23 | 3.6kb | 2023-05-03 14:18:14 | 2023-05-03 14:18:16 |
Judging History
answer
//这回只花了114514min就打完了。
//真好。记得多手造几组。
#include <bits/stdc++.h>
//#define int long long
#define pai 3.141592653589793238462643383279502884197169399375105820
#define MOD 1000000007
#define eps 0.1
inline int min(int a,int b) {return a<b?a:b;}
inline int max(int a,int b) {return a>b?a:b;}
#define ULL unsigned long long
#define LL long long
#define INF 0x3f3f3f3f
#define INF_LL 0x3f3f3f3f3f3f3f3f
static char buf[1000000],*paa=buf,*pd=buf;
static char buf2[1000000],*pp=buf2;
#define getchar() paa==pd&&(pd=(paa=buf)+fread(buf,1,1000000,stdin),paa==pd)?EOF:*paa++
inline void pc(char ch){
if(pp-buf2==1000000) fwrite(buf2,1,1000000,stdout),pp=buf2;
*pp++=ch;
}
inline void pcc(){
fwrite(buf2,1,pp-buf2,stdout);
pp=buf2;
}
inline int read(void) {
int w=1;
register int x(0);register char c(getchar());
while(c<'0'||c>'9') {if(c=='-') w=-1;c=getchar();}
while(c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
return w*x;
}
void write(int x) {
if(x<0) pc('-'),x=-x;
static int sta[20];
int top=0;
do {
sta[top++]=x%10,x/=10;
} while(x);
while(top) pc(sta[--top]+48);
}
void we(int x) {
write(x);
pc('\n');
}
inline bool cmp_xi(int a,int b) {return a<b;}
inline bool cmp_da(int a,int b) {return a>b;}
int N,M,k;
long double l1[100009],l2[100009],r1[100009],r2[100009],a[400009],su,su2[400009],qw2[400009],qw1[400009],ans;
bool v1[400009],v2[400009];
signed main(void) {
//freopen("m.in","r",stdin);
//freopen("m.out","w",stdout);
N=read();M=read();
a[++k]=-INF;
a[++k]=INF;
for(int i=1;i<=N;i++) {
l1[i]=read();
r1[i]=read();
a[++k]=l1[i];
a[++k]=r1[i]+eps;
}
for(int i=1;i<=M;i++) {
l2[i]=read();
r2[i]=read();
a[++k]=l2[i];
a[++k]=r2[i]+eps;
}
std::sort(a+1,a+k+1);
k=std::unique(a+1,a+k+1)-a-1;
for(int i=1;i<=N;i++) {
int l=1,r=k,md;
while(l<r) {
md=(l+r)>>1;
if(a[md]>=l1[i]) {
r=md;
} else {
l=md+1;
}
}
while(r1[i]>=a[l]) {
v1[l]=1;
l++;
}
}
for(int i=1;i<=M;i++) {
int l=1,r=k,md;
while(l<r) {
md=(l+r)>>1;
if(a[md]>=l2[i]) {
r=md;
} else {
l=md+1;
}
}
while(r2[i]>=a[l]) {
v2[l]=1;
l++;
}
}
for(int i=1;i<k;i++) {
if(v1[i]) su+=a[i+1]-a[i];
}
if(su<1) {
su=0;
for(int i=1;i<k;i++) {
su2[i]=su2[i-1];
qw2[i]=qw2[i-1];
if(v2[i]) su2[i]+=a[i+1]-a[i],qw2[i]+=(a[i+1]-a[i])*(a[i+1]+a[i]-eps)/2;
if(v1[i]) su+=a[i+1]-a[i];
}
for(int i=k-1;i>=1;i--) {
qw1[i]=qw1[i+1];
if(v2[i]) qw1[i]+=(a[i+1]-a[i])*(a[i+1]+a[i]-eps)/2;
}
for(int i=1;i<k;i++) {
if(v1[i]) {
if(su2[i-1]) ans+=((a[i]+a[i+1]-eps)/2-qw2[i-1]/su2[i-1])*(a[i+1]-a[i])/su*su2[i-1]/su2[k-1];
if(su2[k-1]!=su2[i]) ans+=(qw1[i+1]/(su2[k-1]-su2[i])-(a[i]+a[i+1]-eps)/2)*(a[i+1]-a[i])/su*(su2[k-1]-su2[i])/su2[k-1];
}
}
} else {
su=0;
for(int i=1;i<k;i++) {
su2[i]=su2[i-1];
qw2[i]=qw2[i-1];
if(v2[i]) su2[i]+=a[i+1]-a[i]-eps,qw2[i]+=(a[i+1]-a[i]-eps)*(a[i+1]+a[i]-eps)/2;
if(v1[i]) su+=a[i+1]-a[i]-eps;
}
for(int i=k-1;i>=1;i--) {
qw1[i]=qw1[i+1];
if(v2[i]) qw1[i]+=(a[i+1]-a[i]-eps)*(a[i+1]+a[i]-eps)/2;
}
for(int i=1;i<k;i++) {
if(v1[i]) {
if(v2[i]) ans+=(a[i+1]-a[i]-eps)/3*(a[i+1]-a[i]-eps)*(a[i+1]-a[i]-eps);
if(su2[i-1]) ans+=((a[i]+a[i+1]-eps)/2-qw2[i-1]/su2[i-1])*(a[i+1]-a[i]-eps)*su2[i-1];
if(su2[k-1]!=su2[i]) ans+=(qw1[i+1]/(su2[k-1]-su2[i])-(a[i]+a[i+1]-eps)/2)*(a[i+1]-a[i]-eps)*(su2[k-1]-su2[i]);
}
}
ans/=su;
ans/=su2[k-1];
}
printf("%.12Lf",ans);
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 2ms
memory: 14020kb
input:
1 1 0 1 0 1
output:
0.333333333333
result:
ok found '0.333333333', expected '0.333333333', error '0.000000000'
Test #2:
score: -100
Wrong Answer
time: 2ms
memory: 13844kb
input:
1 1 0 1 1 1
output:
-nan
result:
wrong output format Expected double, but "-nan" found