Submission #32186

# Submission time Handle Problem Language Result Execution time Memory
32186 2017-10-03T11:22:28 Z dongwon0427 Aliens (IOI16_aliens) C++11
4 / 100
0 ms 10224 KB
#include "aliens.h"
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<int,ld> pil;
struct _tuple {
    int x,y;
};
_tuple A[100005],B[100005];
int n;
ld dp[100005];
int gaesu[100005];
struct line {
    ld tan,y;
    int idx;
};
inline ld inter(line a,line b) {
    return (b.y-a.y)/(a.tan-b.tan);
}
line L[100005];
pil CHT(ld c) {
    gaesu[0]=0;
    dp[1] = (ld)(A[1].y-A[1].x+1)*(ld)(A[1].y-A[1].x+1)-c;
    gaesu[1]=1;
    int sz = 0;
    int p=0;
    for(int i=2;i<=n;i++) {
        line ins;
        ins.idx = i-1;
        ins.y = dp[i-1] - (ld)max(0,-A[i].x+A[i-1].y+1)*(ld)max(0,-A[i].x+A[i-1].y+1) + (ld)(A[i].x)*(ld)(A[i].x) - (ld)2*(ld)(A[i].x);
        ins.tan = (ld)(-2)*(ld)(A[i].x);
        while(sz>=2 && inter(L[sz-2],L[sz-1]) >= inter(L[sz-1],ins) ) {
            if(p==sz) p--;
            sz--;
        }
        L[sz++] = ins;

        while(p<sz-1 && inter(L[p],L[p+1]) <= (ld)A[i].y) p++;

        dp[i] = (ld)(1) - c + (ld)(2)*(ld)(A[i].y) + (ld)(A[i].y)*(ld)(A[i].y) + L[p].tan*(ld)(A[i].y) + L[p].y;
        int minidx = L[p].idx;
        gaesu[i] = gaesu[minidx] + 1;
    }

    return pil(gaesu[n],dp[n]);
    /*gaesu[0]=0;
    for(int i=1;i<=n;i++) {
        dp[i] = (ld)(A[i].y-A[1].x+1)*(ld)(A[i].y-A[1].x+1)-c;
        int minidx = 0;
        for(int j=1;j<i;j++) {
            if( dp[i] > dp[j] + (ld)(A[i].y-A[j+1].x+1)*(ld)(A[i].y-A[j+1].x+1) - (ld)max(0,-A[j+1].x+A[j].y+1)*(ld)max(0,-A[j+1].x+A[j].y+1) - c ) {
                minidx = j;
                dp[i] = dp[j] + (ld)(A[i].y-A[j+1].x+1)*(ld)(A[i].y-A[j+1].x+1) - (ld)max(0,-A[j+1].x+A[j].y+1)*(ld)max(0,-A[j+1].x+A[j].y+1) - c;
            }
        }
        gaesu[i] = gaesu[minidx]+1;
    }
    return pil(gaesu[n],dp[n]);*/
}
ll take_photos(int N, int m, int k, std::vector<int> r, std::vector<int> c) {
    for(int i=0;i<N;i++) {
        B[i].x=r[i];
        B[i].y=c[i];
        if(B[i].x > B[i].y) swap(B[i].x,B[i].y);
    }
    sort(B,B+N,[](_tuple a,_tuple b){return a.x>b.x;});
    for(int i=0;i<N;i++) {
        while(n!=0 && A[n].y <= B[i].y) {
            n--;
        }
        A[++n] = B[i];
    }
    reverse(A+1,A+1+n);
    if(k>n) k = n;
    ll bs = -(ll)m*(ll)m-1ll, be = 0;
    int lft=1,rht=n;
    ll L=(ll)m*(ll)m,R=(ll)m*(ll)m;
    while(bs<=be) {
        ll bm = (bs+be)/2;
        //printf("%d\n",bm);
        ld gi = (ld)bm + (ld)0.5;
        pil tmp = CHT(gi);
        if(tmp.first==k) {
            return (long long)(tmp.second + gi*(ld)(k));
        }
        if(tmp.first > k) {
            be=bm-1;
            if(rht >= tmp.first) {
                rht = tmp.first;
                R = (ll)(tmp.second + gi*(ld)(tmp.first));
            }
        } else {
            bs=bm+1;
            if(lft <= tmp.first) {
                lft = tmp.first;
                L = (ll)(tmp.second + gi*(ld)(tmp.first));
            }
        }
    }
    return (L*(ll)(rht-k) + R*(ll)(k-lft)) / (ll)(-lft+rht);
}
# Verdict Execution time Memory Grader output
1 Correct 0 ms 10224 KB Correct answer: answer = 4
2 Correct 0 ms 10224 KB Correct answer: answer = 4
3 Correct 0 ms 10224 KB Correct answer: answer = 4
4 Correct 0 ms 10224 KB Correct answer: answer = 12
5 Correct 0 ms 10224 KB Correct answer: answer = 52
6 Correct 0 ms 10224 KB Correct answer: answer = 210
7 Correct 0 ms 10224 KB Correct answer: answer = 88
8 Correct 0 ms 10224 KB Correct answer: answer = 7696
9 Correct 0 ms 10224 KB Correct answer: answer = 1
10 Correct 0 ms 10224 KB Correct answer: answer = 2374
11 Correct 0 ms 10224 KB Correct answer: answer = 9502
12 Correct 0 ms 10224 KB Correct answer: answer = 49
13 Correct 0 ms 10224 KB Correct answer: answer = 151
14 Correct 0 ms 10224 KB Correct answer: answer = 7550
15 Correct 0 ms 10224 KB Correct answer: answer = 7220
16 Correct 0 ms 10224 KB Correct answer: answer = 7550
17 Correct 0 ms 10224 KB Correct answer: answer = 10000
18 Correct 0 ms 10224 KB Correct answer: answer = 10000
19 Correct 0 ms 10224 KB Correct answer: answer = 624
20 Correct 0 ms 10224 KB Correct answer: answer = 10000
# Verdict Execution time Memory Grader output
1 Correct 0 ms 10224 KB Correct answer: answer = 1
2 Correct 0 ms 10224 KB Correct answer: answer = 4
3 Correct 0 ms 10224 KB Correct answer: answer = 1
4 Correct 0 ms 10224 KB Correct answer: answer = 5
5 Incorrect 0 ms 10224 KB Wrong answer: output = 81, expected = 41
6 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 0 ms 10224 KB Correct answer: answer = 4
2 Correct 0 ms 10224 KB Correct answer: answer = 4
3 Correct 0 ms 10224 KB Correct answer: answer = 4
4 Correct 0 ms 10224 KB Correct answer: answer = 12
5 Correct 0 ms 10224 KB Correct answer: answer = 52
6 Correct 0 ms 10224 KB Correct answer: answer = 210
7 Correct 0 ms 10224 KB Correct answer: answer = 88
8 Correct 0 ms 10224 KB Correct answer: answer = 7696
9 Correct 0 ms 10224 KB Correct answer: answer = 1
10 Correct 0 ms 10224 KB Correct answer: answer = 2374
11 Correct 0 ms 10224 KB Correct answer: answer = 9502
12 Correct 0 ms 10224 KB Correct answer: answer = 49
13 Correct 0 ms 10224 KB Correct answer: answer = 151
14 Correct 0 ms 10224 KB Correct answer: answer = 7550
15 Correct 0 ms 10224 KB Correct answer: answer = 7220
16 Correct 0 ms 10224 KB Correct answer: answer = 7550
17 Correct 0 ms 10224 KB Correct answer: answer = 10000
18 Correct 0 ms 10224 KB Correct answer: answer = 10000
19 Correct 0 ms 10224 KB Correct answer: answer = 624
20 Correct 0 ms 10224 KB Correct answer: answer = 10000
21 Correct 0 ms 10224 KB Correct answer: answer = 1
22 Correct 0 ms 10224 KB Correct answer: answer = 4
23 Correct 0 ms 10224 KB Correct answer: answer = 1
24 Correct 0 ms 10224 KB Correct answer: answer = 5
25 Incorrect 0 ms 10224 KB Wrong answer: output = 81, expected = 41
26 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 0 ms 10224 KB Correct answer: answer = 4
2 Correct 0 ms 10224 KB Correct answer: answer = 4
3 Correct 0 ms 10224 KB Correct answer: answer = 4
4 Correct 0 ms 10224 KB Correct answer: answer = 12
5 Correct 0 ms 10224 KB Correct answer: answer = 52
6 Correct 0 ms 10224 KB Correct answer: answer = 210
7 Correct 0 ms 10224 KB Correct answer: answer = 88
8 Correct 0 ms 10224 KB Correct answer: answer = 7696
9 Correct 0 ms 10224 KB Correct answer: answer = 1
10 Correct 0 ms 10224 KB Correct answer: answer = 2374
11 Correct 0 ms 10224 KB Correct answer: answer = 9502
12 Correct 0 ms 10224 KB Correct answer: answer = 49
13 Correct 0 ms 10224 KB Correct answer: answer = 151
14 Correct 0 ms 10224 KB Correct answer: answer = 7550
15 Correct 0 ms 10224 KB Correct answer: answer = 7220
16 Correct 0 ms 10224 KB Correct answer: answer = 7550
17 Correct 0 ms 10224 KB Correct answer: answer = 10000
18 Correct 0 ms 10224 KB Correct answer: answer = 10000
19 Correct 0 ms 10224 KB Correct answer: answer = 624
20 Correct 0 ms 10224 KB Correct answer: answer = 10000
21 Correct 0 ms 10224 KB Correct answer: answer = 1
22 Correct 0 ms 10224 KB Correct answer: answer = 4
23 Correct 0 ms 10224 KB Correct answer: answer = 1
24 Correct 0 ms 10224 KB Correct answer: answer = 5
25 Incorrect 0 ms 10224 KB Wrong answer: output = 81, expected = 41
26 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 0 ms 10224 KB Correct answer: answer = 4
2 Correct 0 ms 10224 KB Correct answer: answer = 4
3 Correct 0 ms 10224 KB Correct answer: answer = 4
4 Correct 0 ms 10224 KB Correct answer: answer = 12
5 Correct 0 ms 10224 KB Correct answer: answer = 52
6 Correct 0 ms 10224 KB Correct answer: answer = 210
7 Correct 0 ms 10224 KB Correct answer: answer = 88
8 Correct 0 ms 10224 KB Correct answer: answer = 7696
9 Correct 0 ms 10224 KB Correct answer: answer = 1
10 Correct 0 ms 10224 KB Correct answer: answer = 2374
11 Correct 0 ms 10224 KB Correct answer: answer = 9502
12 Correct 0 ms 10224 KB Correct answer: answer = 49
13 Correct 0 ms 10224 KB Correct answer: answer = 151
14 Correct 0 ms 10224 KB Correct answer: answer = 7550
15 Correct 0 ms 10224 KB Correct answer: answer = 7220
16 Correct 0 ms 10224 KB Correct answer: answer = 7550
17 Correct 0 ms 10224 KB Correct answer: answer = 10000
18 Correct 0 ms 10224 KB Correct answer: answer = 10000
19 Correct 0 ms 10224 KB Correct answer: answer = 624
20 Correct 0 ms 10224 KB Correct answer: answer = 10000
21 Correct 0 ms 10224 KB Correct answer: answer = 1
22 Correct 0 ms 10224 KB Correct answer: answer = 4
23 Correct 0 ms 10224 KB Correct answer: answer = 1
24 Correct 0 ms 10224 KB Correct answer: answer = 5
25 Incorrect 0 ms 10224 KB Wrong answer: output = 81, expected = 41
26 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 0 ms 10224 KB Correct answer: answer = 4
2 Correct 0 ms 10224 KB Correct answer: answer = 4
3 Correct 0 ms 10224 KB Correct answer: answer = 4
4 Correct 0 ms 10224 KB Correct answer: answer = 12
5 Correct 0 ms 10224 KB Correct answer: answer = 52
6 Correct 0 ms 10224 KB Correct answer: answer = 210
7 Correct 0 ms 10224 KB Correct answer: answer = 88
8 Correct 0 ms 10224 KB Correct answer: answer = 7696
9 Correct 0 ms 10224 KB Correct answer: answer = 1
10 Correct 0 ms 10224 KB Correct answer: answer = 2374
11 Correct 0 ms 10224 KB Correct answer: answer = 9502
12 Correct 0 ms 10224 KB Correct answer: answer = 49
13 Correct 0 ms 10224 KB Correct answer: answer = 151
14 Correct 0 ms 10224 KB Correct answer: answer = 7550
15 Correct 0 ms 10224 KB Correct answer: answer = 7220
16 Correct 0 ms 10224 KB Correct answer: answer = 7550
17 Correct 0 ms 10224 KB Correct answer: answer = 10000
18 Correct 0 ms 10224 KB Correct answer: answer = 10000
19 Correct 0 ms 10224 KB Correct answer: answer = 624
20 Correct 0 ms 10224 KB Correct answer: answer = 10000
21 Correct 0 ms 10224 KB Correct answer: answer = 1
22 Correct 0 ms 10224 KB Correct answer: answer = 4
23 Correct 0 ms 10224 KB Correct answer: answer = 1
24 Correct 0 ms 10224 KB Correct answer: answer = 5
25 Incorrect 0 ms 10224 KB Wrong answer: output = 81, expected = 41
26 Halted 0 ms 0 KB -