QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#102486#5433. Absolute Differencelgvc#WA 2ms5912kbC++232.8kb2023-05-03 13:51:032023-05-03 13:51:04

Judging History

你现在查看的是最新测评结果

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-05-03 13:51:04]
  • 评测
  • 测评结果:WA
  • 用时:2ms
  • 内存:5912kb
  • [2023-05-03 13:51:03]
  • 提交

answer

//这回只花了114514min就打完了。
//真好。记得多手造几组。
#include <bits/stdc++.h>
//#define int long long
#define pai 3.141592653589793238462643383279502884197169399375105820
#define MOD 1000000007
#define eps 0.0000000001
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]=0;
	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++) {
		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])/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])/2;
	}
	for(int i=1;i<k;i++) {
		if(v1[i]) {
			if(v2[i]) ans+=(a[i+1]-a[i])/3*(a[i+1]-a[i])/su;
			if(su2[i-1]) ans+=((a[i]+a[i+1])/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])/2)*(a[i+1]-a[i])/su*(su2[k-1]-su2[i])/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: 5800kb

input:

1 1
0 1
0 1

output:

0.333333333367

result:

ok found '0.333333333', expected '0.333333333', error '0.000000000'

Test #2:

score: 0
Accepted
time: 2ms
memory: 5800kb

input:

1 1
0 1
1 1

output:

0.500000000000

result:

ok found '0.500000000', expected '0.500000000', error '0.000000000'

Test #3:

score: -100
Wrong Answer
time: 2ms
memory: 5912kb

input:

1 1
-1000000000 1000000000
-1000000000 1000000000

output:

833333333.333333333372

result:

wrong answer 1st numbers differ - expected: '666666666.6666666', found: '833333333.3333334', error = '0.2500000'