\(\mathcal{Description}\)

  Link.

  给定一个含 \(n\) 个结点 \(m\) 条边的简单无向图,每条边的边权是一个常数项为 \(0\) 的 \(T\) 次多项式,求所有从 \(1\) 结点出发回到 \(1\) 结点的环路中,边权之积的 \(T\) 次项系数和。

  \(n,m\le10\),\(T\le4\times10^4\)。

\(\mathcal{Solution}\)

  令 \(f_i(x)=\sum_{j\ge0}f_{i,j}x^j\),从 \(1\) 出发到 \(i\) 的所有路径边权积的和,那么对于一条边 \(e(u,v)\),设其边权为 \(g_e\),它会为 \(f\) 提供转移:

\[f_{u,i}\leftarrow f_{u,i}+\sum_{d=1}^T [x^d]g_ef_{v,i-d}
\]

\[f_u(x)\leftarrow f_u(x)+g_e(x)f_v(x)
\]

所以整体做一个分治 FFT 就能求出所有 \(f\)。复杂度 \(\mathcal O(mT\log^2T)\)。

\(\mathcal{Code}\)

/*~Rainybunny~*/

#include <cstdio>
#include <vector>
#include <cassert>
#include <algorithm> #define rep( i, l, r ) for ( int i = l, rep##i = r; i <= rep##i; ++i )
#define per( i, r, l ) for ( int i = r, per##i = l; i >= per##i; --i ) typedef std::vector<int> Poly; const int MAXN = 10, MAXL = 1 << 17, MOD = 998244353;
int N, M, T, eu[MAXN + 5], ev[MAXN + 5];
Poly E[MAXN + 5], F[MAXN + 5]; inline void subeq( int& a, const int b ) { ( a -= b ) < 0 && ( a += MOD ); }
inline int sub( int a, const int b ) { return ( a -= b ) < 0 ? a + MOD : a; }
inline int add( int a, const int b ) { return ( a += b ) < MOD ? a : a - MOD; }
inline void addeq( int& a, const int b ) { ( a += b ) >= MOD && ( a -= MOD ); }
inline int mul( const int a, const int b ) { return int( 1ll * a * b % MOD ); }
inline int mpow( int a, int b ) {
int ret = 1;
for ( ; b; a = mul( a, a ), b >>= 1 ) ret = mul( ret, b & 1 ? a : 1 );
return ret;
} namespace PolyOper { const int MG = 3;
int omega[17][MAXL]; inline void init() {
rep ( i, 0, 16 ) {
int* oi = omega[i];
oi[0] = 1, oi[1] = mpow( MG, MOD - 1 >> i >> 1 );
rep ( j, 2, ( 1 << i ) - 1 ) oi[j] = mul( oi[j - 1], oi[1] );
}
} inline void ntt( Poly& u, const int type ) {
static int rev[MAXL]; rev[0] = 0;
int n = int( u.size() ), lgn = 1; for ( ; 1 << lgn < n; ++lgn );
rep ( i, 1, n - 1 ) rev[i] = rev[i >> 1] >> 1 | ( i & 1 ) << lgn >> 1;
rep ( i, 1, n - 1 ) if ( i < rev[i] ) {
u[i] ^= u[rev[i]] ^= u[i] ^= u[rev[i]];
}
for ( int i = 0, stp = 1; stp < n; ++i, stp <<= 1 ) {
int* oi = omega[i];
for ( int j = 0; j < n; j += stp << 1 ) {
rep ( k, j, j + stp - 1 ) {
int ev = u[k], ov = mul( oi[k - j], u[k + stp] );
u[k] = add( ev, ov ), u[k + stp] = sub( ev, ov );
}
}
}
if ( !~type ) {
int ivn = MOD - ( MOD - 1 ) / n;
rep ( i, 0, n - 1 ) u[i] = mul( u[i], ivn );
std::reverse( u.begin() + 1, u.end() );
}
} } // namespace PolyOper. inline Poly operator * ( Poly u, Poly v ) {
assert( u.size() && v.size() );
int su = int( u.size() ), sv = int( v.size() ), len = 1;
for ( ; len < su + sv - 1; len <<= 1 ); u.resize( len ), v.resize( len );
PolyOper::ntt( u, 1 ), PolyOper::ntt( v, 1 );
rep ( i, 0, len - 1 ) u[i] = mul( u[i], v[i] );
PolyOper::ntt( u, -1 ); return u.resize( su + sv - 1 ), u;
} inline void solve( const int l, const int r ) {
if ( l == r ) return ;
int mid = l + r >> 1; solve( l, mid ); rep ( i, 1, M ) {
int u = eu[i], v = ev[i];
static Poly A, B, R; A = { F[u].begin() + l, F[u].begin() + mid + 1 };
B = { E[i].begin() + 1, E[i].begin() + r - l + 1 };
R = A * B;
rep ( j, mid + 1, r ) addeq( F[v][j], R[j - l - 1] ); A = { F[v].begin() + l, F[v].begin() + mid + 1 };
B = { E[i].begin() + 1, E[i].begin() + r - l + 1 };
R = A * B;
rep ( j, mid + 1, r ) addeq( F[u][j], R[j - l - 1] );
} solve( mid + 1, r );
} int main() {
PolyOper::init();
scanf( "%d %d %d", &N, &M, &T );
rep ( i, 1, M ) {
scanf( "%d %d", &eu[i], &ev[i] ), --eu[i], --ev[i];
E[i].resize( T + 1 );
rep ( j, 1, T ) scanf( "%d", &E[i][j] );
} rep ( i, 0, N - 1 ) F[i].resize( T + 1 );
F[0][0] = 1, solve( 0, T );
printf( "%d\n", F[0][T] );
return 0;
}

Solution -「ABC 213H」Stroll的更多相关文章

  1. Solution -「ABC 219H」Candles

    \(\mathcal{Description}\)   Link.   有 \(n\) 支蜡烛,第 \(i\) 支的坐标为 \(x_i\),初始长度为 \(a_i\),每单位时间燃烧变短 \(1\) ...

  2. Solution -「ABC 215H」Cabbage Master

    \(\mathcal{Description}\)   Link.   有 \(n\) 种颜色的,第 \(i\) 种有 \(a_i\) 个,任意两球互不相同.还有 \(m\) 个盒子,每个盒子可以被放 ...

  3. Solution -「ABC 213G」Connectivity 2

    \(\mathcal{Description}\)   Link.   给定简单无向图 \(G=(V,E)\),点的编号从 \(1\) 到 \(|V|=n\).对于 \(k=2..n\),求 \(H= ...

  4. Solution -「ABC 217」题解

    D - Cutting Woods 记录每一个切割点,每次求前驱后驱就好了,注意简单判断一下开闭区间. 考场上采用的 FHQ_Treap 无脑莽. #include <cstdio> #i ...

  5. Solution -「ARC 104E」Random LIS

    \(\mathcal{Description}\)   Link.   给定整数序列 \(\{a_n\}\),对于整数序列 \(\{b_n\}\),\(b_i\) 在 \([1,a_i]\) 中等概率 ...

  6. Solution Set -「ABC 217」

      大家好屑兔子又来啦! [A - Lexicographic Order]   说个笑话,\(\color{black}{\text{W}}\color{red}{\text{alkingDead} ...

  7. Solution -「ARC 110E」Shorten ABC

    \(\mathcal{Description}\)   Link.   给定长度为 \(n\),包含 A, B, C 三种字符的字符串 \(S\),定义一次操作为将其中相邻两个不相同的字符替换为字符集 ...

  8. Solution -「CTS 2019」「洛谷 P5404」氪金手游

    \(\mathcal{Description}\)   Link.   有 \(n\) 张卡牌,第 \(i\) 张的权值 \(w_i\in\{1,2,3\}\),且取值为 \(k\) 的概率正比于 \ ...

  9. Solution -「BZOJ 3812」主旋律

    \(\mathcal{Description}\)   Link.   给定含 \(n\) 个点 \(m\) 条边的简单有向图 \(G=(V,E)\),求 \(H=(V,E'\subseteq E)\ ...

随机推荐

  1. vs2017 winform 组件 -- 总结

    1.ComboBox  [下拉框] (1) 添加选项 this.[控件名].Items.Add("内容") (2)设置下拉框 自动完成 模式 和 数据源 this.[控件名].Au ...

  2. tomcat启动卡在了 At least one JAR was scanned for TLDs yet contained no TLDs 的根本原因与解决办法

    1.前言 有时候服务器开启时启动不了,卡在了 org.apache.jasper.servlet.TldScanner.scanJars At least one JAR was scanned fo ...

  3. nuxt写路由接口

    //在server/interface/city.js import Router from 'koa-router'; const router = new Router({ prefix:'/ci ...

  4. MapReduce和Hive学习文档链接学习顺序

    1.<CentOS6.5下安装Hadoop-2.7.3(图解教程)> https://www.toutiao.com/i6627365258090512909/ 2.<CentOS6 ...

  5. docker安装easymock

    一.准备 参考docker安装mongodb与redis文章 二.拉取 docker pull docker.io/easymock/easymock 三.启动 1.创建配置目录 mkdir -f / ...

  6. docker镜像制作及发布

    以centos为例. 主要内容:安装docker,制作镜像,发布镜像. 安装docker 1. 安装docker yum install -y docker 等待一会,安装成功,查看安装列表 2. 启 ...

  7. Windows和Linux关闭占用端口

    关闭端口的方式有很多种,但是常用的就是这种比较来的快一点 如果通过以下方式解决不了,可以通过关闭服务来解决 Windows 1.查看端口占用的进程 netstat -ano | findstr 800 ...

  8. 程序员必备的编程助手!SmartCoder助你轻松集成HMS Core

    当开发者在集成HMS Core遇到一些疑问时,需要翻阅官网文档,反复查看集成说明或者API调用说明,或者研究GitHub上的开源示例代码,花费较多的时间,在IDE环境和网页浏览器之间反复切换也会耗费很 ...

  9. collate utf8_bin是什么意思

    创建数据库时 collate utf_bin是 以二进制值比较,也就是区分大小写,collate是核对的意思 uft-8_general_ci  一般比较,不区分大小写

  10. JavaScript之详述闭包导致的内存泄露

    一.内存泄露 1. 定义:一块被分配的内存既不能使用,也不能回收.从而影响性能,甚至导致程序崩溃. 2. 起因:JavaScript的垃圾自动回收机制会按一定的策略找出那些不再继续使用的变量,释放其占 ...