【学习笔记】[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";
}
相关推荐
晴空对晚照27 分钟前
[动手学习深度学习]12.权重衰退
人工智能·深度学习·学习
宇寒风暖34 分钟前
HTML嵌入CSS样式超详解(尊享)
前端·css·笔记·学习·html
Ting丶丶1 小时前
Electron入门笔记
javascript·笔记·electron
虾球xz2 小时前
游戏引擎学习第149天
人工智能·学习·游戏引擎
奕天者2 小时前
C++学习笔记(十七)——类之封装
c++·笔记·学习
亭墨2 小时前
linux0.11内核源码修仙传第六章——中断初始化
linux·c语言·驱动开发·学习·面试
好多大米2 小时前
[资源分享]-web3/区块链/学习路线/资料/找工作方式/水龙头
学习·web3·区块链
李长渊哦3 小时前
学习文章:Spring Boot 中如何使用 `@Async` 实现异步处理
java·spring boot·学习
李长渊哦4 小时前
学习文章:Spring Boot 中 Redis 配置与序列化管理
spring boot·redis·学习
你的微笑像拥抱4 小时前
HTTP 强 Etag 和 弱 Etag
笔记