$bzoj2560$ 串珠子 容斥+$dp$
正解:容斥+$dp$
解题报告:
$umm$虽然题目蛮简练的了但还是有点难理解,,,我再抽象一点儿,就说有$n$个点,点$i$和点$j$之间有$a_{i,j}$条无向边可以连,问有多少种方案可以连成一张联通图
显然考虑容斥呗?设$f_i$表示状态为$i$的点连成联通图的合法方案,$g_i$表示状态为$i$的点随便连边的所有方案
显然$g_i$可以先预处理出来?就等于$\prod_{u,v\in i}a_{u,v}$.然后$f_i$就等于$g_i$减去不合法的数量.不合法数量显然就考虑枚举子集${i}'$,就等于$\sum f_{{i}'}\cdot g_{i-{i}'}$.
但是这样显然依然会有锅,即一个不合法方案会被枚举其包含的联通块次.为了保证不重不漏,就只用枚指定点的联通块大小,比较通常的做法就枚举最大/最小点的联通块大小,也就钦定${i}'$中包含了最大/最小的点
然后就做完了$QwQ$
$over$
因为一些不知名原因我本机$AC$,$BZOJ$上$WA$了(事实上是,$emacs\ AC$,$lemon\ WA$,$darkbzoj\ WA$,$QAQ$
但是我暂时懒得搞了先把代码放上来趴$kk$
#include<bits/stdc++.h>
using namespace std;
#define il inline
#define lf double
#define int long long
#define ll long long
#define gc getchar()
#define ri register int
#define rc register char
#define rb register bool
#define lowbit(x) (x&(-x))
#define rp(i,x,y) for(ri i=x;i<=y;++i)
#define my(i,x,y) for(ri i=x;i>=y;--i)
#define gdgs(i,x) for(ri i=x-lowbit(x);i;i-=lowbit(i)) const int N=,mod=;
int n,a[N][N],lg[<<N],tot,d[N],cnt;
ll g[<<N],f[<<N],re[N]; il int read()
{
rc ch=gc;ri x=;rb y=;
while(ch!='-' && (ch>'' || ch<''))ch=gc;
if(ch=='-')ch=gc,y=;
while(ch>='' && ch<='')x=(x<<)+(x<<)+(ch^''),ch=gc;
return y?x:-x;
} signed main()
{
freopen("2560.in","r",stdin);freopen("2560.out","w",stdout);
n=read();rp(i,,n-)rp(j,,n-)a[i][j]=read();tot=(<<n)-;rp(i,,n-)lg[<<i]=i;g[]=;
rp(i,,tot)
{ll tmp=;gdgs(j,i)tmp=1ll*tmp*(a[lg[lowbit(i)]][lg[lowbit(j)]]+)%mod;g[i]=g[i-lowbit(i)]*tmp%mod;}
rp(i,,tot)
{
cnt=;gdgs(j,i)d[cnt++]=lowbit(j);
rp(j,,(<<cnt)-)re[j]=re[j-lowbit(j)]|d[lg[lowbit(j)]],f[i]=(f[i]+f[i^re[j]]*g[re[j]])%mod;
f[i]=(g[i]-f[i]+mod)%mod;
}
printf("%lld\n",f[tot]);
return ;
}
随机推荐
- Celery后台任务
Celery 在程序运行过程中,经常会遇到一些耗时耗资源的任务,为了避免这些任务阻塞主进程的运行,我们会采用多线程或异步任务去处理.比如在Web中需要对新注册的用户发一封激活邮件来验证账户,而发邮件本 ...
- 如何安装 btsync
本文告诉大家如何在 windows 和 Linux 安装使用 Btsync 而且分享一些小东西给大家 btsync 是分布式网盘 在这高速运作的信息化时代,使用云端来衔接工作和生活的点滴已是寻常事.可 ...
- JQuery完整验证&密码的显示与隐藏&验证码
HTML <link href="bootstrap.css" rel="stylesheet"> <link href="gloa ...
- oracle WHERE子句中的连接顺序
ORACLE采用自下而上的顺序解析WHERE子句,根据这个原理,表之间的连接必须写在其他WHERE条件之前, 那些可以过滤掉最大数量记录的条件必须写在WHERE子句的末尾. 例如: (低效,执行时间1 ...
- @省选模拟赛03/16 - T3@ 超级树
目录 @description@ @solution@ @accepted code@ @details@ @description@ 一棵 k-超级树(k-SuperTree) 可按如下方法得到:取 ...
- ∆ (triangle)
2.1 题目描述 给定一个无自环重边的无向图,求这个图的三元环1的个数以及补图2的三元环个数. 2.2 输入格式 第一行 2 个数 n, m ,分别表示图的点数.边数. 接下来 m 行,每行两个数 u ...
- hdu 1556 Color the ball(区间更新,单点求值)
Color the ball Time Limit: 9000/3000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)To ...
- seo优化:302跳转变为301跳转
<?php header("Location: http://www.XXX.com",TRUE,301); exit; ?> 要加exit;不加exit;还是会是30 ...
- Spring Boot JPA 懒加载
最近在使用spring jpa 的过程中经常遇到懒加载的错误:"` org.hibernate.LazyInitializationException: could not initiali ...
- pycharm下的多个python版本共存(一)
经历过IDLE,anaconda,和pycharn的编程环境,并进行了一段时间的项目编程后,决定使用pycharm作为以后的工作环境. 一方面因为项目组其他人推荐,另一方面在使用过程中比较顺手.当然很 ...