首页 > 其他分享> 文章详细

2022江苏 B. Prime Ring Plus

2022-07-28 12:03:35 阅读：303 来源： 互联网

标签：Prime 匹配 int 奇点建图 Plus 偶点 Ring

分析

网络流，直接说建图了。

简单分析后可以发现，长度为奇数的环一定不合法，因为相邻的奇偶性一定要不同
可以发现是经典模型，一个奇点可以与两个偶点匹配，求最大匹配，最大匹配必须将所有的点都包括
则我们可以分奇偶建图，从源点向所有的奇点连一条容量为2的边，从所有偶点向汇点连一条容量为2的边，奇点与偶点相加为素数的连接一条容量为1的边

确定方案

这个也比较简单了，我们可以从1开始，找到一条匹配的边，则代表两者相邻，然后转移到匹配的边的终点，然后不断递归下去，就能找到一整个环

一点碎碎念，其实这题最难的是想到，只要有合法方案，那我们一定可以将其合成一个环

AC_code

#include<bits/stdc++.h>

using namespace std;

const int N = 1e4 + 10,M = 2e7,INF = 0x3f3f3f3f;

int h[N],ne[M],e[M],f[M],idx;
int cur[N],d[N],q[N];
bool st[N<<1];
int primes[N<<1],cnt;
int n,S,T;

void init(int x)
{
    st[x] = 1;
    for(int i=2;i<x;i++)
    {
        if(!st[i]) primes[cnt++] = i;
        for(int j=0;primes[j]*i<x;j++)
        {
            st[primes[j]*i] = 1;
            if(i%primes[j]==0) break;
        }
    }
}

void add(int a,int b,int c)
{
    e[idx] = b,ne[idx] = h[a],f[idx] = c,h[a] = idx++;
    e[idx] = a,ne[idx] = h[b],f[idx] = 0,h[b] = idx++;
} 

bool bfs()
{
    int hh = 0,tt = -1;
    memset(d,0x3f,sizeof d);
    q[++tt] = S,d[S] = 1,cur[S] = h[S];
    while(hh<=tt)
    {
        int t = q[hh++];
        for(int i=h[t];~i;i=ne[i])
        {
            int j = e[i];
            if(d[j]==INF&&f[i])
            {
                d[j] = d[t] + 1;
                cur[j] = h[j];
                if(j==T) return 1;
                q[++tt] = j;
            }
        }
    }
    return 0;
}

int find(int u,int limit)
{
    if(u==T) return limit;
    int flow = 0;
    for(int i=cur[u];~i&&flow<limit;i=ne[i])
    {
        int j = e[i];
        cur[u] = i;
        if(d[j]==d[u]+1&&f[i])
        {
            int t = find(j,min(limit-flow,f[i]));
            if(!t) d[j] = -1;
            f[i] -= t,f[i^1] += t,flow += t;
        }
    }
    return flow;
}

int dinic()
{
    int r = 0,flow;
    while(bfs()) if(flow = find(S,INF)) r += flow;
    return r;
}

int main()
{
    init(N<<1);
    scanf("%d",&n);
    memset(h,-1,sizeof h);
    S = 0,T = n + 1;
    for(int i=1;i<=n;i+=2) add(S,i,2);
    for(int i=1;i<=n;i+=2)
        for(int j=2;j<=n;j+=2)
            if(!st[i+j])
                add(i,j,1);
    for(int i=2;i<=n;i+=2) add(i,T,2);
    if(dinic()!=n) cout<<"-1\n";
    else 
    {
        cout<<"1\n"<<n<<" ";
        memset(st,0,sizeof st);
        int t = 1;
        while(1)
        {
            cout<<t;
            st[t] = 1;
            bool F = 0;
            for(int i=h[t];~i;i=ne[i])
            {
                int j = e[i];
                if(j==T||j==S) continue;
                if(st[j]) continue;
                F = 1;
                t = j;
                break;
            }
            if(!F) break;
            else cout<<" ";
        }
    }
    return 0;
}

标签：Prime,匹配,int,奇点,建图,Plus,偶点,Ring
来源： https://www.cnblogs.com/aitejiu/p/16528129.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

2022江苏 B. Prime Ring Plus

分析

AC_code