Day19 二分图染色法，匈牙利算法（二分图匹配问题） 质数知识（试除法、分解质因数、埃及筛、线性筛）

二分图染色法
二分图就是说可以把图的点分成两部分，这两部分的点没有相连的边
就是说二分图就是不包含奇数条边的环的图
传送门
这个dfs也没啥好说的
要记得用链表要初始化，一定要写，写add前就要初始化h数组

#include<iostream>
#include<cstring>
using namespace std;
const int N=1e5+10,M=2e5+10;
int h[N],ne[M],e[M],idx,n,m,color[N];
void add(int a,int b){
    e[idx]=b;
    ne[idx]=h[a];
    h[a]=idx++;
}
bool dfs(int u,int c){
    color[u]=c;                             //染色
    for(int i=h[u];i!=-1;i=ne[i]){          //读取单链表
        int j=e[i];                         //读取下一个下标
        if(!color[j]){                      //如果未被染色
            if(!dfs(j,3-c))return 0;        //就染成3-c
        }
        else if(color[j]==c)return 0;       //如果下一个已经被染色且和自己相同颜色，就不是二分图
    }
    return 1;
}

int main(){
    cin>>n>>m;
    memset(h,-1,sizeof h);              //要初始化。。。
    while(m--){
        int a,b;
        cin>>a>>b;
        add(a,b),add(b,a);
    }
    bool flag=1;
    for(int i=1;i<=n;i++){              //逐个染色
        if(!color[i]){
            if(!dfs(i,1)){              //如果染色失败就不是二分图
                flag=0;
                break;
            }
        }
    }
    if(flag)puts("Yes");
    else puts("No");
    return 0;
}

---

二分图匹配问题
传送门

就是匈牙利算法
可以想象成配对问题
计算一个二分图能匹配的对的最大值

#include<iostream>
#include<cstring>
using namespace std;
const int N=510,M=1e5+10;

int n1,n2,m,h[N],e[M],ne[M],idx;
bool st[N];
int match[N];
void add(int a,int b){
    e[idx]=b;
    ne[idx]=h[a];
    h[a]=idx++;
}

bool find(int x){
    for(int i=h[x];i!=-1;i=ne[i]){
        int j=e[i];
        if(!st[j]){
            st[j]=1;
            if(match[j]==0||find(match[j])){            //右边的点没被匹配或者  匹配到的右边的点可以换
                match[j]=x;
                return 1;
            }
        }
    }
    return 0;
}

int main(){
    memset(h,-1,sizeof h);
    
    cin>>n1>>n2>>m;
    
    while(m--){
        int a,b;
        cin>>a>>b;
        add(a,b);
    }
    int res=0;
    for(int i=1;i<=n1;i++){
        memset(st,0,sizeof st);
        if(find(i))res++;
    }
    cout<<res;
    return 0;
}

---

试除法
要注意小于2的情况
传送门

#include<iostream>
using namespace std;

bool check(int x){
    if(x<=1)return 0;
    for(int i=2;i<=x/i;i++){
        if(x%i==0)return 0;
      //  cout<<i<<" "<<x<<endl;
    }
    return 1;
}

int main(){
    int t;
    cin>>t;
    
    while(t--){
        int n;
        cin>>n;
        if(check(n))puts("Yes");
        else puts("No");
    }
    return 0;
}

---

分解质因数

注意处理完后如果x仍然大于1那其就是一个质因数

#include<iostream>
using namespace std;

void fun(int x){
    for(int i=2;i<=x/i;i++){
        int res=0;
        while(x%i==0){
            res++;
            x/=i;
        }
        if(res)cout<<i<<" "<<res<<endl;
    }
    if(x>1){
        cout<<x<<" "<<1<<endl;
    }
    cout<<endl;
}

int main(){
    int t;
    cin>>t;
    while(t--){
        int n;
        cin>>n;
        fun(n);
    }
    return 0;
}

---

埃及筛
这个是只用质数筛选后方数

#include<iostream>
using namespace std;
const int N=1e6+10;
int prime[N];
bool st[N];
int ai(int n){
    int res=0;
    for(int i=2;i<=n;i++){
       if(!st[i]){
           prime[res++]=i;
           for(int j=2*i;j<=n;j+=i)st[j]=1;
       }
    } 
    return res;
}


int main(){
    int n;
    cin>>n;
    
    cout<<ai(n);
    return 0;
}

---

线性筛
线性筛就是每个数都最多只筛一次
所以是线性的。。。

#include<iostream>
using namespace std;
const int N=1e6+10;
int prime[N];
bool st[N];

int oula(int n){
    int res=0;
    for(int i=2;i<=n;i++){              //线性筛到n 
        if(!st[i])prime[res++]=i;       //如果没被筛就加入
        for(int j=0;prime[j]<=n/i;j++){ //记住i*prime[j]<=n
            st[i*prime[j]]=1;           //prime[j]是i*prime[j]的最小质因数
            if(i%prime[j]==0)break;     //如果i的最小质因数就是prime[j]了，那么此循环的最小质因数就永远不是prime[j]了
        }
    }
    return res;
}

int main(){
    int n;
    cin>>n;
    
    cout<<oula(n);
    return 0;
}