【学习笔记】[COCI2018-2019#1] Teoretičar

首先,可以发现 C C C等于所有点度数的最大值 ,我们能用到的颜色数目为 2 x ≥ C 2^x\ge C 2x≥C。

考虑分治,将边集划分为 E = E 1 + E 2 E=E_1+E_2 E=E1+E2,使得 E 1 , E 2 E_1,E_2 E1,E2中点度数的最大值都不超过 2 x − 1 2^{x-1} 2x−1。这样,我们将 [ 1 , 2 x − 1 ] [1,2^{x-1}] [1,2x−1]中的颜色分配给 E 1 E_1 E1, [ 2 x − 1 + 1 , 2 x ] [2^{x-1}+1,2^x] [2x−1+1,2x]中的颜色分配给 E 2 E_2 E2,就可以递归到子问题。显然当 C = 1 C=1 C=1时,可以将所有边染成同一个颜色。

显然,我们可以用欧拉回路解决这个问题。但是这个做法常数比较大。考虑 CF1499G Graph Coloring 的做法,维护若干个环和若干条端点互不相同的链(因为是二分图,所以环的长度一定是偶数),这可以用带权并查集维护。注意并查集的维护对象是每条边,合并两条边的时候限制等价于这两条边所属集合不同。

复杂度 O ( m log ⁡ m ) O(m\log m) O(mlogm)。

cpp 复制代码
#include<bits/stdc++.h>
#define ll long long
#define fi first
#define se second
#define pb push_back
#define inf 0x3f3f3f3f
using namespace std;
const int N=2e5+5;
const int M=5e5+5;
int L,R,m,U[M],V[M],du[N],st[N];
int fa[M],val[M],tag[M],res[M];
vector<int>v;
int find(int x){
    if(fa[x]==x)return x;
    int tmp=fa[x];
    fa[x]=find(fa[x]);
    val[x]^=val[tmp];
    return fa[x];
}
void upd(int x){
    find(x);
    if(!(val[x]^tag[fa[x]])){
        tag[fa[x]]^=1;
    }
}
void unionset(int x,int y){
    int u=find(x),v=find(y);
    if(u!=v){
        val[u]=val[x]^val[y]^1,fa[u]=v;
    }
}
void solve(vector<int>&v,int cur){
    if(v.size()==0)return;
    for(auto e:v){
        fa[e]=e,val[e]=0,tag[e]=0,st[U[e]]=st[V[e]]=0;
    }
    int f=0;
    for(auto e:v){
        int x=U[e],y=V[e];
        if(st[x]||st[y])f=1;
        if(!st[x]&&!st[y]){
            st[x]=st[y]=e;
        }
        else if(!st[x]){
            upd(st[y]);
            unionset(e,st[y]);
            st[x]=e,st[y]=0;
        }
        else if(!st[y]){
            upd(st[x]);
            unionset(e,st[x]);
            st[y]=e,st[x]=0;
        }
        else{
            upd(st[x]),upd(st[y]);
            unionset(e,st[x]),unionset(e,st[y]);
            st[x]=st[y]=0;
        }
    }
    if(f==0){
        for(auto e:v)res[e]=1;
    }
    else{
        vector<int>vl,vr;
        for(auto e:v){
            find(e);
            if(val[e]^tag[fa[e]]){
                vr.pb(e);
            }
            else{
                vl.pb(e);
            }
        }
        solve(vl,cur>>1),solve(vr,cur>>1);
        for(auto e:vr)res[e]+=cur>>1;
    }
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    cin>>L>>R>>m;
    for(int i=1;i<=m;i++){
        cin>>U[i]>>V[i],V[i]+=L;
        du[U[i]]++,du[V[i]]++,v.pb(i);
    }
    int dmax=0;
    for(int i=1;i<=L+R;i++)dmax=max(dmax,du[i]);
    int x=2;while(x<dmax)x<<=1;
    solve(v,x);
    cout<<x<<"\n";
    for(int i=1;i<=m;i++)cout<<res[i]<<"\n";
}
相关推荐
孞㐑¥12 分钟前
算法——BFS
开发语言·c++·经验分享·笔记·算法
Gary Studio1 小时前
rk芯片驱动编写
linux·学习
mango_mangojuice1 小时前
Linux学习笔记(make/Makefile)1.23
java·linux·前端·笔记·学习
工程师老罗2 小时前
YOLOv1 核心知识点笔记
笔记·yolo
lingggggaaaa2 小时前
安全工具篇&动态绕过&DumpLsass凭据&Certutil下载&变异替换&打乱源头特征
学习·安全·web安全·免杀对抗
PP东2 小时前
Flowable学习(二)——Flowable概念学习
java·后端·学习·flowable
学电子她就能回来吗2 小时前
深度学习速成:损失函数与反向传播
人工智能·深度学习·学习·计算机视觉·github
半壶清水3 小时前
[软考网规考点笔记]-软件开发、项目管理与知识产权核心知识与真题解析
网络·笔记·压力测试
tq10863 小时前
先探索,后设计
笔记
hnult4 小时前
2026 在线培训考试系统选型指南:核心功能拆解与选型逻辑
人工智能·笔记·课程设计