热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

【UOJ】#37.【清华集训2014】主旋律

题解一道,神奇的题我们考虑正难则反,我们求去掉这些边后有多少图不是强连通的怎么求呢,不是强连通的图缩点后一定是一个DAG,并

题解

一道,神奇的题= =

我们考虑正难则反,我们求去掉这些边后有多少图不是强连通的

怎么求呢,不是强连通的图缩点后一定是一个DAG,并且这个DAG里面有两个点

1005017-20180528161738173-663740640.png

我们想一下,如果我们把1当成入度为0的点,随便造出个图,可以是这个图吧
如果把2当成入度为0的点,随便造出个图,也可以是这个图吧
把1和2当成入度为0的点,随便造出个图,还可以是这个图吧……

那么这像什么,容斥啊

以下的点说的都是缩点后的点

奇数个入度为0的点就是+,偶数个入度为0的点就是-

那么我们就有了一个精妙的容斥!
\(f[S]\)为点集S是强连通分量的方案数
\(g[S]\) ……是……
定义可能很奇怪,是容斥出来的,点集S的系数

嗯,这么考虑,我们知道了S这个点全是入度为0的点,那么剩下的造图的过程,就是S这个点集向外延伸的边,和除S中的点之外的点之间的边,有或没有,乘上一个2的指数幂就好

那么,怎么知道S这个系数呢???

我们考虑S中包含一个编号最小的点u的强连通分量,且入度为0

假如这个点的集合是T,T包含u

那么\(g[S] = g[S] - g[S ^ T] * f[T]\),咦,为什么是减法
因为\(g[S ^ T]\)容斥的系数,在多了一个联通块后,所有的奇偶性都改变了,那么+-号也取反了,所以是减法

最后,我们用每个S的子集(包括S)算一遍不合法的图,用总共的图减掉就行,就是\(f[S]\)的值
还有\(g[S] += f[S]\)这样的话,代表整个S是一个强连通分量,S缩点后构成的入度为0的点是1个

\(h[S]\)表示S这个点集中的边
\(w[T]\)表示在全集是S的情况下,选择T这个点集作为缩点后入度为0的点,能向外延伸的边集

代码

#include
#include
#include
#include
#include
#include
#include
//#define ivorysi
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define mo 974711
#define MAXN 100005
#define RG register
using namespace std;
typedef long long int64;
typedef double db;
template
void read(T &res) {res &#61; 0;char c &#61; getchar();T f &#61; 1;while(c <&#39;0&#39; || c > &#39;9&#39;) {if(c &#61;&#61; &#39;-&#39;) f &#61; -1;c &#61; getchar();}while(c >&#61; &#39;0&#39; && c <&#61; &#39;9&#39;) {res &#61; res * 10 &#43; c - &#39;0&#39;;c &#61; getchar();}res *&#61; f;
}
template
void out(T x) {if(x <0) {putchar(&#39;-&#39;);x &#61; -x;}if(x >&#61; 10) {out(x / 10);}putchar(&#39;0&#39; &#43; x % 10);
}
const int MAXS &#61; 1 <<15;
const int MOD &#61; 1000000007;
int N,M;
int cnt[MAXS &#43; 5],In[MAXS &#43; 5],Out[MAXS &#43; 5],h[MAXS &#43; 5],w[MAXS &#43; 5],g[MAXS &#43; 5],f[MAXS &#43; 5],pow2[405];
inline int inc(int a,int b) {a &#61; a &#43; b;if(a >&#61; MOD) a -&#61; MOD;return a;}
inline int mul(int a,int b) {return 1LL * a * b % MOD;}
void Init() {read(N);read(M);pow2[0] &#61; 1;for(int i &#61; 1 ; i <&#61; M ; &#43;&#43;i) {pow2[i] &#61; pow2[i - 1] <<1;if(pow2[i] >&#61; MOD) pow2[i] -&#61; MOD;}for(int i &#61; 1 ; i }
void Solve() {for(int S &#61; 1 ; S <(1 <}
int main() {
#ifdef ivorysifreopen("f1.in","r",stdin);
#endifInit();Solve();return 0;
}

转载于:https://www.cnblogs.com/ivorysi/p/9100856.html


推荐阅读
  • 本题涉及一棵由N个节点组成的树(共有N-1条边),初始时所有节点均为白色。题目要求处理两种操作:一是改变某个节点的颜色(从白变黑或从黑变白);二是查询从根节点到指定节点路径上的第一个黑色节点,若无则输出-1。 ... [详细]
  • 扫描线三巨头 hdu1928hdu 1255  hdu 1542 [POJ 1151]
    学习链接:http:blog.csdn.netlwt36articledetails48908031学习扫描线主要学习的是一种扫描的思想,后期可以求解很 ... [详细]
  • Splay Tree 区间操作优化
    本文详细介绍了使用Splay Tree进行区间操作的实现方法,包括插入、删除、修改、翻转和求和等操作。通过这些操作,可以高效地处理动态序列问题,并且代码实现具有一定的挑战性,有助于编程能力的提升。 ... [详细]
  • 题目Link题目学习link1题目学习link2题目学习link3%%%受益匪浅!-----&# ... [详细]
  • 本文详细探讨了KMP算法中next数组的构建及其应用,重点分析了未改良和改良后的next数组在字符串匹配中的作用。通过具体实例和代码实现,帮助读者更好地理解KMP算法的核心原理。 ... [详细]
  • 技术分享:从动态网站提取站点密钥的解决方案
    本文探讨了如何从动态网站中提取站点密钥,特别是针对验证码(reCAPTCHA)的处理方法。通过结合Selenium和requests库,提供了详细的代码示例和优化建议。 ... [详细]
  • Java 中的 BigDecimal pow()方法,示例 ... [详细]
  • 1:有如下一段程序:packagea.b.c;publicclassTest{privatestaticinti0;publicintgetNext(){return ... [详细]
  • 本文介绍了Java并发库中的阻塞队列(BlockingQueue)及其典型应用场景。通过具体实例,展示了如何利用LinkedBlockingQueue实现线程间高效、安全的数据传递,并结合线程池和原子类优化性能。 ... [详细]
  • 深入理解Cookie与Session会话管理
    本文详细介绍了如何通过HTTP响应和请求处理浏览器的Cookie信息,以及如何创建、设置和管理Cookie。同时探讨了会话跟踪技术中的Session机制,解释其原理及应用场景。 ... [详细]
  • 使用 Azure Service Principal 和 Microsoft Graph API 获取 AAD 用户列表
    本文介绍了一段通用代码示例,该代码不仅能够操作 Azure Active Directory (AAD),还可以通过 Azure Service Principal 的授权访问和管理 Azure 订阅资源。Azure 的架构可以分为两个层级:AAD 和 Subscription。 ... [详细]
  • XNA 3.0 游戏编程:从 XML 文件加载数据
    本文介绍如何在 XNA 3.0 游戏项目中从 XML 文件加载数据。我们将探讨如何将 XML 数据序列化为二进制文件,并通过内容管道加载到游戏中。此外,还会涉及自定义类型读取器和写入器的实现。 ... [详细]
  • 本文详细介绍了如何构建一个高效的UI管理系统,集中处理UI页面的打开、关闭、层级管理和页面跳转等问题。通过UIManager统一管理外部切换逻辑,实现功能逻辑分散化和代码复用,支持多人协作开发。 ... [详细]
  • 本文探讨了如何在模运算下高效计算组合数C(n, m),并详细介绍了乘法逆元的应用。通过扩展欧几里得算法求解乘法逆元,从而实现除法取余的计算。 ... [详细]
  • 本文介绍如何使用 NSTimer 实现倒计时功能,详细讲解了初始化方法、参数配置以及具体实现步骤。通过示例代码展示如何创建和管理定时器,确保在指定时间间隔内执行特定任务。 ... [详细]
author-avatar
燕门雪_346
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有