BZOJ4784 [Zjoi2017]仙人掌
标签:树形DP,tarjan,仙人掌
题目
题目传送门
Description
如果一个无自环无重边无向连通图的任意一条边最多属于一个简单环,我们就称之为仙人掌。所谓简单环即不经过重复的结点的环。
现在九条可怜手上有一张无自环无重边的无向连通图,但是她觉得这张图中的边数太少了,所以她想要在图上连上一些新的边。同时为了方便的存储这张无向图,图中的边数又不能太多。经过权衡,她想要加边后得到的图为一棵仙人掌。不难发现合法的加边方案有很多,可怜想要知道总共有多少不同的加边方案。两个加边方案是不同的当且仅当一个方案中存在一条另一个方案中没有的边。
Input
多组数据,第一行输入一个整数T表示数据组数。
每组数据第一行输入两个整数n,m,表示图中的点数与边数。
接下来m行,每行两个整数u,v(1≤u,v≤n,u!=v)表示图中的一条边。保证输入的图联通且没有自环与重边
Sigma(n)<=5*10^5,m<=10^6,1<=m<=n * (n-1)/2
Output
对于每组数据,输出一个整数表示方案数,当然方案数可能很大,请对998244353取模后
输出。
Sample Input
2
3 2
1 2
1 3
5 4
1 2
2 3
2 4
1 5
Sample Output
2
8
对于第一组样例合法加边的方案有 {}, {(2,3)},共 2 种。
题意
给定一个无向连通图,添加一些边,使其变为仙人掌
问方案数
分析
我大概。。。只会写10分暴力+20分特判吧
正解如下
首先特判这张图是否是仙人掌,如果是,那么输出0
显然环是对答案不产生贡献的,那么破环为链
这样原图转化成一棵树,部分分设置的用意如此
用h[i]表示有i个点,它们之间匹配的方案数
h[i]=h[i-1]+(i-1)*h[i-2]
- 如果i节点不与其他节点相连,那么方案就是h[i−1]
- 如果与其他节点连接,那么有(i−1)中选择方式,而当选择一个点以后,有两个点不能再连接,那么方案就是(i−1)∗h[i−2]
用f[i]表示子树i的方案数,g[i]表示可向上拓展的子树的方案数,num表示子树的节点个数(即大小)
f[i]=h[num]*\prod g[son](son\in Child(i))
每一个子节点都可以向上扩展并相对独立,然后一共有h[num]种儿子的匹配方案
g[i]=f[i]+h[num-1]*num*\prod g[son](son\in Child(i))
该节点可以自己想上扩展为f[i],并且有num个子节点,每个子节点还可以选择一个儿子,并且一共有h[num-1]种匹配方案
code
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#define rep(i,a,b) for(register int i=a;i<=b;i++)
#define dep(i,a,b) for(register int i=a;i>=b;i--)
#define ll long long
#define mem(x,num) memset(x,num,sizeof x)
#define reg(x) for(int i=last[x];i;i=e[i].next)
using namespace std;
inline ll read()
{ll f=1,x=0;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;
}
const int mod=998244353,maxn=5e5+6,maxm=1e6+6;
struct edge{int to,next,from;}e[maxm*2];
bool isca,vis[maxn],flag[maxm*2];ll f[maxn],g[maxn],h[maxn],ans;
int que[maxn],dfn[maxn],low[maxn],col[maxn],tim,top,T,n,m,cnt,last[maxn];
void insert(int u,int v){e[++cnt]=(edge){v,last[u],u};last[u]=cnt;flag[cnt]=1;e[++cnt]=(edge){u,last[v],v};last[v]=cnt;flag[cnt]=1;
}
void tarjan(int x,int fa){bool flag=0;dfn[x]=low[x]=++tim;que[++top]=x;reg(x){if(e[i].to==fa)continue;if(!dfn[e[i].to]){tarjan(e[i].to,x);low[x]=min(low[x],low[e[i].to]);if(low[e[i].to]<dfn[x]){if(flag){isca=0;return;}flag=1;}}else{low[x]=min(low[x],dfn[e[i].to]);if(dfn[e[i].to]<dfn[x]){if(flag){isca=0;return;}flag=1;}}}if(low[x]==dfn[x])while(true){col[que[top--]]=x;if(que[top+1]==x)break;}
}
void dfs(int x,int fa){vis[x]=1;f[x]=1,g[x]=0;int num=0;reg(x){if(!flag[i]||e[i].to==fa)continue;dfs(e[i].to,x);f[x]=f[x]*g[e[i].to]%mod;num++;}g[x]=(f[x]*h[num]%mod+f[x]*h[num-1]%mod*num%mod)%mod;f[x]=f[x]*h[num]%mod;
}
int main()
{h[0]=h[1]=1;rep(i,2,maxn-5)h[i]=(h[i-1]+(h[i-2]*(i-1)%mod))%mod;T=read();while(T--){isca=1,ans=1;cnt=tim=top=0;n=read(),m=read();rep(i,1,n){vis[i]=0;last[i]=col[i]=dfn[i]=low[i]=0;}rep(i,1,m){int u=read(),v=read();insert(u,v);}tarjan(1,0);if(!isca){printf("0\n");continue;}rep(i,1,cnt)if(col[e[i].from]==col[e[i].to])flag[i]=0;rep(i,1,n){if(vis[i])continue;dfs(i,0);ans=ans*f[i]%mod;}printf("%lld\n",ans);}return 0;
}
BZOJ4784 [Zjoi2017]仙人掌相关推荐
- [zjoi2017]仙人掌
前言 谨以此题纪念我第一次参加省选时刚了5h这一题得到0分的经历 题目相关 链接 题目大意 给出仙人掌定义:如果一个无自环无重边无向连通图的任意一条边最多属于一个简单环,我们就称之为仙人掌 给出一个图 ...
- ●洛谷P3687 [ZJOI2017]仙人掌
题链: https://www.luogu.org/problemnew/show/P3687 题解: 计数DP,树形DP. (首先对于这个图来说,如果初始就不是仙人掌,那么就直接输出0) 然后由于本 ...
- 【题解】ZJOI2017仙人掌
感觉这题很厉害啊,虽然想了一天多但还是失败了--(:д:) 这题首先注意到给定图中如果存在环其实对于答案是没有影响的.然后关键之处就在于两个 \(dp\) 数组,其中 \(f[u]\) 表示以 \(u ...
- ZJOI2017 仙人掌
题目大概是给个无向图,求添加边使其变为仙人掌的方案数. 直接判断是否仙人掌,特判输出0即可. 否则的话,把环拆开成链变成一个树,考虑暴力计算儿子的配对方案数,打表规律:\(f[i] = f[i - 1 ...
- [BZOJ4784][UOJ290][ZJOI017]仙人掌-树形DP
仙人掌 Description 如果一个无自环无重边无向连通图的任意一条边最多属于一个简单环,我们就称之为仙人掌.所谓简单环即不经过重复的结点的环. 现在九条可怜手上有一张无自环无重边的无向连通图,但 ...
- 【ZJOI2017】仙人掌 题解
题目大意 ~~~~~~给出一个无重边无自环的无向连通图(n 个点 m 条边),问有多少种再往上加边的方案,使得新图是仙人掌. ~~~~~~多组数据, n<=5e5, ∑m\sum m ...
- 【ZJOI2017】仙人掌
题面 https://www.luogu.org/problem/P3687 题解 如果原图不是仙人掌(在这里,我们认为一棵树也是一个仙人掌),则无解,输出$0$. 如果原图是一个仙人掌,加的边是不能 ...
- Forest Program dfs+tanjar仙人掌
题目链接 CCPC2019 F题. 题意:给一颗仙人掌树,让你求每一个小环的边的个数,用快速幂即可求解. 思路:第一反应是tanjar乱搞,把每个环上的点取出来,类似于缩点的方法.但是忽然感觉dfs能 ...
- 基于卷积神经网络(CNN)的仙人掌图像分类
点击上方"小白学视觉",选择加"星标"或"置顶" 重磅干货,第一时间送达 今天我们的目标是建立一个分类器,将图像分类为"仙人掌&q ...
最新文章
- 在 Windows XP 中,无法使用 Windows 图片和传真查看器来查看图片
- javamail 收邮件并解析附件
- canvas系列教程02-直线和曲线
- JS回调函数(callback)
- fputc会覆盖吗_深覆盖 VS 深覆合如何区分?这一篇看完就都明白了!
- 【Vue.js 3.0】不使用根标签,Error:[vue/no-multiple-template-root] The template root requires exactly one ...
- PAT1070. 结绳
- 学一学番茄风控独创的体系内容
- java 设计模式学习笔记四 prototype模式
- 【转】otool 与dylib
- Android科大讯飞语音合成
- UVALive3713 Astronauts
- JS自动弹出广告窗口
- 【Android -- 开源库】BottomNavigation 的基本使用
- codeforces A. Parity
- Ubuntu系统 不使用snap安装firefox,使用firefox官网下载的安装包安装firefox
- java内网推送app,java – 使用来自app-engine的Parse发送推送通...
- 【计算机网络】因特网和互联网的区别
- 【用复制粘贴让脚本动起来(二)】根据画面点击特定坐标
- axure 图片切换图片的交互_Axure教程:在动态面板里面设置图片自动/手动切换
热门文章
- CXO是什么?医药板块龙头都有哪些?
- export,sorce
- 通信原理包络是什么意思_通信原理调制:频谱搬移 先冲模拟调制中调幅说起...
- 针对“使用U盘提示使用前需要将其格式化,无法访问U盘情况”的解决办法
- matlab上的神经网络模型导出成C文件
- 游戏运营数据分析指标
- 简化是解决复杂性问题的首要前提---真正的管理者是把复杂的问题变得简单化
- blender简单模型上色
- 12123查询违章显示服务器异常,我的交管12123以前可以处理违章 现在处理不了是什么情况 。提示服务异常...
- 解析三极管放大电路设计技巧