[BZOJ1051] [HAOI2006]受欢迎的牛(强连通)

传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1051


一开始想网络流去了(跪。。。
后来发现肥肠简单,只用强连通缩点之后看一下出度为0的分量个数,如果只有一个那就合法,这个强连通分量里面的所有点都可以被其他点到达,所以输出这个分量里面的点数,否则(出度为0的连通分量大于1)就输出0


code:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=10010;
const int maxm=50010;
struct node
{
    int x,y,next;
}a[maxm*4],b[maxm*4]; 
int len,last[maxn];
int elen,elast[maxn];
int st,ed;
void ins(int x,int y)
{
    len++; 
    a[len].x=x;a[len].y=y;
    a[len].next=last[x];last[x]=len;
}/*
void eins(int x,int y)
{
    elen++;
    b[elen].x=x;b[elen].y=y;
    b[elen].next=elast[x]; elast[x]=elen;
}*/
int low[maxn],dfn[maxn],id;
int sta[maxn],top; bool insta[maxn];
int bl[maxn],scc; 
int list[maxm*4],head,tail;
int sum[maxn];
void dfs(int x)
{
    dfn[x]=low[x]=++id;
    sta[++top]=x; insta[x]=true;
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(dfn[y]==-1)
        {
            dfs(y);
            low[x]=min(low[x],low[y]);
        }
        else if(insta[y]) low[x]=min(low[x],dfn[y]);
    }
    if(low[x]==dfn[x])
    {
        scc++; int i;
        do{
            i=sta[top--];
            bl[i]=scc;
            insta[i]=false;
            sum[scc]++;
        }while(i!=x);
    }
}
int chu[maxn];
void rebuild()
{
    for(int i=1;i<=len;i++)
    {
        if(bl[a[i].x]!=bl[a[i].y])
        {
            //eins(bl[a[i].x],bl[a[i].y]);
            chu[bl[a[i].x]]++;
        }
    }
}
int main()
{
    int n,m;scanf("%d%d",&n,&m);
    memset(last,0,sizeof(last)); len=0;
    for(int i=1;i<=m;i++)
    {
        int x,y;scanf("%d%d",&x,&y);
        ins(x,y);
    }
    memset(dfn,-1,sizeof(dfn));
    id=top=scc=0;
    for(int i=1;i<=n;i++)
        if(dfn[i]==-1)
            dfs(i);
    //memset(elast,0,sizeof(elast)); elen=0;
    rebuild();
    int cnt=0,tmp=0;
    for(int i=1;i<=scc;i++)
    {
        if(chu[i]==0)
        {
            cnt++; if(cnt>1){printf("0\n"); return 0;}
            tmp=sum[i];
        }
    }
    printf("%d\n",tmp);
    return 0;
}
### 关于 BZOJ1728 Two-Headed Cows (双头) 的算法解析 此问题的核心在于如何通过有效的图论方法解决给定约束下的最大独立集问题。以下是详细的分析和解答。 #### 问题描述 题目要求在一个无向图中找到最大的一组节点集合,使得这些节点之间满足特定的颜色匹配条件。具体来说,每条边连接两个节点,并附带一种颜色标记(A 或 B)。对于任意一条边 \(u-v\) 和其对应的颜色 \(c\),如果这条边属于最终选取的子集中,则必须有至少一个端点未被选入该子集或者两端点均符合指定颜色关系。 #### 解决方案概述 本题可以通过 **二分枚举 + 图染色验证** 来实现高效求解。核心思想如下: 1. 假设当前最优解大小为 \(k\),即尝试寻找是否存在一个大小为 \(k\) 的合法子集。 2. 枚举每一个可能作为起点的节点并将其加入候选子集。 3. 对剩余部分执行基于 BFS/DFS 的图遍历操作,在过程中动态调整其他节点的状态以确保整体合法性。 4. 如果某次试探能够成功构建符合条件的大规模子集,则更新答案;反之则降低目标值重新测试直至收敛至最佳结果。 这种方法利用了贪心策略配合回溯机制来逐步逼近全局最优点[^1]。 #### 实现细节说明 ##### 数据结构设计 定义三个主要数组用于记录状态信息: - `color[]` : 存储每个顶点所分配到的具体色彩编号; - `used[]`: 表示某个定点是否已经被处理过; - `adjList[][]`: 记录邻接表形式表示的原始输入数据结构便于后续访问关联元素。 ##### 主要逻辑流程 ```python from collections import deque def check(k, n): def bfs(start_node): queue = deque([start_node]) used[start_node] = True while queue: u = queue.popleft() for v, c in adjList[u]: if not used[v]: # Assign opposite color based on edge constraint 'c' target_color = ('B' if c == 'A' else 'A') if color[u]==c else c if color[v]!=target_color and color[v]!='?': return False elif color[v]=='?': color[v]=target_color queue.append(v) used[v] =True elif ((color[u]==c)==(color[v]==('B'if c=='A'else'A'))): continue return True count=0 success=True for i in range(n): if not used[i]: temp_count=count+int(color[i]=='?' or color[i]=='A') if k<=temp_count: color_copy=color[:] if bfs(i): count=temp_count break else : success=False return success n,m=list(map(int,input().split())) colors=[['?']*m]*n for _ in range(m): a,b,c=input().strip().split() colors[int(a)-1].append((int(b),c)) low ,high,res=0,n,-1 while low<=high: mid=(low+high)//2 color=['?']*n used=[False]*n if check(mid,n): res=mid low=mid+1 else : high=mid-1 print(res) ``` 上述代码片段展示了完整的程序框架以及关键函数 `check()` 的内部运作方式。它接受参数 \(k\) 并返回布尔值指示是否有可行配置支持如此规模的选择[^2]。 #### 复杂度分析 由于采用了二分查找技术缩小搜索空间范围再加上单轮 DFS/BFS 时间复杂度 O(V+E),总体性能表现良好适合大规模实例运行需求。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值