最近突然有很多人来问我这些题目怎么做OwO

然而并不是我出的,结论我也不是很懂

研究了一下觉得非常的一颗赛艇,于是就打算写这样一篇题解

DAG 1

我们考虑DAG的性质,枚举出度为0的点

设出度为0的点有k个,则一共有C(n,k)种方案

对于剩下的(n-k)个点和这k个点之间可以任意连边,方案为2^(k*(n-k))

去掉这k个点之后剩下(n-k)个点仍然是DAG,方案为f(n-k)

则方案数为C(n,k)*2^(k*(n-k))*f(n-k)

考虑到剩下的(n-k)个点中也可能有出度为0的点,即方案可能会有重复

考虑算重的部分,用容斥原理解决即可

最后得到f(n)=sigma((-1)^(k-1)*C(n,k)*2^(k*(n-k))*f(n-k))

这个递推式是O(n^2)的,可以通过DAG 1

DAG 2

不难想到要用FFT来优化上面的递推式

瓶颈在于2^(k*(n-k))的拆分

我们发现(n-k)^2=n^2+2nk+k^2

不难构造出k*(n-k)= n^2/2 - k^2/2 - (n-k)^2/2

代入之后得到卷积形式,直接做CDQ+FFT即可

定义g(n)=2^(n^2/2)*n!

注意到f(n)/g(n)=sigma( (-1)^(k-1)/g(k) * f(n-k)/g(n-k) )

定义多项式h,h(n)=(-1)^(n-1)/g(n)

定义多项式ans,ans(n)=f(n)/g(n)

变形之后得ans-ans*h=1

之后得到ans=1/(1-h)

多项式求逆即可

UPD:注意一下这里的实现,由于n^2/2有可能不是整数,而且(mod-1)和2不是互素的

但是因为在模意义下我们可以找到x^2=2(mod 998244353)

所以2^(n^2/2)=x^(n^2)

至于怎么求解x,因为可以预处理,所以直接枚举算出结果在代码里直接用就可以了

正常向的做法是用原根取对数之后BSGS搞一搞

DAG 3

设f(n)为n个点的DAG的个数(可以不连通)

设g(n)为n个点的连通DAG的个数

不难想到用f(n)减去不连通的就是连通的

计算不连通图常见技巧是枚举特定的某个点所在的联通块的大小

则得到g(n)=f(n)-sigma( C(n-1,k-1)*g(k)*f(n-k) )

这样我们就在O(n^2)的时间内求出来了

DAG 4

注意到上面的那个式子自然就是卷积形式

一发CDQ+FFT就可以了

但是我们可以做到更好

变形之后得到sigma( g(k)/(k-1)! *f(n-k)/(n-k)! ) = f(n)/(n-1)!

我们要求g,那么多项式求逆即可

我们注意到实际上不连通的DAG是由多个连通的DAG组成的

设可以不连通的DAG的多项式为f,连通DAG的多项式为g

不难得到f=e^g

那么g=ln(f)

多项式求ln即可

cojs DAG计数问题1-4 题解报告的更多相关文章

  1. cojs QAQ的图论题 题解报告

    话说这个题目应该叫做 斯特林数的逆袭 QAQ 先说一说部分分的算法 1.n<=5 直接暴力搜索就可以了 2.k=0的时候不难发现任意一张图的价值都是n,问题转化为计算有多少种图,显然是2^C(n ...

  2. cojs 简单的01串 题解报告

    题意显然是求n位二进制串中不大于其逆序串,取反串,逆序取反串的所有串按字典序排序后的第k个 由于n很小,k很大所以我们可以考虑逐位确定 问题转化为了求方案数,这显然是可以用数位DP做的 设f[len] ...

  3. cojs 简单的数位DP 题解报告

    首先这道题真的是个数位DP 我们考虑所有的限制: 首先第六个限制和第二个限制是重复的,保留第二个限制即可 第五个限制在转移中可以判断,不用放在状态里 对于第一个限制,我们可以增加一维表示余数即可 对于 ...

  4. cojs 二分图计数问题1-3 题解报告

    OwO 良心的FFT练手题,包含了所有的多项式基本运算呢 其中一部分解法参考了myy的uoj的blog 二分图计数 1: 实际是求所有图的二分图染色方案和 我们不妨枚举这个图中有多少个黑点 在n个点中 ...

  5. cojs 强连通图计数1-2 题解报告

    OwO 题目含义都是一样的,只是数据范围扩大了 对于n<=7的问题,我们直接暴力搜索就可以了 对于n<=1000的问题,我们不难联想到<主旋律>这一道题 没错,只需要把方程改一 ...

  6. 2015浙江财经大学ACM有奖周赛(一) 题解报告

    2015浙江财经大学ACM有奖周赛(一) 题解报告 命题:丽丽&&黑鸡 这是命题者原话. 题目涉及的知识面比较广泛,有深度优先搜索.广度优先搜索.数学题.几何题.贪心算法.枚举.二进制 ...

  7. 题解报告:hdu 1398 Square Coins(母函数或dp)

    Problem Description People in Silverland use square coins. Not only they have square shapes but also ...

  8. 题解报告:hdu 2069 Coin Change(暴力orDP)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2069 Problem Description Suppose there are 5 types of ...

  9. 题解报告:hdu 1028 Ignatius and the Princess III(母函数or计数DP)

    Problem Description "Well, it seems the first problem is too easy. I will let you know how fool ...

随机推荐

  1. pandas生成时间列表(某段连续时间或者固定间隔时间段)

    python生成一个日期列表 首先导入pandas import pandas as pd def get_date_list(begin_date,end_date): date_list = [x ...

  2. Docker添加镜像加速器

    Docker默认pull连接镜像为国外镜像,速度较慢,注册阿里云可以生成一个镜像加速器 登录阿里云 https://cr.console.aliyun.com获取私有加速地址 修改配置文件/etc/d ...

  3. Scikit-learn技巧(拓展)总结

    Scikit-learn技巧(拓展)总结 本文转载自:http://www.jianshu.com/p/516f009c0875 最近看了<Python数据挖掘入门与实战>,网上有说翻译地 ...

  4. iOS判断UIView是否显示在屏幕上

    @interface - (BOOL)isDisplayedInScreen; @end @implementation UIView(UIScreenDisplaying) //判断View是否显示 ...

  5. qt——QT中QWidget、QDialog及QMainWindow的区别

    QWidget类是所有用户界面对象的基类. 窗口部件是用户界面的一个基本单元:它从窗口系统接收鼠标.键盘和其它事件,并且在屏幕上绘制自己.每一个窗口部件都是矩形的,并且它们按Z轴顺序排列.一个窗口部件 ...

  6. vscode主题配色

    https://www.cnblogs.com/garvenc/p/vscode_customize_color_theme.html

  7. android 打开系统相机,

    1.第一步在androidmanifest.xml中注册 <uses-permission android:name="android.permission.WRITE_EXTERNA ...

  8. 3.cassandra遇到内存占用过高的问题

    目前cssandra的内存分配如下: https://docs.datastax.com/en/cassandra/2.1/cassandra/operations/ops_tune_jvm_c.ht ...

  9. golang 的精髓--pipeline流水线,对现实世界的完美模拟

    https://blog.golang.org/pipelines https://www.cnblogs.com/junneyang/p/6215785.html 简介 Go语言的并发原语允许开发者 ...

  10. Scala集合类详解

    对scala中的集合类虽然有使用,但是一直处于一知半解的状态.尤其是与java中各种集合类的混合使用,虽然用过很多次,但是一直也没有做比较深入的了解与分析.正好趁着最近项目的需要,加上稍微有点时间,特 ...