【HNOI/AHOI2018】毒瘤【容斥】【虚树/动态dp】

题意:nnn 个点 mmm 条边的连通无向图的独立集个数模 998244353998244353998244353

n≤105,m≤n+10n\leq 10^5,m\leq n+10n105,mn+10

为什么标题要把两个算法写一起?因为这两个东西在这类问题上是本质相同的,这也是写这篇博客的原因。

显然是随便找棵树出来 dp,然后枚举非树边容斥。

dp(u,0/1)dp(u,0/1)dp(u,0/1) 表示考虑 uuu 的子树,uuu 不选/选的方案

dp(u,0)=∏v∈son(u)(dp(v,0)+dp(v,1))dp(u,0)=\prod_{v\in son(u)}(dp(v,0)+dp(v,1))dp(u,0)=vson(u)(dp(v,0)+dp(v,1))

dp(u,1)=∏v∈son(u)dp(v,0)dp(u,1)=\prod_{v\in son(u)}dp(v,0)dp(u,1)=vson(u)dp(v,0)

用动态 dp 维护就可以了。但修改操作是乘 000 和撤销乘 000,你需要记录乘的 000 的个数,不是很好搞。这里介绍等价的虚树的做法。

考虑加入一个儿子 vvv,写成矩阵的形式

[dp′(u,0)dp′(u,1)]=[dp(v,0)+dp(v,1)00dp(v,0)][dp(u,0)dp(u,1)]\left[ \begin{matrix} dp'(u,0) \\ dp'(u,1) \end{matrix}\right]= \left[ \begin{matrix} dp(v,0)+dp(v,1) & 0\\ 0&dp(v,0) \end{matrix} \right] \left[ \begin{matrix} dp(u,0)\\ dp(u,1) \end{matrix} \right] [dp(u,0)dp(u,1)]=[dp(v,0)+dp(v,1)00dp(v,0)][dp(u,0)dp(u,1)]

这其实可以不用矩阵直接算出来。

因为有关于关键点的修改,我们希望维护这些信息来快速得到需要的 dp 值。所以我们建出关键点的虚树,维护虚树上每条边的转移系数。因为儿子的信息是不确定的(或者说是动态的),并且中途有加的操作,所以这里必须用矩阵。

对于这条边对应的原树路径,每向上走一步 [dp(u,0)dp(u,1)]T\left[\begin{matrix}dp(u,0)&dp(u,1)\end{matrix}\right]^T[dp(u,0)dp(u,1)]T 产生的贡献就会变成 [dp(u,0)+dp(u,1)dp(u,0)]T\left[\begin{matrix}dp(u,0)+dp(u,1)&dp(u,0)\end{matrix}\right]^T[dp(u,0)+dp(u,1)dp(u,0)]T ,对应转移矩阵为

S=[1110]S=\left[\begin{matrix}1&1\\1&0\end{matrix}\right]S=[1110]

然后每走一步再加上不在虚树上的儿子的贡献。

设得到的转移矩阵为 trans(u)trans(u)trans(u),每次更新的答案为 ans(u)ans(u)ans(u) 那么

ans(u)=f(u)∏(u,v)∈VirtualTreetrans(v)ans(v)ans(u)=f(u)\prod_{(u,v)\in VirtualTree}trans(v)ans(v)ans(u)=f(u)(u,v)VirtualTreetrans(v)ans(v)

其中 f(u)f(u)f(u)uuu 的不在虚树上的儿子的贡献。

开始听信题解区 O(n)\Omicron(n)O(n) 建虚树挂了,怎么都调不出来,然后极度愤怒的情况下换成 O(klog⁡k)\Omicron(k\log k)O(klogk) 就过了……这个故事告诉我们考场上不要自己 yy 奇怪的算法。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <utility>
#include <algorithm>
#include <vector>
#define MAXN 100005
#define MAXM 200005
using namespace std;
inline int read()
{int ans=0;char c=getchar();while (!isdigit(c)) c=getchar();while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();return ans;
}
const int MOD=998244353;
inline int add(const int& x,const int& y){return x+y>=MOD? x+y-MOD:x+y;}
inline int dec(const int& x,const int& y){return x<y? x-y+MOD:x-y;}
typedef long long ll;
inline int qpow(int a,int p)
{int ans=1;while (p){if (p&1) ans=(ll)ans*a%MOD;a=(ll)a*a%MOD,p>>=1;}return ans;
}
#define inv(x) qpow(x,MOD-2)
struct mat
{int e[2][2];inline mat(const int& v=0){e[0][0]=e[1][1]=v,e[0][1]=e[1][0]=0;}inline int* operator [](const int& i){return e[i];}inline const int* operator [](const int& i)const{return e[i];}	
};
mat S;
inline mat operator *(const mat& a,const mat& b)
{mat c;for (int i=0;i<2;i++)for (int k=0;k<2;k++)for (int j=0;j<2;j++)c[i][j]=(c[i][j]+(ll)a[i][k]*b[k][j])%MOD;return c;
}
typedef pair<int,int> pi;
vector<int> e[MAXN];
vector<pi> ex;
int dep[MAXN],fa[MAXN][20],siz[MAXN],key[MAXN],dfn[MAXN],tim;
int F[MAXN][2];
void dfs(int u,int f)
{dfn[u]=++tim,dep[u]=dep[fa[u][0]=f]+1;for (int i=1;i<20;i++) fa[u][i]=fa[fa[u][i-1]][i-1];F[u][0]=F[u][1]=1;for (int i=0;i<(int)e[u].size();i++)if (e[u][i]!=f){if (dfn[e[u][i]]) {if (dfn[e[u][i]]<dfn[u]) continue;ex.push_back(make_pair(u,e[u][i]));e[u][i]=f;continue;}dfs(e[u][i],u);F[u][0]=(ll)add(F[e[u][i]][0],F[e[u][i]][1])*F[u][0]%MOD;F[u][1]=(ll)F[e[u][i]][0]*F[u][1]%MOD;}
}
inline int lca(int x,int y)
{if (dep[x]<dep[y]) swap(x,y);int t=dep[x]-dep[y];for (int i=0;(1<<i)<=t;i++) if (t&(1<<i)) x=fa[x][i];if (x==y) return x;for (int i=19;i>=0;i--)if (fa[x][i]!=fa[y][i])x=fa[x][i],y=fa[y][i];return fa[x][0];
}
mat dp[MAXN];
vector<int> lis,son[MAXN];
int col[MAXN],tmp[MAXN][2];
inline bool cmp(const int& x,const int& y){return dfn[x]<dfn[y];}
const int Rt=1;
int main()
{S[0][0]=S[0][1]=S[1][0]=1;int n,m;n=read(),m=read();for (int i=1;i<=m;i++){int u,v;u=read(),v=read();e[u].push_back(v),e[v].push_back(u);}dfs(Rt,0);for (int i=0;i<(int)ex.size();i++)	lis.push_back(ex[i].first),lis.push_back(ex[i].second);lis.push_back(Rt);sort(lis.begin(),lis.end(),cmp);int s=(int)lis.size(); for (int i=0;i<s-1;i++) lis.push_back(lca(lis[i],lis[i+1]));sort(lis.begin(),lis.end(),cmp);lis.erase(unique(lis.begin(),lis.end()),lis.end());	for (int i=0;i<(int)lis.size();i++) key[lis[i]]=1;for (int i=0;i<(int)lis.size();i++)	for (int u=lis[i];u;++siz[u],u=fa[u][0]);for (int T=0;T<(int)lis.size();T++){mat cur=mat(1);for (int v=lis[T],u=fa[v][0];u;v=u,u=fa[u][0]){if (key[u]){son[u].push_back(lis[T]);dp[lis[T]]=cur;break;}cur=S*cur;for (int i=0;i<(int)e[u].size();i++)if (!siz[e[u][i]]&&e[u][i]!=fa[u][0]){int a=add(F[e[u][i]][0],F[e[u][i]][1]),b=F[e[u][i]][0];cur[0][0]=(ll)cur[0][0]*a%MOD,cur[0][1]=(ll)cur[0][1]*a%MOD;cur[1][0]=(ll)cur[1][0]*b%MOD,cur[1][1]=(ll)cur[1][1]*b%MOD;}}}for (int t=0;t<(int)lis.size();t++){int u=lis[t];F[u][0]=F[u][1]=1;for (int i=0;i<(int)e[u].size();i++)if (!siz[e[u][i]]&&e[u][i]!=fa[u][0]){F[u][0]=(ll)add(F[e[u][i]][0],F[e[u][i]][1])*F[u][0]%MOD;F[u][1]=(ll)F[e[u][i]][0]*F[u][1]%MOD;}}int k=(int)ex.size();int ans=0;
//	for (int i=0;i<k;i++) printf("%d %d\n",ex[i].first,ex[i].second);
//	for (int i=1;i<=n;i++) printf("%d %d\n",F[i][0],F[i][1]);for (int S=0;S<(1<<k);S++){int w=0;for (int i=0;i<k;i++) col[ex[i].first]=col[ex[i].second]=0;for (int i=0;i<k;i++) if (S&(1<<i)) col[ex[i].first]=col[ex[i].second]=1,w^=1;for (int T=(int)lis.size()-1;T>=0;T--){int u=lis[T];tmp[u][0]=F[u][0],tmp[u][1]=F[u][1];for (int i=0;i<(int)son[u].size();i++){tmp[u][0]=(ll)add(tmp[son[u][i]][0],tmp[son[u][i]][1])*tmp[u][0]%MOD;tmp[u][1]=(ll)tmp[son[u][i]][0]*tmp[u][1]%MOD;}if (col[u]) tmp[u][0]=0;if (u!=Rt){mat t;t[0][0]=tmp[u][0],t[1][0]=tmp[u][1];t=dp[u]*t;tmp[u][0]=t[0][0],tmp[u][1]=t[1][0];}}int t=add(tmp[Rt][0],tmp[Rt][1]);if (w) ans=dec(ans,t);else ans=add(ans,t);
//		printf("%d\n",t);}cout<<ans;return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314802.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

MediatR-进程内的消息通信框架

MediatR是一款进程内的消息订阅、发布框架&#xff0c;提供了Send方法用于发布到单个处理程序、Publish方法发布到多个处理程序&#xff0c;使用起来非常方便。目前支持 .NET Framework4.5、.NET Stardand1.3、.NET Stardand2.0等版本&#xff0c;可跨平台使用。要在项目中使用…

Codeforces Round #586 (Div. 1 + Div. 2) D. Alex and Julian 数学 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个无限个点的坐标轴&#xff0c;一个集合BBB&#xff0c;如果存在∣i−j∣bk|i-j|b_k∣i−j∣bk​的话&#xff0c;那么i,ji,ji,j之间就连边。现在问你至少要从集合BBB中去掉多少个数才能使得连完边之…

【十二省联考2019】字符串问题【后缀自动机】【拓扑排序】

题意&#xff1a;给一个字符串 SSS&#xff0c;以子串的形式给出一些 A 类串和 B 类串以及 mmm 对 A 类串支配 B 类串的关系。求一个总长度最长的 A 类串序列&#xff0c;使得每个串都存在一个 B 类串前缀被后一个串支配。无穷输出 −1-1−1。 ∣S∣,m≤2105|S|,m\leq 2\times …

Codeforces Round #586 (Div. 1 + Div. 2) B. Multiplication Table 思维 + 公式

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个n∗nn*nn∗n的矩阵&#xff0c;每个位置由ai∗aja_i*a_jai​∗aj​得来&#xff0c;主对角线为000&#xff0c;让你求出来aia_iai​。 n≤1e3n\le1e3n≤1e3 思路&#xff1a; 由公式ai,j∗ai,kaj,…

不好意思,这么久没有更新《从零开始掌握ASP.NET Core 》

点击上方蓝字&#xff0c;关注「我们」等了快个月了&#xff0c;终于开始更新了。因为感冒&#xff0c;弄的嗓子有点沙哑。所以停了半个月才是更新&#xff0c;目前一口气更新了12个章节&#xff0c;大家可以耐心观看内容了。《从零开始学ASP.NET Core 》-- 更新通知视频课程更…

【NOI2018】你的名字【后缀自动机】【可持久化线段树合并】【乱搞】

题意&#xff1a;给一个串 SSS&#xff0c;qqq 次询问&#xff0c;每次给定串 TTT 和 l,rl,rl,r &#xff0c;求有多少个本质不同的串是 TTT 的子串而不是 Sl…rS_{l\dots r}Sl…r​ 的子串。 ∣S∣≤5105,q≤105,∑∣T∣≤106|S|\leq 5\times 10^5,q\leq 10^5,\sum|T|\leq 10^…

2021牛客暑期多校训练营1 A.Alice and Bob 博弈 SG函数

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 有两堆石子&#xff0c;两个人每次可以进行如下操作&#xff1a;从某一堆狮子中拿出x(x>0)x(x>0)x(x>0)个&#xff0c;从另一堆石子中拿出s∗x(s>0)s*x(s>0)s∗x(s>0)个。谁不能操作谁输&…

【CTSC2010】珠宝商【后缀自动机】【点分治】【根号分治】

题意&#xff1a;给一棵 nnn 个点的树&#xff0c;每个点有个字符&#xff0c;另给一个长度为 mmm 的特征串&#xff0c;求树上 n2n^2n2 条有向路径在特征串中出现的次数之和。 n,m≤5104n,m\leq 5\times 10^4n,m≤5104 看到母串先建 SAM &#xff08;bushi 树上路径统计问题…

使用Azure云原生构建博客是怎样一种体验?(下篇)

点击上方蓝字关注“汪宇杰博客”接上篇《使用Azure云原生构建博客是怎样一种体验&#xff1f;&#xff08;上篇&#xff09;》DNSAzure DNS 是一套分布全球的域名解析服务。具有超高可用性和接近实时的记录更新及生效速度。我的博客也使用了这项服务。Azure 现在可以提供域名注…

2021牛客暑期多校训练营1 G Game of Swapping Numbers 思维 + 巧妙的转换

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个数组A,BA,BA,B&#xff0c;你可以选择AAA的两个位置i,j,i<ji,j,i<ji,j,i<j交换Ai,AjA_i,A_jAi​,Aj​&#xff0c;需要交换正好kkk次&#xff0c;问你最大的∑i1n∣Ai−Bi∣\sum_{i1}^n|A_…

.NET Core 3.0之深入源码理解HttpClientFactory(一)

写在前面创建HttpClient实例的时候&#xff0c;在内部会创建HttpMessageHandler链&#xff0c;我们知道HttpMessageHandler是负责建立连接的抽象处理程序&#xff0c;所以HttpClient的维护实际上就是维护HttpMessageHandler的使用&#xff0c;释放HttpClient并不会及时释放连接…

WTM 构建DotNetCore开源生态,坐而论道不如起而行之

作为一个8岁开始学习编程&#xff0c;至今40岁的老程序员&#xff0c;这辈子使用过无数种语言&#xff0c;从basic开始&#xff0c;到pascal, C, C&#xff0c;到后来的 java, c#,perl,php,再到现在流行的python。小时候的我总觉得多掌握一门语言&#xff0c;我的技术能力就又前…

【WC2014】时空穿梭【组合数】【莫比乌斯反演】【整除分块】【暴力多项式】

题意&#xff1a;TTT 组数据&#xff0c;给一个 nnn 维空间&#xff0c;第 iii 维大小为 [1,mi]∩Z[1,m_i]\cap \Z[1,mi​]∩Z&#xff0c;求大小为 ccc 的严格偏序上升的共线点集个数。答案模 100071000710007。 T≤100,n≤11,m≤105,c≤20T\leq 100,n\leq 11,m\leq 10^5,c\le…

2021牛客暑期多校训练营1 H Hash Function FFT\NTT

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个数组aaa&#xff0c;你需要找一个最小的模数xxx&#xff0c;使得aaa中每个数都模上xxx之后互不相同。 n≤5e5,ai≤5e5,ai!ajn\le5e5,a_i\le5e5,a_i!a_jn≤5e5,ai​≤5e5,ai​!aj​ 思路&#xff1a…

架构杂谈《六》

超时处理模式在服务化或者微服务架构里&#xff0c;传统的整体应用拆分成多个职责单一的微服务&#xff0c;微服务之间通过某种网络通信协议互相通信和交互&#xff0c;完成特定的功能&#xff0c;然而由于网络通信的不稳定&#xff0c;在设计系统时必须考虑到对网络通信的容错…

【BZOJ4543】Hotel加强版【神仙树形dp】【长链剖分】

题意&#xff1a;给一棵 nnn 个点的树&#xff0c;求两两距离相等的三元组个数。 n≤105n\leq 10^5n≤105 显然相当于是找一个点到这三个点距离相等。子树内和子树外到当前点的距离为某个值的点的个数可以长链剖分快速得到&#xff0c;但统计答案非常棘手。 接下来是个鬼才想…

hdu 6962 I love tree 线段树维护二次函数

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个点的一颗树&#xff0c;有mmm次询问&#xff0c;每次询问有两个操作&#xff1a; (1)(1)(1)将[a,b][a,b][a,b]路径上的点依次加上12,22,32,...,len2,lenpath(a,b)1^2,2^2,3^2,...,len^2,lenpath(a,…

基于surging 的stage组件设计,谈谈我眼中的微服务

一、前言surging 开源地址&#xff1a;https://github.com/dotnetcore/surging随着业务的发展&#xff0c;并发量的增多&#xff0c;业务的复杂度越来越大&#xff0c;对于系统架构能力要求越来越高&#xff0c;这时候微服务的设计思想应运而生&#xff0c;但是对于微服务需要引…

【PKUSC2018】星际穿越【结论】【倍增dp】

题意&#xff1a;有一张边权为 111 的无向图&#xff0c;对 i∈[2,n]i\in [2,n]i∈[2,n]&#xff0c;iii 与 [li,i−1][l_i,i-1][li​,i−1] 间有边。 qqq 次询问 l,r,xl,r,xl,r,x&#xff0c;表示 xxx 与 [l,r][l,r][l,r] 中的所有点的最短路长度的平均值&#xff0c;其中 l&l…

HDU - 6971 K - I love max and multiply sosdp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 直接求i&j>ki\And j>ki&j>k不是很好求&#xff0c;所以转换成i&jki\And jki&jk的情况。 考虑对a,ba,ba,b求一遍超集&#xff0c;让后从[0,n−1][0,n-1][0,n−1]扫…