NC20116. [HNOI2015]落忆枫音
描述
「恒逸,你相信灵魂的存在吗?」 郭恒逸和姚枫茜漫步在枫音乡的街道上。望着漫天飞舞的红枫,枫茜突然问出这样一个问题。
「相信吧。不然我们是什么,一团肉吗?要不是有灵魂......我们也不可能再见到你姐姐吧。」 恒逸给出了一个略微无厘头的回答。枫茜听后笑了笑。 「那你仔细观察过枫叶吗?」 说罢,枫茜伸手,接住了一片飘落的枫叶。
「其实每一片枫叶都是有灵魂的。你看,枫叶上不是有这么多脉络吗?我听说,枫叶上有一些特殊的位置,就和人的穴位一样。脉络都是连接在这些穴位之间的。枫树的灵魂流过每片枫叶的根部,沿着这些脉络,慢慢漫进穴位,沁入整片枫叶。也是因为这个原因,脉络才都是单向的,灵魂可不能倒着溜回来呢。」 恒逸似懂非懂地点了点头。枫茜接着说了下去。
「正是因为有了灵魂,每片枫叶才会与众不同。也正是因为有了灵魂,每片枫叶也都神似其源本的枫树,就连脉络也形成了一棵树的样子。但如果仔细看的话,会发现,在脉络树之外,还存在其它的非常细的脉络。虽然这些脉络并不在树上,但他们的方向也同样顺着灵魂流淌的方向,绝不会出现可能使灵魂倒流的回路。」 恒逸好像突然想到了什么。 「那这些脉络岂不是可以取代已有的脉络,出现在脉络树上?」 枫茜闭上了眼睛。
「是啊,就是这样。脉络树并不是唯一的。只要有一些微小的偏差,脉络树就可能差之万里,哪怕是在这同一片枫叶上。就像我们的故事,结局也不是唯一的。只要改变一个小小的选项,故事流程可能就会被彻底扭转。」
「真是深奥啊......」 恒逸盯着这片红枫,若有所思地说。枫茜继续说道。
「还不止如此呢。所有的脉络都不会永恒存在,也不会永恒消失。不管是脉络树上的脉络,还是之外的细小脉络,都是如此。存在的脉络可能断开消失,消失的脉络也可能再次连接。万物皆处在永恒的变化之中,人与人之间的羁绊也是。或许有一天,我们与大家的羁绊也会如同脉络一样,被无情地斩断。或许我们也终将成为”枫音乡的过客“。或许这一切都会是必然,是枫树的灵魂所决定的......」
枫茜的眼角泛起了几滴晶莹剔透的泪珠。恒逸看着这样的枫茜,将她抱入怀中。
「别这样想,枫茜。就算脉络断开,也有可能还会有新的脉络树,也还会与枫树的根相连。这样的话,我们的羁绊仍然存在,只是稍微绕了一些远路而已。无论如何,我都不会离开你的。因为你是我穷尽一生所寻找的,我的真恋啊!」
两人的目光对上了。枫茜幸福地笑了,把头埋进了恒逸的怀抱。从远方山上的枫林中,传来了枫的声音。
脉络树的形式化定义为:以穴位 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
输入:
4 4 4 3 1 2 1 3 2 4 3 2
输出:
3
C++14(g++5.4) 解法, 执行用时: 284ms, 内存消耗: 9108K, 提交时间: 2020-03-29 12:52:42
#include <iostream> #include <algorithm> #include <cstring> using namespace std; typedef long long ll; const ll MAXN=400007; const ll MOD=1e9+7; ll n,m,x,y,rt[MAXN],rd[MAXN],start[MAXN],rs[MAXN],u[MAXN]; struct ed{ ll s,e; }ge[MAXN]; ll pw(ll a,int p) { if (p==0) return 1; ll ab=pw(a,p/2); ab=ab*ab%MOD; if (p%2==1) return ab*a%MOD; else return ab; } bool operator<(ed a,ed b) { return a.s<b.s; } void tps(ll x) { // cout<<x<<endl; for (ll i=start[x];ge[i].s==x;i++) { // cout<<ge[i].e<<' '<<rd[ge[i].e]<<endl; rt[ge[i].e]+=rt[x]*pw(u[ge[i].e],MOD-2);//*(start[x-1]-start[x]); rt[ge[i].e]%=MOD; rd[ge[i].e]--; if (rd[ge[i].e]==0) tps(ge[i].e); } } int main() { cin>>n>>m>>x>>y; for (ll i=1;i<=m;i++) { cin>>ge[i].s>>ge[i].e; } sort(ge+1,ge+m+1); for (ll i=1;i<=m;i++) { if (ge[i].s>ge[i-1].s) start[ge[i].s]=i; rd[ge[i].e]++; } memcpy(u,rd,sizeof(rd)); start[n+1]=m+1; ll ans=1; //rd[y]++; //rd[y]--; rt[y]=pw(rd[y]+1,MOD-2); tps(1);if (y!=1){ ans=(MOD+1-rt[x]); u[y]++;} else ans=1; for (int i=2;i<=n;i++){ ans*=u[i]; // cout<<u[i]<<' '; ans%=MOD; } cout<<ans<<endl; // for (int i=1;i<=n;i++) // cout<<i<<' '<<rt[i]<<endl; }
C++(clang++11) 解法, 执行用时: 100ms, 内存消耗: 6776K, 提交时间: 2021-05-20 08:53:49
#include <bits/stdc++.h> typedef long long ll; const int N=100005,mu=1e9+7; int inv[N],dp[N],x,y,n,m,s,t,du[N]; std::queue<int> q; std::vector<int> e[N]; void reduce(int &x){ x+=x>>31μ } int ksm(int x,int y){ int ans=1; for (;y;y>>=1,x=(ll)x*x%mu) if (y&1) ans=(ll)ans*x%mu; return ans; } int main(){ scanf("%d%d%d%d",&n,&m,&s,&t); if (n==1) return puts("1"),0; du[t]++; for (int i=1;i<=m;i++){ scanf("%d%d",&x,&y); du[y]++; e[x].push_back(y); } int ans=1; for (int i=2;i<=n;i++){ ans=(ll)ans*du[i]%mu; inv[i]=ksm(du[i],mu-2); } dp[t]=1; du[t]--; for (int i=1;i<=n;i++) if (!du[i]) q.push(i); while (q.size()){ int x=q.front(); //printf("%d\n",x); q.pop(); dp[x]=(ll)dp[x]*inv[x]%mu; //printf("%d\n",dp[x]); for (int u:e[x]){ reduce(dp[u]+=dp[x]-mu); du[u]--; if (!du[u]) q.push(u); } } //for (int i=1;i<=n;i++) printf("%d\n",dp[i]); printf("%d\n",(ll)ans*(mu+1-dp[s])%mu); }