【JZOJ 省选模拟】向日葵人生

Description

Input

Output

Sample Input

样例 1 输入:
2
3 2
1 2
2 3

Sample Output

样例 1 输出:
665496241

Data Constraint

Hint

思路

树上一条路径可行的概率是1/len.（len个点，第一个点要第一个选到）
环上做一个容斥：路径A可行+路径B可行-路径AB均可行。
仙人掌上做dp容斥即可。

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int mod = 998244353,N = 477;
int n,m;
ll fac[N],unfac[N],inv[N];
ll power(ll x,ll y)
{ll ret = 1;for(; y; y >>= 1){if (y & 1) ret = ret * x % mod;x = x * x % mod;}return ret;
}int len[N];
ll f[N][N];
int no[N],cir[N],cnt,cirbel[N];
int ls[N],nx[N * 8],to[N * 8],tot = 1;void add(int x,int y) {to[++tot] = y,nx[tot] = ls[x],ls[x] = tot;
}int dfn[N],stm,p[N],e[N * 8];
int cirfa[N],dep[N];
void dfs(int x) {dfn[x] = ++ stm;for(int i = ls[x]; i; i = nx[i]) {if (i == p[x]) continue;int y = to[i];if (dfn[y] == 0) {p[y] = i ^ 1;dep[y] = dep[x] + 1;dfs(y);} else if(dfn[y] < dfn[x]) {e[i] = e[i ^ 1] = 1;++cnt;for(int z = x,w = 1; z != y; z = to[p[z]],w++) {no[z] = w;len[z] = dep[x] - dep[y] + 1;add(z,y),add(y,z);e[p[z]] = e[p[z] ^ 1] = 1;cirfa[z] = y;cirbel[z] = cnt;}}}
}ll ans;
void work(int x,int from) {for(int i = 1; i <= n; i ++) {ans = (ans + inv[i] * f[x][i]) % mod;}for(int i = ls[x]; i; i = nx[i]) if (e[i] == 0) {int y = to[i];if (y == from) continue;if (y != cirfa[x] && cirfa[y] != x) { for(int i = 1; i <= n; i++) f[y][i] = f[x][i - 1];work(y,x);} else {int ori = 0,a = 0,b = 0,cirlen = 0;if (cirfa[x] == y) {ori = x,a = no[x],b = len[x] - no[x],cirlen = len[x];} else {if (cirbel[y] == cirbel[from]) {cirlen = len[y];ori = from,a = (no[from] - no[y] + len[y]) % cirlen,b = cirlen - a;} else {ori = x,a = no[y],b = len[y] - no[y],cirlen = len[y];}}for(int i = 1; i <= n; i++) {if (i - a > 0) f[y][i] = f[ori][i - a];if (i - b > 0) f[y][i] = (f[y][i] + f[ori][i - b]) % mod;if (i - cirlen + 1 > 0) f[y][i] = (f[y][i] - f[ori][i - cirlen + 1]) % mod;}work(y,x);}}
}int main() {freopen("falldream.in","r",stdin);freopen("falldream.out","w",stdout);cin >> n >> n >> m;for(int i = 1; i <= m; i++) {int x,y; scanf("%d %d",&x,&y);add(x,y),add(y,x);}dfs(1);fac[0] = 1;for(int i = 1; i <= n; i++) fac[i] = fac[i - 1] * i % mod;unfac[n] = power(fac[n],mod - 2);for(int i = n - 1; ~i; i--) unfac[i] = unfac[i + 1] * (i + 1) % mod;for(int i = 1; i <= n; i++) inv[i] = unfac[i] * fac[i - 1] % mod;for(int r = 1; r <= n; r++){memset(f,0,sizeof f);f[r][1] = 1;work(r,0);}cout << (ans + mod) % mod << endl;
}

【JZOJ 省选模拟】向日葵人生相关推荐

【JZOJ 省选模拟】人生(life)
题目 Description Input 第一行一个正整数n. 第二行n个整数,其中第i个整数为1表示点i的颜色为黑色,为0表示点i的颜色为白色,为-1 则不确定. Output 输出一行一个整数表示 ...
6687. 【JZOJ 省选模拟】树没了(tree)
题目 Description Input Output 2操作的个数行,每行代表一个询问的答案. Sample Input 5 5 1 1 2 3 4 5 1 2 2 3 2 4 1 5 Q 1 M ...
【JZOJ 省选模拟】多项式（poly）
题目 Description Input Output Sample Input 样例输入 1: 2 7 -7 14 样例输入 2: 3 1 4 1 5 样例输入 3: 0 998244353 Sam ...
【JZOJ 省选模拟】铺路（road ）
题目 Description Input Output Sample Input 样例输入 1: 4 4 1 3 4 2 4 8 1 2 2 3 4 3 样例输入 2: 3 2 1 2 3 2 3 4 ...
【JZOJ 省选模拟】我的朋友们
题目 Description Input Output 一行一个整数,表示题目描述中的期望对998244353取模的值. Sample Input 样例 1输入: 4 3 1 5 4 6 3 7 1 ...
【JZOJ 省选模拟】感受清风
题目 Description Input Output 对于每个 qu 事件或者 qd 事件输出一个整数表示答案,答案对1000001取模(提示:这个数不是质数,有101×9901=1000001) ...
【JZOJ 省选模拟】6681.图
题目 Description 给定一张无重边.自环.割点的平面图,你需要回答 Q 次询问,每次询问会给出一个简单环,你需要回答在由这个简单环围成的多边形内部(包括边界上)的点有多少个. 保证图中每条线 ...
【JZOJ 省选模拟】6707.异或树（xor）
题目 Description Input Output Sample Input 样例输入1: 2 2 998244353 样例输入2: 10 3 998244353 Sample Output 样例 ...
【JZOJ 省选模拟】迷宫（maze）
题目 Description Input Output Sample Input 样例输入 4 8 0 A 1 B 1 A 0 B 2 A 3 B 3 A 2 B 0 C 2 D 2 C 0 D 1 ...

【JZOJ 省选模拟】向日葵人生

Description

Input

Output

Sample Input

Sample Output

Data Constraint

Hint

思路

代码

【JZOJ 省选模拟】向日葵人生相关推荐

最新文章

热门文章