[我可能是个SB啊] UOJ #290 【ZJOI2017】仙人掌

最新推荐文章于 2020-03-29 11:58:00 发布
最新推荐文章于 2020-03-29 11:58:00 发布
阅读量823 收藏
点赞数
分类专栏： DP
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u014609452/article/details/65644103
版权

这好像不难啊 考场上智商掉线了吧 可能是脑袋被吃掉了吧

先说正经的
这个题 要是不是棵仙人掌 那么就无解
要是是棵树 直接树形DP
那仙人掌呢 把所有环上的边拆掉 对剩下的森林中的每棵树DP 然后乘法原理就好了
然后 做完了…

还是不知道考场上在干什么

我可能忍不住分享下自己的悲惨经历啊
考场上想出了树形DP 然后就没有然后了？ 尼玛仙人掌这么好想 就被我弃了？！
行吧 考挂不如狗 下午切这道题 结果先是常量 N 和 M 不分 数组大小乱开？
然后是数据范围看错？
UOJ上目前全部70分都是我交的？
我真是不知道自己在干什么

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<set>
using namespace std;
typedef long long ll;
typedef pair<int,int> abcd;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc());
}

const int N=500005;
const int M=1000005;
const int P=998244353;
const int INV2=(P+1)>>1;

ll fac[N],inv[N],pw[N];

inline void Pre(int n){
  fac[0]=1; for (int i=1;i<=n;i++) fac[i]=fac[i-1]*i%P;
  inv[1]=1; for (int i=2;i<=n;i++) inv[i]=(ll)(P-P/i)*inv[P%i]%P;
  inv[0]=1; for (int i=1;i<=n;i++) inv[i]=inv[i-1]*inv[i]%P;
  pw[0]=1; for (int i=1;i<=n;i++) pw[i]=pw[i-1]*INV2%P;
}

struct edge{
  int u,v,next;
}G[M<<1];
int head[N],inum;
inline void add(int u,int v,int p){
  G[p].u=u; G[p].v=v; G[p].next=head[u]; head[u]=p;
}

int tag[N];
int fat[N],depth[N];

#define V G[p].v

inline void dfs_j(int u,int fa){
  fat[u]=fa; depth[u]=depth[fa]+1;
  for (int p=head[u];p;p=G[p].next)
    if (V!=fa)
      if (!depth[V])
    dfs_j(V,u),tag[u]+=tag[V];
      else if (depth[V]<depth[u])
    tag[u]++,tag[V]--;
}

int fap[N]; bool ed[M];

inline void find(int u,int fa){
  fat[u]=fa; depth[u]=depth[fa]+1;
  for (int p=head[u];p;p=G[p].next)
    if (V!=fa)
      if (!depth[V])
    fap[V]=p>>1,find(V,u);
      else if (depth[V]<depth[u]){
    int t=u; ed[p>>1]=1;
    while (t!=V)
      ed[fap[t]]=1,t=fat[t];
      } 
}

bool vst[N]; ll f[N];

inline ll C(int n,int a){
  return fac[n]*inv[n-2*a]%P*inv[a]%P*pw[a]%P;
}

inline void dfs(int u,int fa){
  vst[u]=1; int cnt=0;
  f[u]=1;
  for (int p=head[u];p;p=G[p].next)
    if (V!=fa)
      dfs(V,u),f[u]=f[u]*f[V]%P,cnt++;
  ll tmp=0;
  for (int i=0;i<=cnt/2;i++)
    tmp+=C(cnt,i)*(!fa?1:(cnt-2*i+1))%P;
  tmp%=P;
  f[u]=f[u]*tmp%P;
}

int n,m;
int x[M],y[M];

int main(){
  int T;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  read(T); Pre(500000);
  while (T--){
    read(n); read(m);
    for (int i=1;i<=n;i++) head[i]=0; inum=1;
    for (int i=1;i<=m;i++){
      read(x[i]),read(y[i]); if (x[i]>y[i]) swap(x[i],y[i]);
      add(x[i],y[i],++inum),add(y[i],x[i],++inum);
    }

    for (int i=1;i<=n;i++) depth[i]=fat[i]=tag[i]=0;
    dfs_j(1,0);
    int flag=0; for (int i=1;i<=n;i++) flag|=tag[i]>1;
    if (flag){
      printf("0\n"); continue;
    }

    for (int i=1;i<=n;i++) depth[i]=fat[i]=0;
    for (int i=1;i<=m;i++) ed[i]=0;
    find(1,0);

    for (int i=1;i<=n;i++) head[i]=0; inum=1;
    for (int i=1;i<=m;i++)
      if (!ed[i])
    add(x[i],y[i],++inum),add(y[i],x[i],++inum);

    for (int i=1;i<=n;i++) f[i]=vst[i]=0;
    ll ans=1;
    for (int i=1;i<=n;i++)
      if (!vst[i])
    dfs(i,0),ans=ans*f[i]%P;

    printf("%lld\n",ans);
  }
  return 0;
}

UPD 再次被自己傻哭 根本不用树形DP啊
直接 ans=Π Fdegi 就好了 F <script type="math/tex" id="MathJax-Element-51">F</script>表示配对的方案数 因为每个点的度怎么配是不相互影响的 反正就是怎么连都只可能是链

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<set>
using namespace std;
typedef long long ll;
typedef pair<int,int> abcd;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc());
}

const int N=500005;
const int M=1000005;
const int P=998244353;
const int INV2=(P+1)>>1;

ll fac[N],inv[N],pw[N];
ll F[N];

inline void Pre(int n){
  fac[0]=1; for (int i=1;i<=n;i++) fac[i]=fac[i-1]*i%P;
  inv[1]=1; for (int i=2;i<=n;i++) inv[i]=(ll)(P-P/i)*inv[P%i]%P;
  inv[0]=1; for (int i=1;i<=n;i++) inv[i]=inv[i-1]*inv[i]%P;
  pw[0]=1; for (int i=1;i<=n;i++) pw[i]=pw[i-1]*INV2%P;
  F[0]=F[1]=1; for (int i=2;i<=n;i++) F[i]=(F[i-1]+F[i-2]*(i-1)%P)%P;
}

struct edge{
  int u,v,next;
}G[M<<1];
int head[N],inum;
inline void add(int u,int v,int p){
  G[p].u=u; G[p].v=v; G[p].next=head[u]; head[u]=p;
}

int tag[N];
int fat[N],depth[N];

#define V G[p].v

inline void dfs_j(int u,int fa){
  fat[u]=fa; depth[u]=depth[fa]+1;
  for (int p=head[u];p;p=G[p].next)
    if (V!=fa)
      if (!depth[V])
    dfs_j(V,u),tag[u]+=tag[V];
      else if (depth[V]<depth[u])
    tag[u]++,tag[V]--;
}

int fap[N]; bool ed[M];

inline void find(int u,int fa){
  fat[u]=fa; depth[u]=depth[fa]+1;
  for (int p=head[u];p;p=G[p].next)
    if (V!=fa)
      if (!depth[V])
    fap[V]=p>>1,find(V,u);
      else if (depth[V]<depth[u]){
    int t=u; ed[p>>1]=1;
    while (t!=V)
      ed[fap[t]]=1,t=fat[t];
      } 
}

int deg[N];

int n,m;
int x[M],y[M];

int main(){
  int T;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  read(T); Pre(500000);
  while (T--){
    read(n); read(m);
    for (int i=1;i<=n;i++) head[i]=0; inum=1;
    for (int i=1;i<=m;i++){
      read(x[i]),read(y[i]); if (x[i]>y[i]) swap(x[i],y[i]);
      add(x[i],y[i],++inum),add(y[i],x[i],++inum);
    }

    for (int i=1;i<=n;i++) depth[i]=fat[i]=tag[i]=0;
    dfs_j(1,0);
    int flag=0; for (int i=1;i<=n;i++) flag|=tag[i]>1;
    if (flag){
      printf("0\n"); continue;
    }

    for (int i=1;i<=n;i++) depth[i]=fat[i]=0;
    for (int i=1;i<=m;i++) ed[i]=0;
    find(1,0);

    for (int i=1;i<=n;i++) deg[i]=0;
    for (int i=1;i<=m;i++)
      if (!ed[i])
    deg[x[i]]++,deg[y[i]]++;

    ll ans=1;
    for (int i=1;i<=n;i++)
      ans=ans*F[deg[i]]%P;

    printf("%lld\n",ans);
  }
  return 0;
}
里阿奴摩西
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
LOJ #2252. 「ZJOI2017」多项式（倍增）
Freopen的博客
10-26 408
题目 思路简单代码毒瘤的经典简单OI题，真的只需要会倍增 题解主要就是一句话：一个(mod2)\pmod 2(mod2)意义下的多项式的平方就是每一项的指数×=2\times =2×=2，那么可以考虑使用递归形式的快速幂来利用可快速平方和乘一个原多项式的性质，又因为需求串长度<=18<=18<=18，所以我们可以考虑存长度为需求串长的所有串的个数，并在倍增时简单转移所有串的个数，...
【LOJ2250】「ZJOI2017」仙人掌
cz_xuyixuan的博客
06-23 300
【题目链接】 点击打开链接 【思路要点】 考虑给定的图为大小为 N+1N+1N+1 的点的菊花图的情况，记方案数为 dpNdp_NdpN​ ，考虑最后一个儿子是否连边，则有转移 dpi=dpi−1+(i−1)dpi−2dp_i=dp_{i-1}+(i-1)dp_{i-2}dpi​=dpi−1​+(i−1)dpi−2​ 断开环边，剩余问题即为如何解决图是一棵树的情况。 注意到由于不允许存...
【BZOJ4784】[ZJOI2017]仙人掌（Tarjan，动态规划）
weixin_30414245的博客
02-17 110
【BZOJ4784】[ZJOI2017]仙人掌（Tarjan，动态规划） 题面 BZOJ 洛谷 题解 显然如果原图不是仙人掌就无解。 如果原图是仙人掌，显然就是把环上的边给去掉，变成若干森林连边成为仙人掌的方案数。 那么对于一棵树而言，考虑其变成仙人掌的方案数。 设\(a_i\)表示匹配\(i\)个儿子的方案数，显然转移时\(a_i=a_{i-1}+(i-1)*a_{i-2}\)，即考虑新加入的儿...
【ZJOI2017】仙人掌 题解
Four's
03-24 2369
题目大意      ~~~~~~给出一个无重边无自环的无向连通图（n 个点 m 条边），问有多少种再往上加边的方案，使得新图是仙人掌。       ~~~~~~多组数据， n<=5e5， ∑m\sum m<=1e6
ZJOI 2017 仙人掌 题解
a_forever_dream的博客
03-04 234
题目传送门 题目大意： 给出一张无向无自环无重边图，你可以加任意条边，问有多少种加边方式使得最后的图是一张仙人掌。 题解 首先判断这张图是不是仙人掌，如果不是那么可以直接输出 000 了。判断方法：先跑一次 tarjantarjantarjan，然后利用 dfsdfsdfs 树可以统计出图中每条边在多少个环里面，每个环会被正反两次统计，也就是说，假如一条边被统计了超过两次，就一定在一个以上的环内了...
BZOJ 2125: 最短路 （仙人掌，树链剖分）
EM LGH
03-29 170
第一道仙人掌题. 由于仙人掌中每条边最多只属于一个环，所以两个在环中的点的最短距离是好算的. code: #include <bits/stdc++.h> #define N 200006 #define setIO(s) freopen(s".in","r",stdin) using namespace std; int edges,n,m,Q,...
【ZJOI2017】仙人掌
psc233的博客
03-21 151
题目大意：给你一张联通图，问你加边使得他是个仙人掌的方案数 解析： 首先如果一开始就不是个仙人掌，那方案数显然是0 否则，我们发现由于每条边都只能被至多覆盖一次（仙人掌的性质），所以多出来的边必定可以把这棵树分割成一些独立的树，那么我们把问题转化为对于一颗树来说，求出他加边使得他是个仙人掌的方案数 那么接下来我们思考一下dp方程 设dp[i][0]表示当前这个点i的子树内部不存在向上连的...
uoj-data-converter.zip
09-21
UOJ Data Converter
南京大学服务器上搭建UOJ系统的保姆级步骤
最新发布
09-30
从申请南京大学服务器开始，然后搭建docker系统，之后创建UOJ系统用于发布课堂作业，用于帮助助教进行UOJ系统搭建，帮助助教更好的完成老师布置的任务
UOJ-System:通用在线法官（社区版）
03-11
支持添加多个评估机。 过渡搜索，局部放置，任意页面至少快速到达。 所有译文从编译，运行到评分，都可以由出题人自定义。 ♡额外测试和Hack机制，更严谨，更多乐趣。 支持OI / IOI / ACM等比赛模式；比赛内部设有...
UOJ228基础数据结构练习题.cpp
01-11
UOJ228基础数据结构练习题.cpp
unh6uoj.github.io
03-07
欢迎来到GitHub Pages 您可以使用的来维护和预览Markdown文件中网站的内容。 每当您提交到该存储库时，GitHub Pages都将运行从Markdown文件中的内容重建站点中的页面。 降价促销 Markdown是一种轻巧且易于使用的...
省选知识点总结（咕咕列表）
romiqi
02-21 388
一.数据结构 学习按重要程度降序 加强按熟练程度升序 待学习： 可并堆 树上莫队 可持久化trie 树套树 KD-tree 块状链表 舞蹈链 待加强： 虚树 线段树合并 单调队列单调栈 线段树高级操作 可持久化一系列 树链剖分（重，长） 二.图论 待学习： 上下界费用流 kruscal重构树 仙人掌 生成树拓展（度限制，最小树形图） 支配树 斯坦纳树...
【LOJ2252】「ZJOI2017」多项式
cz_xuyixuan的博客
06-26 939
【题目链接】 点击打开链接 【思路要点】 模 222 意义下，多项式平方相等于将每一项的指数翻倍。 我们需要维护一些关于字符串 SSS 的信息，使得其能够进行 “指数翻倍” 和 “指数翻倍并乘以原多项式” 的操作。 首先考虑测试点 6,7,86,7,86,7,8 。 考虑维护字符串 SSS 中 ∣S∣−17|S|-17∣S∣−17 个长度为 181818 的字符串中，每一种字符串的个数...
「ZJOI2017」字符串 （Hash）（分块）（线段树）
FSYo 的博客
02-15 375
LOJLOJLOJ 传送门 题解： 我们考虑用线段树维护最小后缀的出现位置，那么需要考虑如何合并两个区间 发现最小后缀可能是前面区间的某一个后缀加上后面的一整个字符串，注意这个某一个后缀并不一定是前面区间的最小后缀，于是我们需要维护可能后缀的集合。 定义：一个 iii 的 “kkk-后缀” 指的是字符串 S[i...k]S[i...k]S[i...k]（以下用 SiS_iSi​ 表示），一个 i...
UOJ#290. 【ZJOI2017】仙人掌
wxh010910
03-24 890
链接:http://uoj.ac/problem/290 树形DP (披着仙人掌的外衣) 去掉环就变成了若干棵树,分别DP即可 没有重边的性质非常好,这样我们就可以对于没有被覆盖的边,变成重边,每条边都要被覆盖,就可以DP了 复杂度是线性的 #include #define xx first #define yy second #defin
ZJOI2017 仙人掌 转化模型后的简单树形dp
YxuanwKeith的博客
03-26 1588
题目大意给定一个nn个点，mm条变的无向无自环的连通图，问都多少种加边方案使得加完边的图是一幅没有重边仙人掌。(即满足任意一条边只属于一个简单环中的无向无自环图的连通图) 多组数据。∑n≤5∗105\sum_n \leq 5*10^5 ∑m≤106\sum_m \leq 10^6解题思路首先讨论给定的图是树的情况。（转化一） 要求不能有边存在与两个简单环内相当于我们要加入新的边去覆盖这个树，是
ZJOI2017 仙人掌
Akakii的博客
03-25 1453
题解如果一开始的图就不是仙人掌，答案显然为0，可以Tarjan判断。 环显然不能产生贡献，所以可以把环边都断开。 现在模型转化为，给定一棵树，用路径去覆盖树上的每一条边，且路径不能相交，求方案数。 设fif_i表示做完了ii的子树，且没有路径可以向上扩展。 设gig_i表示做完了ii的子树，且有路径可以向上扩展。 设hih_i表示有ii个点，它们之间匹配的方案数。 记numnum为点xx
[BZOJ4784][Zjoi2017]仙人掌（树形DP）
女装的你，如此好看！
10-14 238
Address 洛谷P3687 BZOJ4784 UOJ#290 LOJ#2250 Solution 首先，如果原图不是仙人掌，就直接输出 000 。 否则对图进行一遍 DFS ，找出所有的环并去掉，原图变成森林。答案显然是每棵树的答案的乘积。 考虑到在如果在树上的点对 (u,v)(u,v)(u,v) 之间连边，就相当于把 uuu 到 vvv 的路径上的所有边标记为环边。 但存在两个问题： （1）...
UOJ#435题目“Simple Tree”
05-30
UOJ#435题目“Simple Tree”要求实现一种数据结构，支持以下两种操作： 1. 修改某个节点的权值； 2. 查询某个节点到根节点路径上权值的最大值。 因为需要查询某个节点到根节点的路径，所以可以考虑使用树链剖分...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值