Submission #804892

#	Time	Username	Problem	Language	Result	Execution time	Memory
804892		Boas	Detecting Molecules (IOI16_molecules)	C++17	46 / 100	53 ms	948 KiB

molecules

#include "molecules.h"

using namespace std;
#include <bits/stdc++.h>

#define MAX_N 200001
#define ALL(x) x.begin(), x.end()
#define loop(x) for (int i = 0; i < (x); i++)
typedef pair<int, int> ii;
typedef vector<ii> vii;
typedef vector<int> vi;
typedef set<int> si;

typedef int64_t in64;

int diff(in64 i, in64 l, in64 u)
{
    if (i < l)
        return l - i;
    if (i > u)
        return i - u;
    return 0;
}

vi find_subset(int l, int u, vi w)
{
    in64 alles = accumulate(ALL(w), 0);
    if (alles < l)
        return {};
    if (alles <= u)
    {
        vi ret;
        loop(w.size()) ret.push_back(i);
        return ret;
    }
    vii sortedMolecules; // first weight, second real index
    loop(w.size())
    {
        if (w[i] > u)
            continue;
        sortedMolecules.push_back({w[i], i});
    }
    sort(ALL(sortedMolecules));
    si used;
    bitset<MAX_N> usin;
    in64 totalW = 0;
    while (totalW < l)
    {
        bool added = false;
        for (int i = sortedMolecules.size() - 1; i >= 0; i--)
        {
            if (usin[i])
                continue;
            in64 W = sortedMolecules[i].first; // weight, real index
            if (totalW + W > u)
            {
                continue;
            }
            added = true;
            used.insert(i);
            usin[i] = true;
            totalW += W;
        }
        if (!added)
            break;
    }
    if (totalW >= l)
    {
        vi indexesOfUsedMols;
        for (int i : used)
        {
            indexesOfUsedMols.push_back(sortedMolecules[i].second);
        }
        return indexesOfUsedMols;
    }
    if (sortedMolecules.size() < 1)
        return {};
    in64 smallestW = sortedMolecules[0].first;
    used.insert(0);
    totalW += smallestW;
    usin[0] = true;
    si oldUsed(used);            // sorted indexes
    bitset<MAX_N> oldUsin(usin); // sorted indexes
    in64 oldW = totalW;
    for (int i1 = 0; i1 < sortedMolecules.size(); i1++)
    {
        used = oldUsed;
        usin = oldUsin;
        if (usin[i1])
            continue;
        ii p1 = sortedMolecules[i1];
        for (int i2 : oldUsed)
        {
            used = oldUsed;
            usin = oldUsin;
            ii p2 = sortedMolecules[i2];
            used = oldUsed;
            used.erase(i2);
            usin[i2] = false;
            used.insert(i1);
            usin[i1] = true;
            totalW = oldW - p2.first + p1.first;
            if (diff(totalW, l, u) > diff(oldW, l, u))
                continue;
            while (totalW < l)
            {
                bool added = false;
                for (int i = sortedMolecules.size() - 1; i >= 0; i--)
                {
                    if (usin[i])
                        continue;
                    int W = sortedMolecules[i].first;
                    if (totalW + W > u)
                    {
                        continue;
                    }
                    added = true;
                    used.insert(i);
                    usin[i] = true;
                    totalW += W;
                }
                if (!added)
                    break;
            }
            if (totalW >= l && totalW <= u)
            {
                vi indexesOfUsedMols;
                for (int i : used)
                {
                    indexesOfUsedMols.push_back(sortedMolecules[i].second);
                }
                return indexesOfUsedMols;
            }
            if (diff(totalW, l, u) < diff(oldW, l, u))
            {
                oldW = totalW;
                oldUsed = used;
                oldUsin = usin;
                i1 = 0;
                break;
            }
        }
    }
    return {};
}

Compilation message (stderr)

molecules.cpp: In function 'vi find_subset(int, int, vi)':
molecules.cpp:8:35: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
    8 | #define loop(x) for (int i = 0; i < (x); i++)
      |                                 ~~^~~~~
molecules.cpp:33:9: note: in expansion of macro 'loop'
   33 |         loop(w.size()) ret.push_back(i);
      |         ^~~~
molecules.cpp:8:35: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
    8 | #define loop(x) for (int i = 0; i < (x); i++)
      |                                 ~~^~~~~
molecules.cpp:37:5: note: in expansion of macro 'loop'
   37 |     loop(w.size())
      |     ^~~~
molecules.cpp:85:25: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<int, int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   85 |     for (int i1 = 0; i1 < sortedMolecules.size(); i1++)
      |                      ~~~^~~~~~~~~~~~~~~~~~~~~~~~

• Subtask 1: 9 / 9
• Subtask 2: 10 / 10
• Subtask 3: 12 / 12
• Subtask 4: 15 / 15
• Subtask 5: 0 / 23
• Subtask 6: 0 / 31