Submission #1030535

#	Time	Username	Problem	Language	Result	Execution time	Memory
1030535		12345678	Graph (BOI20_graph)	C++17	0 / 100	1 ms	2652 KiB

Graph

#include <bits/stdc++.h>

using namespace std;

const int nx=1e5+5;

int n, m, u, v, w, a[nx], b[nx], f, vs[nx];
double res[nx];
vector<tuple<int, int, int>> d[nx];
vector<int> comp;
set<double> s;

void dfs(int u, int p)
{
    vs[u]=1;
    comp.push_back(u);
    for (auto [v, w, idx]:d[u])
    {
        if (idx==p) continue;
        if (vs[v])
        {
            int sma=a[u]+a[v], smb=b[u]+b[v];
            if (sma==0&&b[u]+b[v]==w) continue;
            else if (sma==0&&b[u]+b[v]!=w)
            {
                f=1;
                break;
            }
            s.insert(((double)(w-smb))/((double)sma));
        }
        else
        {
            a[v]=-a[u];
            b[v]=w-b[u];
            dfs(v, idx);
        }
    }
}
int main()
{
    cin.tie(NULL)->sync_with_stdio(false);
    cin>>n>>m;
    for (int i=1; i<=m; i++) cin>>u>>v>>w, d[u].push_back({v, w, i}), d[v].push_back({u, w, i});
    for (int i=1; i<=n; i++)
    {
        if (!vs[i])
        {
            comp.clear();
            s.clear();
            a[i]=1;
            dfs(i, -1);
            if (s.size()>1) f=1;
            else if (s.size()==1)
            {
                double vl=*s.begin();
                for (auto x:comp) res[x]=a[x]*vl+b[x];
            }
            else
            {
                vector<double> v;
                for (auto x:comp) if (a[x]!=0) v.push_back((double)(-b[x])/(double)a[x]);
                sort(v.begin(), v.end());
                cout<<"v ";
                for (auto x:v) cout<<x<<' ';
                cout<<'\n';
                int l=0, r=(int)v.size()-1;
                while (l<r)
                {
                    int md=(l+r+1)/2;
                    //cout<<"md "<<md<<'\n';
                    double sm1=0, sm2=0;
                    for (auto x:comp) sm1+=abs((double)a[x]*v[md]+(double)b[x]), sm2+=abs((double)a[x]*v[md-1]+(double)b[x]);
                    //cout<<"sm1 "<<sm1<<' '<<"sm2 "<<sm2<<'\n';
                    if (sm1<sm2) l=md;
                    else r=md-1;
                }
                //cout<<"debug "<<l<<' '<<v[l]<<'\n';
                double t=0;
                if (v.empty()) t=0;
                else t=v[l];
                //cout<<"mn "<<mn.second<<'\n';
                //if (t!=mn.second) cout<<1/0;
                for (auto x:comp) res[x]=t*a[x]+b[x];
            }
        }
    }
    if (f) return cout<<"NO\n", 0;
    cout<<"YES\n";
    for (int i=1; i<=n; i++) cout<<setprecision(8)<<res[i]<<' ';
}

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