【学习笔记】[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";
}
相关推荐
dengqingrui12344 分钟前
【树形DP】AT_dp_p Independent Set 题解
c++·学习·算法·深度优先·图论·dp
我的心永远是冰冰哒1 小时前
ad.concat()学习
学习
ZZZ_O^O1 小时前
二分查找算法——寻找旋转排序数组中的最小值&点名
数据结构·c++·学习·算法·二叉树
slomay3 小时前
关于对比学习(简单整理
经验分享·深度学习·学习·机器学习
hengzhepa3 小时前
ElasticSearch备考 -- Async search
大数据·学习·elasticsearch·搜索引擎·es
小小洋洋5 小时前
BLE MESH学习1-基于沁恒CH582学习
学习
m0_689618286 小时前
水凝胶发生器,不对称设计妙,医电应用前景广
笔记
Ace'6 小时前
每日一题&&学习笔记
笔记·学习
IM_DALLA6 小时前
【Verilog学习日常】—牛客网刷题—Verilog进阶挑战—VL25
学习·fpga开发·verilog学习
挥剑决浮云 -6 小时前
Linux 之 安装软件、GCC编译器、Linux 操作系统基础
linux·服务器·c语言·c++·经验分享·笔记