QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#417235#5433. Absolute DifferenceocharinWA 0ms4064kbC++232.4kb2024-05-22 16:34:042024-05-22 16:34:04

Judging History

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

  • [2024-05-22 16:34:04]
  • 评测
  • 测评结果:WA
  • 用时:0ms
  • 内存:4064kb
  • [2024-05-22 16:34:04]
  • 提交

answer

/*

                                _/                                      _/
       _/_/        _/_/_/      _/_/_/        _/_/_/      _/  _/_/               _/_/_/
    _/    _/    _/            _/    _/    _/    _/      _/_/          _/       _/    _/
   _/    _/    _/            _/    _/    _/    _/      _/            _/       _/    _/
    _/_/        _/_/_/      _/    _/      _/_/_/      _/            _/       _/    _/

*/
#include<bits/stdc++.h>
#define int long long

using namespace std;
const int inf=1e18;

signed main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    int n,m;
    cin>>n>>m;
    vector<array<int,2>>a(n),b(m);
    for(auto &[x,y]:a) cin>>x>>y;
    for(auto &[x,y]:b) cin>>x>>y;
    sort(a.begin(),a.end());
    sort(b.begin(),b.end());
    double res=0;
    auto cal1=[&]()->void{
        int j=0;
        double c0=0,c1=0,c2=0;
        for(int i=0;i<a.size();++i){
            while(j<b.size() && b[j][1]<=a[i][0]){
                if(b[j][0]==b[j][1]) c0+=b[j][0];
                else{
                    c1+=b[j][1]-b[j][0];
                    c2+=b[j][1]*b[j][1]-b[j][0]*b[j][0];
                }
                j++;
            }
            if(a[i][0]==a[i][1]) res+=a[i][0]*c1-0.5l*c2;
            else{
                res+=(a[i][1]*a[i][1]-a[i][0]*a[i][0])*1l/2;
                res+=(a[i][1]*a[i][1]-a[i][0]*a[i][0])*1l/2*c1;
                res-=(a[i][1]-a[i][0])*(c0+c2*0.5l);
            }
        }
    };
    cal1();
    swap(a,b);
    cal1();
    swap(a,b);
    vector<array<int,2>>c;
    for(int i=0;i<n;++i){
        int x=lower_bound(b.begin(),b.end(),a[i])-b.begin();
        if(x){
            if(b[x-1][1]>a[i][0]) c.push_back({i,x-1});
        }
        while(x<m && b[x][0]<a[i][1]){
            c.push_back({i,x});
            x++;
        }
    }
    auto cal2=[&](array<int,2>a,array<int,2>b)->void{
        if(a>b) swap(a,b);
        long double l1=a[0],r1=a[1],l2=b[0],r2=b[1];
        if(l1<l2) res+=(l2-l1)*0.5l*(r2*r2-l2*l2)-(l2*l2-l1*l1)*0.5l*(r2-l2);
        if(r2<r1) res+=(r2-l2)*0.5l*(r1*r1-r2*r2)-(r1-r2)*0.5l*(r2*r2-l2*l2);
        if(r2<=r1) res+=r2*r2*r2*1l/3-l2*l2*l2*1l/3-(l2+r2)*1l/2*(r2*r2-l2*l2)+(r2-l2)*1l/2*(r2*r2+l2*l2);
        if(r2>r1) res+=r1*r1*r1*1l/3-l2*l2*l2*1l/3-(l2+r2)*1l/2*(r1*r1-l2*l2)+(r1-l2)*1l/2*(r2*r2+l2*l2);
    };
    for(auto [x,y]:c) cal2(a[x],b[y]);
    cout<<fixed<<setprecision(15)<<res<<"\n";
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 0ms
memory: 4064kb

input:

1 1
0 1
0 1

output:

0.333333333333333

result:

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

Test #2:

score: 0
Accepted
time: 0ms
memory: 3836kb

input:

1 1
0 1
1 1

output:

0.500000000000000

result:

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

Test #3:

score: -100
Wrong Answer
time: 0ms
memory: 3836kb

input:

1 1
-1000000000 1000000000
-1000000000 1000000000

output:

2666666666666666518848208896.000000000000000

result:

wrong answer 1st numbers differ - expected: '666666666.6666666', found: '2666666666666666518848208896.0000000', error = '4000000000000000000.0000000'