QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#559376#5433. Absolute DifferenceWuyanruWA 1ms7808kbC++142.5kb2024-09-11 21:46:072024-09-11 21:46:07

Judging History

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

  • [2024-09-11 21:46:07]
  • 评测
  • 测评结果:WA
  • 用时:1ms
  • 内存:7808kb
  • [2024-09-11 21:46:07]
  • 提交

answer

#include<bits/stdc++.h>
#define inf 0x3f3f3f3f3f3f3f3fll
#define debug(x) cerr<<#x<<"="<<x<<endl
using namespace std;
using ll=long long;
using ld=long double;
using pli=pair<ll,int>;
using pi=pair<int,int>;
template<typename A>
using vc=vector<A>;
using pd=pair<ld,ld>;
inline int read()
{
    int s=0,w=1;char ch;
    while((ch=getchar())>'9'||ch<'0') if(ch=='-') w=-1;
    while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
    return s*w;
}
inline ll lread()
{
    ll s=0,w=1;char ch;
    while((ch=getchar())>'9'||ch<'0') if(ch=='-') w=-1;
    while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
    return s*w;
}
const ld L=9000000000000000ll;
ld num[400005];
int numa[400005];
int numb[400005];
int len;
inline ld get(int *a,int *b)
{
    ld ans=0,lena=0,suma=0,lenb=0,sumb=0;
    for(int i=1;i<=len;i++)
    {
        ld l=num[i],r=num[i+1];
        if(a[i]) lena+=r-l,suma+=(r-l)*(l+r)/2;
        if(b[i]) lenb+=r-l,sumb+=(r-l)*(l+r)/2;
    }
    suma/=lena,sumb/=lenb;
    // printf("%Lf %Lf %Lf %Lf\n",lena,suma,lenb,sumb);

    ld P=0;
    for(int i=1;i<=len;i++)
    {
        ld l=num[i],r=num[i+1];
        ld Pb=(r-l)/lenb,Pa=(r-l)/lena;
        if(a[i])
        {
            ans+=P*Pa*(l+r)/2;
            // printf("i=%d ans=%Lf\n",i,ans);
            if(b[i])
            {
                ld t=r*r*r/3-l*l*l/3-l*r*r/2+l*l*l/2;
                ans+=t/lena/lenb;
            }
            // printf("i=%d ans=%Lf\n",i,ans);
        }
        if(b[i]) P+=Pb;
    }

    return 2*ans-suma;
}
int main()
{
    int n=read(),m=read();
    vc<pd>Sa,Sb;
    for(int i=1;i<=n;i++){ ld l=num[++len]=lread()*L,r=num[++len]=lread()*L+1;Sa.push_back(pd(l,r));}
    for(int i=1;i<=m;i++){ ll l=num[++len]=lread()*L,r=num[++len]=lread()*L+1;Sb.push_back(pd(l,r));}
    sort(num+1,num+len+1),len=unique(num+1,num+len+1)-num-1;
    for(auto i:Sa)
    {
        int l=lower_bound(num+1,num+len+1,i.first)-num;
        int r=lower_bound(num+1,num+len+1,i.second)-num;
        for(int j=l;j<r;j++) numa[j]=1;
    }
    for(auto i:Sb)
    {
        int l=lower_bound(num+1,num+len+1,i.first)-num;
        int r=lower_bound(num+1,num+len+1,i.second)-num;
        for(int j=l;j<r;j++) numb[j]=1;
    }
    // for(int i=1;i<=len;i++) printf("%lld ",num[i]);;putchar('\n');
    // for(int i=1;i<=len;i++) printf("%d ",numa[i]);;putchar('\n');
    // for(int i=1;i<=len;i++) printf("%d ",numb[i]);;putchar('\n');
    printf("%.20Lf\n",(get(numa,numb)+get(numb,numa))/L);
    return 0;
}

詳細信息

Test #1:

score: 100
Accepted
time: 1ms
memory: 7808kb

input:

1 1
0 1
0 1

output:

0.33333333333333337034

result:

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

Test #2:

score: 0
Accepted
time: 1ms
memory: 5984kb

input:

1 1
0 1
1 1

output:

0.49999999999999977777

result:

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

Test #3:

score: -100
Wrong Answer
time: 1ms
memory: 7808kb

input:

1 1
-1000000000 1000000000
-1000000000 1000000000

output:

-nan

result:

wrong output format Expected double, but "-nan" found