Submission #1094707

#	Time	Username	Problem	Language	Result	Execution time	Memory
1094707		buzdi	Graph (BOI20_graph)	C++17	0 / 100	2 ms	2788 KiB

Graph

#include <iostream>
#include <algorithm>
#include <cassert>
#include <cmath>
#include <iomanip>
#include <vector>
#include <set>
#define ll long long

using namespace std;

const int NMAX = 1e5;
const int INF = 1e9;
const long double EPS = 1e-6;

struct Value {
    int a, b;
}values[NMAX + 1];

int n, m;
vector<pair<int, int>> G[NMAX + 1];
bool visited[NMAX + 1];
bool ok;
long double unique_solution;

int GCD(int a, int b) {
    while (b) {
        int r = a % b;
        a = b;
        b = r;
    }
    return a;
}

Value make_value(int a, int b) {
    Value answer;
    answer.a = a;
    answer.b = b;
    return answer;
}

void Simplify(pair<int, int> &a) {
    int gcd = GCD(abs(a.first), abs(a.second));
    a.first /= gcd;
    a.second /= gcd;
}

bool GetSolution(Value x, Value y, int which) {
    if (x.a + y.a == 0) {
        return x.b + y.b == which;
    }

    pair<int, int> solution_x = { which - x.b - y.b, x.a + y.a };
    if ((solution_x.first <= 0 && solution_x.second <= 0) || (solution_x.first >= 0 && solution_x.second <= 0)) {
        solution_x.first = -solution_x.first;
        solution_x.second = -solution_x.second;
    }
    Simplify(solution_x);
    long double real_solution_x = (long double) solution_x.first / (long double) solution_x.second;
    if(unique_solution == INF) {
        unique_solution = real_solution_x;
        return true;
    }
    return abs(unique_solution - real_solution_x) <= EPS;
}

void DFS(int node, Value curent_value) {
    values[node] = curent_value;
    visited[node] = 1;
    for (auto [next_node, which] : G[node]) {
        if (!visited[next_node]) {
            DFS(next_node, make_value(-curent_value.a, which - curent_value.b));
        }
        else {
            if(!GetSolution(values[node], values[next_node], which)) {
                ok = false;
            }
        }
    }
}

long double GetValue(long double x) {
    long double answer = 0;
    for (int i = 1; i <= n; i++) {
        answer += abs(values[i].a * x + values[i].b);
    }
    return answer;
}

signed main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);

    cin >> n >> m;
    for (int i = 1; i <= m; i++) {
        int a, b, c;
        cin >> a >> b >> c;
        G[a].push_back({ b, c });
        G[b].push_back({ a, c });
    }

    ok = true;
    unique_solution = INF;
    for (int i = 1; i <= n; i++) {
        if (!visited[i]) {
            DFS(i, make_value(1, 0));
        }
    }

    if (ok) {
        cout << "YES" << '\n';
        if (unique_solution != INF) {
            for (int i = 1; i <= n; i++) {
                cout << fixed << setprecision(6) << values[i].a * unique_solution + values[i].b << ' ';
            }
        }
        else {
            long double left = -INF, right = INF;
            int iterations = 400;
            while (iterations--) {
                long double mid1 = left + (right - left) / 3.0;
                long double mid2 = right - (right - left) / 3.0;

                long double value1 = GetValue(mid1);
                long double value2 = GetValue(mid2);
                if (value1 > value2) {
                    left = mid1;
                }
                else {
                    right = mid2;
                }
            }

            for (int i = 1; i <= n; i++) {
                cout << fixed << setprecision(6) << values[i].a * left + values[i].b << ' ';
            }
        }
    }
    else {
//        assert(false);
        cout << "NO" << '\n';
    }

    return 0;
}

• Subtask 1: 0 / 5
• Subtask 2: 0 / 12
• Subtask 3: 0 / 17
• Subtask 4: 0 / 24
• Subtask 5: 0 / 42