HDU 5486(Difference of Clustering-聚类)

Difference of Clustering

给定10^6个关系
现在要求从左到右1对多,多对1,1:1的情况的数量
这里写图片描述
但是n:n的情况要排除
这里写图片描述

直接离散,求连通块

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=next[p])  
#define Lson (x<<1)
#define Rson ((x<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define MAXN (2000000+10)
#define MAXM (2000000+10) 
#define mp make_pair
#define pb push_back
#define fi first
#define se second
typedef long long ll;

int u[MAXN],v[MAXN],n,a[MAXN];
map<int,int> h;

int sz;
bool b[MAXN];
int cal1,cal2;

class link_table
{
public: 
    void mem()
    {
        MEM(pre) MEM(edge) MEM(pre) size=1;
    }
    int edge[MAXM],next[MAXM],pre[MAXN],size;

    void addedge(int u,int v)
    {
        edge[++size]=v;
        next[size]=pre[u];  
        pre[u]=size;
    }
    void addedge2(int u,int v){addedge(u,v);addedge(v,u);}


    void dfs(int x)
    {
        b[x]=1;
        if (x<=n) ++cal1;else ++cal2;
        Forp(x) {
            int v=edge[p];
            if (!b[v]) dfs(v);
        }

    } 
}St;
int main()
{
    freopen("C.in","r",stdin);
//  freopen(".out","w",stdout);

    int T; cin>>T;
    For(kcase,T) {
        int A=0,B=0,C=0;
        cin>>n;
        For(i,n) {
            int x,y;
            scanf("%d%d",&x,&y);    
            u[i]=x,v[i]=y;
        }
        sz=0;
        For(i,n) a[++sz]=u[i],a[++sz]=v[i];
        sort(a+1,a+1+sz); 
        sz=unique(a+1,a+1+sz)-(a+1);
        h.clear();
        For(i,sz) h[a[i]]=i;
        For(i,n) u[i]=h[u[i]],v[i]=h[v[i]];

        St.mem();
        For(i,n)
        {
            St.addedge2(u[i],v[i]+n);
        }
        For(i,n*2) b[i]=0;
        For(i,n*2) {
            if (!b[i]) 
            {
                cal1=cal2=0;
                St.dfs(i);
                if (cal1==1&&cal2==1) C++;
                else if (cal1==1&&cal2>1) A++;
                else if (cal2==1&&cal1>1) B++;

            }
        }       
        printf("Case #%d: %d %d %d\n",kcase,A,B,C);
    }


    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值