博客
关于我
BZOJ 3812 主旋律
阅读量:283 次
发布时间:2019-03-01

本文共 1822 字,大约阅读时间需要 6 分钟。

题目链接

题解

考虑求非强连通子图的数量,假设为 g g g,那么答案就是 2 m − g 2^m-g 2mg。现在考虑求 g g g

假设 f s f_s fs表示用 s s s这些点能构成的强连通图的个数, g s g_s gs表示用 s s s这些点能构成的非强连通图的方案数,其中构成 i i i个强连通分量则对 g s g_s gs的贡献为 ( − 1 ) i (-1)^i (1)i。容易发现

g s = f s − ∑ t ⊂ s , u ∈ t g t f s − t g_s=f_s-\sum_{t\subset s,u\in t} g_tf_{s-t} gs=fsts,utgtfst
那么
f s = 2 e s − ∑ t ⊆ s , t ̸ = ∅ 2 e s − t + e s − t , t g s f_s=2^{e_s}-\sum_{t\subseteq s,t\not= \varnothing}2^{e_{s-t}+e_{s-t,t}}g_s fs=2ests,t̸=2est+est,tgs
容易发现, f s f_s fs此时需要的是不包含 f s f_s fs g s g_s gs,因此 g s g_s gs在求出 f s f_s fs之前是不能 + f s +f_s +fs的。

代码

#include 
int read(){ int x=0,f=1; char ch=getchar(); while((ch<'0')||(ch>'9')) { if(ch=='-') { f=-f; } ch=getchar(); } while((ch>='0')&&(ch<='9')) { x=x*10+ch-'0'; ch=getchar(); } return x*f;} const int maxn=15;const int maxm=1<
=mod) { pow[i]-=mod; } } f[0]=g[0]=1; for(int s=1; s<=full; ++s) { int sk=s^lowbit(s); for(int t=sk; t; t=sk&(t-1)) { g[s]-=1ll*f[s^t]*g[t]%mod; if(g[s]<0) { g[s]+=mod; } } for(int i=1; i<=n; ++i) { if((1<<(i-1))&s) { in[s]+=ecnt[s][i]; } } f[s]+=pow[in[s]]; for(int t=s; t; t=s&(t-1)) { int e=0; for(int i=1; i<=n; ++i) { if((1<<(i-1))&t) { e+=ecnt[s^t][i]; } } f[s]-=1ll*pow[e+in[s^t]]*g[t]%mod; if(f[s]<0) { f[s]+=mod; } } g[s]+=f[s]; if(g[s]>=mod) { g[s]-=mod; } } printf("%d\n",f[full]); return 0;}

转载地址:http://scwo.baihongyu.com/

你可能感兴趣的文章
Nacos服务注册与发现demo
查看>>
Nacos服务注册与发现的2种实现方法!
查看>>
nacos服务注册和发现原理简单实现案例
查看>>
Nacos服务注册总流程(源码分析)
查看>>
nacos服务注册流程
查看>>
Nacos服务部署安装
查看>>
nacos本地可以,上服务器报错
查看>>
Nacos注册Dubbo(2.7.x)以及namespace配置
查看>>
Nacos注册中心有几种调用方式?
查看>>
nacos注册失败,Feign调用失败,feign无法注入成我们的bean对象
查看>>
nacos源码 nacos注册中心1.4.x 源码 nacos源码如何下载 nacos 客户端源码下载地址 nacos discovery下载地址(一)
查看>>
nacos源码 nacos注册中心1.4.x 源码 spring cloud alibaba 的discovery做了什么 nacos客户端是如何启动的(二)
查看>>
nacos源码 nacos注册中心1.4.x 源码 如何注册服务 发送请求,nacos clinet客户端心跳 nacos 注册中心客户端如何发送的心跳 (三)
查看>>
Nacos源码分析:心跳机制、健康检查、服务发现、AP集群
查看>>
nacos看这一篇文章就够了
查看>>
Nacos简介、下载与配置持久化到Mysql
查看>>
Nacos简介和控制台服务安装
查看>>
Nacos管理界面详细介绍
查看>>
Nacos编译报错NacosException: endpoint is blank
查看>>
nacos自动刷新配置
查看>>