标签：DP，拓扑排序

题目

题目传送门

题目描述

「恒逸，你相信灵魂的存在吗？」郭恒逸和姚枫茜漫步在枫音乡的街道上。望着漫天飞舞的红枫，枫茜突然问出这样一个问题。

「相信吧。不然我们是什么，一团肉吗？要不是有灵魂…我们也不可能再见到你姐姐吧。」恒逸给出了一个略微无厘头的回答。枫茜听后笑了笑。「那你仔细观察过枫叶吗？」说罢，枫茜伸手，接住了一片飘落的枫叶。

「其实每一片枫叶都是有灵魂的。你看，枫叶上不是有这么多脉络吗？我听说，枫叶上有一些特殊的位置，就和人的穴位一样。脉络都是连接在这些穴位之间的。枫树的灵魂流过每片枫叶的根部，沿着这些脉络，慢慢漫进穴位，沁入整片枫叶。也是因为这个原因，脉络才都是单向的，灵魂可不能倒着溜回来呢。」恒逸似懂非懂地点了点头。枫茜接着说了下去。

「正是因为有了灵魂，每片枫叶才会与众不同。也正是因为有了灵魂，每片枫叶也都神似其源本的枫树，就连脉络也形成了一棵树的样子。但如果仔细看的话，会发现，在脉络树之外，还存在其它的非常细的脉络。虽然这些脉络并不在树上，但他们的方向也同样顺着灵魂流淌的方向，绝不会出现可能使灵魂倒流的回路。」恒逸好像突然想到了什么。「那这些脉络岂不是可以取代已有的脉络，出现在脉络树上？」枫茜闭上了眼睛。

「是啊，就是这样。脉络树并不是唯一的。只要有一些微小的偏差，脉络树就可能差之万里，哪怕是在这同一片枫叶上。就像我们的故事，结局也不是唯一的。只要改变一个小小的选项，故事流程可能就会被彻底扭转。」

「真是深奥啊…」恒逸盯着这片红枫，若有所思地说。枫茜继续说道。

「还不止如此呢。所有的脉络都不会永恒存在，也不会永恒消失。不管是脉络树上的脉络，还是之外的细小脉络，都是如此。存在的脉络可能断开消失，消失的脉络也可能再次连接。万物皆处在永恒的变化之中，人与人之间的羁绊也是。或许有一天，我们与大家的羁绊也会如同脉络一样，被无情地斩断。或许我们也终将成为”枫音乡的过客“。或许这一切都会是必然，是枫树的灵魂所决定的…」

枫茜的眼角泛起了几滴晶莹剔透的泪珠。恒逸看着这样的枫茜，将她抱入怀中。

「别这样想，枫茜。就算脉络断开，也有可能还会有新的脉络树，也还会与枫树的根相连。这样的话，我们的羁绊仍然存在，只是稍微绕了一些远路而已。无论如何，我都不会离开你的。因为你是我穷尽一生所寻找的，我的真恋啊！」

两人的目光对上了。枫茜幸福地笑了，把头埋进了恒逸的怀抱。从远方山上的枫林中，传来了枫的声音。

【问题描述】不妨假设枫叶上有 n个穴位，穴位的编号为 1 ~ n。有若干条有向的脉络连接着这些穴位。穴位和脉络组成一个有向无环图——称之为脉络图（例如图 1），穴位的编号使得穴位 1 没有从其他穴位连向它的脉络，即穴位 1 只有连出去的脉络；由上面的故事可知，这个有向无环图存在一个树形子图，它是以穴位 1为根的包含全部n个穴位的一棵树——称之为脉络树（例如图 2和图 3给出的树都是图1给出的脉络图的子图）；值得注意的是，脉络图中的脉络树方案可能有多种可能性，例如图2和图 3就是图 1给出的脉络图的两个脉络树方案。

脉络树的形式化定义为：以穴位 r 为根的脉络树由枫叶上全部 n个穴位以及 n- 1 条脉络组成，脉络树里没有环，亦不存在从一个穴位连向自身的脉络，且对于枫叶上的每个穴位 s，都存在一条唯一的包含于脉络树内的脉络路径，使得从穴位r 出发沿着这条路径可以到达穴位 s。现在向脉络图添加一条与已有脉络不同的脉络（注意：连接 2个穴位但方向不同的脉络是不同的脉络，例如从穴位3到4的脉络与从4到3的脉络是不同的脉络，因此，图 1 中不能添加从 3 到 4 的脉络，但可添加从 4 到 3 的脉络），这条新脉络可以是从一个穴位连向自身的（例如，图 1 中可添加从 4 到 4 的脉络）。原脉络图添加这条新脉络后得到的新脉络图可能会出现脉络构成的环。请你求出添加了这一条脉络之后的新脉络图的以穴位 1 为根的脉络树方案数。

由于方案可能有太多太多，请输出方案数对 1,000,000,007 取模得到的结果。

输入输出格式

输入格式

输入文件的第一行包含四个整数 n、m、x和y，依次代表枫叶上的穴位数、脉络数，以及要添加的脉络是从穴位 x连向穴位y的。接下来 m行，每行两个整数，由空格隔开，代表一条脉络。第 i 行的两个整数为ui和vi，代表第 i 条脉络是从穴位 ui连向穴位vi的。

输出格式

输出一行，为添加了从穴位 x连向穴位 y的脉络后，枫叶上以穴位 1 为根的脉络树的方案数对 1,000,000,007取模得到的结果。

输入输出样例

输入样例#1

输出样例#1

说明

对于所有测试数据，1≤n≤100000,n−1≤m≤min(200000,n(n−1)2),1≤x,y,ui,vi≤n1 \leq n \leq 100000, n - 1 \leq m \leq min(200000, \frac {n(n -1)} {2}),1 \leq x, y, u_i, v_i \leq n1≤n≤100000,n−1≤m≤min(200000,2n(n−1)),1≤x,y,ui,vi≤n

分析

对于原图来说（DAG）

ans=∏i=2ndegreeians=\prod_{i=2}^ndegree_ians=∏i=2ndegreei

新加入一条边，就有可能出现环

把原答案减去环的情况就可以了

del=∑uinloop∑u−>vf[v](vinloop)degreeudel=\sum_{u\ in\ loop} \frac {\sum_{u->v} f[v] (v\ in\ loop)} {degree_u}del=∑u in loopdegreeu∑u−>vf[v](v in loop)

注意要特判y=1y=1y=1的情况

code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define dep(i,a,b) for(int i=a;i>=b;i--)
#define ll long long
#define mem(x,num) memset(x,num,sizeof x)
#define reg(x) for(int i=last[x];i;i=e[i].next)using namespace std;
inline ll read(){ll f=1,x=0;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;
}
//******head by yjjr******
const int maxn=1e5+6,mod=1e9+7;
ll ans=1;
int n,m,x,y,cnt=0,last[maxn],du[maxn],b[maxn],que[maxn<<2];
ll f[maxn],inv[maxn<<1];
struct edge{int to,next;}e[maxn<<1];
void insert(int u,int v){e[++cnt]=(edge){v,last[u]};last[u]=cnt;}
void dp(){f[y]=ans;int head=0,tail=1;rep(i,1,n)if(!du[i])que[tail++]=i;while(head<tail){int now=que[head++];f[now]=f[now]*inv[b[now]]%mod;reg(now){int v=e[i].to;         f[v]=(f[v]+f[now])%mod;du[v]--;if(!du[v])que[tail++]=v;}}
}
int main(){n=read(),m=read(),x=read(),y=read();inv[1]=1;rep(i,2,m+1)inv[i]=(-inv[mod%i]*(mod/i)%mod+mod)%mod;rep(i,1,m){int u=read(),v=read();insert(u,v);du[v]++;}du[y]++;rep(i,1,n)b[i]=du[i];rep(i,2,n)ans=ans*du[i]%mod;du[y]--;if(y==1){cout<<ans<<endl;return 0;}dp();(ans-=f[x]-mod)%=mod;cout<<ans<<endl;return 0;
}

洛谷3244 [HNOI2015]落忆枫音相关推荐

P3244 [HNOI2015]落忆枫音（DAG上的动态规划问题，朱刘定理，乘法逆元）
P3244 [HNOI2015]落忆枫音样例太坑了!竟然和题目描述给的图不一样! 题目描述给定一张有向无环图,这张图满足一个性质:以点1为根节点,保证至少有一棵有向树,连接所有的节点. 现在向这张 ...
[HNOI2015] 落忆枫音
题目描述「恒逸,你相信灵魂的存在吗?」郭恒逸和姚枫茜漫步在枫音乡的街道上.望着漫天飞舞的红枫,枫茜突然问出这样一个问题. 「相信吧.不然我们是什么,一团肉吗?要不是有灵魂......我们也不可能再 ...
Bzoj4011 [HNOI2015]落忆枫音
Time Limit: 10 Sec Memory Limit: 512 MB Submit: 983 Solved: 533 Description 「恒逸,你相信灵魂的存在吗?」郭恒逸和姚枫 ...
BZOJ4011:[HNOI2015]落忆枫音(DP,拓扑排序)
Description 「恒逸,你相信灵魂的存在吗?」郭恒逸和姚枫茜漫步在枫音乡的街道上.望着漫天飞舞的红枫,枫茜突然问出这样一个问题. 「相信吧.不然我们是什么,一团肉吗?要不是有灵魂--我们也 ...
洛谷 P3244 / loj 2115 [HNOI2015] 落忆枫音题解【拓扑排序】【组合】【逆元】
组合计数的一道好题.什么非主流题目题目背景 (背景冗长请到题目页面查看) 题目描述不妨假设枫叶上有 \(n\) 个穴位,穴位的编号为 \(1\sim n\).有若干条有向的脉络连接着这些穴位. ...
[HNOI2015]落忆枫音
题目描述不妨假设枫叶上有 n个穴位,穴位的编号为 1 ~ n.有若干条有向的脉络连接着这些穴位.穴位和脉络组成一个有向无环图--称之为脉络图(例如图 1),穴位的编号使得穴位 1 没有从其他穴 ...
【题解】 [HNOI2015]落忆枫音（拓扑排序+dp+容斥原理）
原题戳我 Solution: (部分复制Navi_Aswon博客) 解释博客中的两个小地方: \[\sum_{\left(S是G中y→x的一条路径的点集\right))}\prod_{2≤j≤n,(j ...
luogu3244 bzoj4011 HNOI2015 落忆枫音
这道题目题面真长,废话一堆. 另外:这大概是我第一道独立做出来的HNOI2011年以后的题目了吧.像我水平这么差的都能做出来,dalao您不妨试一下自己想想? 题目大意:给一个DAG,其中1号点没有入 ...
luogu P3244 [HNOI2015]落忆枫音
传送门 md这题和矩阵树定理没半毛钱关系qwq 首先先不考虑有环,一个\(DAG\)个外向树个数为\(\prod_{i=2}^{n}idg_i(\)就是\(indegree_i)\),因为外向树每个点 ...
BZOJ 4011 HNOI2015 落忆枫音
AC通道:http://www.lydsy.com/JudgeOnline/problem.php?id=4011 题目很长,写得也很有诗意与浪漫色彩,让我们不禁感叹出题人是一个多么英俊潇洒的人. 所 ...

洛谷3244 [HNOI2015]落忆枫音

题目