Submission #429026

#	Time	Username	Problem	Language	Result	Execution time	Memory
429026		JeanBombeur	Aliens (IOI16_aliens)	C++17	4 / 100	1 ms	304 KiB

aliens

#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
#include "aliens.h"
using namespace std;

//    <|°_°|>

const long long INFINI = (1LL << 60);
const int MAX_POINTS = (100 * 1000);
const int LOG = (40);

struct droite {
    long long pente, cst;
    int prix;
};

pair <long long, long long> Points[MAX_POINTS];
pair <int, int> Aliens[MAX_POINTS];

pair <droite, long long> Deque[MAX_POINTS];
int deb = 0, fin = 0;

int nbPoints = 0;

bool Compare(pair <int, int> a, pair <int, int> b) {
    if (a.first != b.first)
        return a.first < b.first;
    return a.second > b.second;
}

long long Inter(droite a, droite b) {
    if (a.pente == b.pente)
        return INFINI;
    return (b.cst - a.cst) / (a.pente - b.pente);
}

long long Eval(droite a, long long absc) {
    return a.pente * absc + a.cst;
}

void Insert(droite a) {
    while (fin > deb && Eval(Deque[fin - 1].first, Deque[fin - 1].second) >= Eval(a, Deque[fin - 1].second))
        fin --;
    if (fin == deb)
        Deque[fin ++] = {a, 0};
    else
    {
        Deque[fin] = {a, Inter(Deque[fin - 1].first, a)};
        fin ++;
    }
    return;
}

pair <long long, int> GetMin(int abcs) {
    while (fin > deb + 1 && abcs >= Deque[deb + 1].second)
        deb ++;
    return {Eval(Deque[deb].first, abcs), Deque[deb].first.prix};
}

pair <long long, int> FindDP(long long lambda) {
    deb = 0, fin = 0;
    pair <long long, int> DP = {0LL, 0};
    Insert({- 2 * Points[0].first, Points[0].first * Points[0].first, 0});
    for (int i = 0; i < nbPoints; i ++)
    {
        DP = GetMin(Points[i].second);
        DP.first += Points[i].second * Points[i].second + lambda;
        DP.second ++;
        if (i + 1 < nbPoints)
        {
            long long d = max(0LL, Points[i].second - Points[i + 1].first);
            Insert({- 2 * Points[i + 1].first,
            Points[i + 1].first * Points[i + 1].first - d * d + DP.first, DP.second});
        }
    }
    return DP;
}

long long take_photos(int nbAliens, int tailleMax, int nbPhotos, vector<int> Lignes, vector<int> Colonnes) {
    
    tailleMax ++;
    
    for (int i = 0; i < nbAliens; i ++)
    {
        if (Lignes[i] > Colonnes[i])
            swap(Lignes[i], Colonnes[i]);
        Aliens[i] = {Lignes[i], Colonnes[i] + 1};
    }
    sort(Aliens, Aliens + nbAliens, Compare);
    
    for (int i = 0; i < nbAliens; i ++)
    {
        if (nbPoints == 0 || Aliens[i].second > (int)Points[nbPoints - 1].second)
            Points[nbPoints ++] = Aliens[i];
    }
    
    long long lambda = -1;
    long long ans = INFINI;
    for (long long i = (1LL << LOG); i > 0; i /= 2)
    {
        pair <long long, int> DP = FindDP(lambda + i);
        if (DP.second <= nbPhotos)
            ans = min(ans, DP.first - (lambda + i) * DP.second);
        else
            lambda += i;
    }
    
    return ans;
}

• Subtask 1: 4 / 4
• Subtask 2: 0 / 12
• Subtask 3: 0 / 9
• Subtask 4: 0 / 16
• Subtask 5: 0 / 19
• Subtask 6: 0 / 40