BZOJ 4205: 卡牌配对

16 篇文章 1 订阅
14 篇文章 0 订阅

难以置信地1A了 挺开心。。
这道题建图太神了 一开始就觉得要在素数搞一些很哲♂学的东西然而想不到QAQ
考虑到200以内只有46个素数 我们可以用一个二元组来表示其实两种属性值的质因子
比如(a,b)表示a是A的质因子 b是B的质因子

其次题目说至多一组属性互质 也就是至少两组有不互质(有公共质因子)
其实上面的二元组应该加多一元 表示是哪两个属性 也就是AB AC BC三种可能
所以就可以挺顺理成章地建图 源点 连到X类卡牌 连到中间至多3*46*46个表示质因子的点 再连到Y类卡牌 and 汇点

(偷黄学长的分析一波)因为200<2*3*5*7,所以两侧每个点至多连出3*3*3条边。于是我们构成了一个70000个点,2000000条边的网络流。
于是乎我们用dicnic愉快地跑完这个分层图 如下图
2333

我觉得我这个代码很优美啊~也没过100行

#include<bits/stdc++.h>
using namespace std;
const int N=70005,M=4000005,inf=1e9+7;
inline int read()
{
    char ch=getchar(); int x=0,f=1;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}
    return x*f;
}
struct edge{int y,c,next;}a[M]; int len,first[N];
void ins(int x,int y,int c){
    a[++len]=(edge){y,c,first[x]},first[x]=len;
    a[++len]=(edge){x,0,first[y]},first[y]=len;
}
queue<int>q; int h[N],t,p[50],pl,n,m; bool v[201];
bool bfs()
{
    for(int i=0;i<=t;i++)h[i]=0;
    q.push(0); h[0]=1;
    while(!q.empty()){
        int x=q.front(); q.pop();
        for(int k=first[x];k;k=a[k].next){
            int y=a[k].y;
            if(a[k].c && !h[y])
                h[y]=h[x]+1,q.push(y);
        }
    }
    return h[t]>0;
}
int dfs(int x,int f)
{
    if(x==t)return f;
    int de=0;
    for(int k=first[x];k;k=a[k].next){
        int y=a[k].y;
        if(h[y]==h[x]+1){
            int u=dfs(y,min(a[k].c,f-de));
            de+=u,a[k].c-=u,a[k^1].c+=u;
            if(f==de)return f;
        }
    }
    if(!de)h[x]=-1;
    return de;
}
vector<int>P[201];
void getp()
{
    for(int i=2;i<200;i++){
        if(!v[i])p[++pl]=i;
        for(int j=1;j<=pl && i*p[j]<200;j++)
        {
            v[i*p[j]]=1;
            if(i%p[j]==0)break;
        }
    }
    for(int i=2;i<=200;i++)for(int j=1;j<=pl && p[j]<=i;j++)
        if(i%p[j]==0)P[i].push_back(j);
}
inline int u(int x1,int x2,int k){return n+m+(x1-1)*pl+x2+k*pl*pl;}
int main()
{
    n=read(),m=read(); int i,j,k,x,y,z;
    getp();
    t=n+m+3*pl*pl+1; len=1;
    for(i=1;i<=n;i++)
    {
        ins(0,i,1);
        x=read(),y=read(),z=read();
        for(j=0;j<P[x].size();j++)for(k=0;k<P[y].size();k++)
            ins(i,u(P[x][j],P[y][k],0),1);
        for(j=0;j<P[x].size();j++)for(k=0;k<P[z].size();k++)
            ins(i,u(P[x][j],P[z][k],1),1);
        for(j=0;j<P[y].size();j++)for(k=0;k<P[z].size();k++)
            ins(i,u(P[y][j],P[z][k],2),1);
    }
    for(i=1;i<=m;i++)
    {
        ins(i+n,t,1);
        x=read(),y=read(),z=read();
        for(j=0;j<P[x].size();j++)for(k=0;k<P[y].size();k++)
            ins(u(P[x][j],P[y][k],0),i+n,1);
        for(j=0;j<P[x].size();j++)for(k=0;k<P[z].size();k++)
            ins(u(P[x][j],P[z][k],1),i+n,1);
        for(j=0;j<P[y].size();j++)for(k=0;k<P[z].size();k++)
            ins(u(P[y][j],P[z][k],2),i+n,1);
    }
    int ans=0;
    while(bfs())ans+=dfs(0,inf);
    printf("%d\n",ans);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值