C.简单环
题目链接:https://ac.nowcoder.com/acm/contest/114/C
题目描述
给定一张n个点m条边的无向图,求出图中所有简单环的数量。(简单环:简单环又称简单回路,图的顶点序列中,除了第一个顶点和最后一个顶点相同外,其余顶点不重复出现的回路叫简单回路。或者说,若通路或回路不重复地包含相同的边,则它是简单的)
n<=20
解题思路
n很小,所以可以状态压缩枚举,设dp[status][i]表示当前枚举的链为集合statucs(注意当前是一条链,枚举的子集不能当成环)
以i这个点为结束点,以status第一个开始点为起始点的方案数,转移的时候枚举在子集中的点j,不在status子集中的点k,看点k是否可以加入这条链,dp[status|(1<<k)][k]+=dp[status][j],然后再看k是否可以和status的起点建边,如果可以,这就是一个合法的简单环方案,加入答案。
#include<bits/stdc++.h>
typedef long long ll;
const int maxn = 1e5 + 10;
const ll mod = 998244353;
using namespace std;
int vis[25][25];
ll dp[1 << 20][20], tot[1 << 20], ans[20];
int main() {
ios::sync_with_stdio(0);
int n, m, k;
ll inv = (mod + 1) / 2;
while(cin >> n >> m >> k) {
memset(vis, 0, sizeof vis);
memset(dp, 0, sizeof dp);
memset(tot, 0, sizeof tot);
memset(ans, 0, sizeof ans);
for(int i = 0; i < m; i++) {
int x, y; cin >> x >> y;
x--; y--;
vis[x][y] = vis[y][x] = 1;
}
for(int i = 0; i < n; i++) dp[1 << i][i] = 1;//单独的一个点也算一种方案
int max_sta = 1 << n;
for(int i = 0; i < max_sta; i++) {
int min_point = -1;
for(int j = 0; j < n; j++) if(((i >> j) & 1)) {
if(min_point == -1) min_point = j; ///固定以标号最小的点为链的开始点
for(int k = min_point + 1; k < n; k++) {
if((i >> k) & 1) continue;
if(!vis[j][k]) continue;
int next_sta = i | (1 << k);
(dp[next_sta][k] += dp[i][j]) %= mod;
if(vis[min_point][k]) (tot[next_sta] += dp[i][j]) %= mod;
}
}
int tx = __builtin_popcount(i);
if(tx >= 3) (ans[tx % k] += tot[i] * inv) %= mod; ///环计算了两次
}
for(int i = 0; i < k; i++) cout << ans[i] << endl;
}
return 0;
}
E.跳格子
链接:https://ac.nowcoder.com/acm/contest/114/E
来源:牛客网
题目描述
sum 个格子排成一排,每次可以往前跳 1-n 格,往回跳 1-m 格,而且在往回跳的时候只能跳在往前跳时踩过的格子。
现在,在格子上跳,问跳到最后一个格子上最后又跳回第一个格子之前的方案数 mod 233333333。
注意:只能一直向前跳,跳到最后一个格子,然后再往回跳,跳回第一个格子之前。
解题思路
跳过去又调回来,可以看成从出发点跳两次,而第二次跳的集合必须是第一次的子集。
所以先求出第一次跳过去的,这很简单,记为a[i],a[i]+=a[i-j];
再求第二次跳过去的,因为第二次跳过去的每个点都必须是第一次跳过的点,所以
b[i]+=b[i-j]*a[j],这样转移就好了.
#include <bits/stdc++.h>
#define int long long
using namespace std;
const int maxn=4e5+10;
const int mod=233333333;
int a[maxn];
int b[maxn];
signed main(){
int T;
scanf("%lld",&T);
while(T--){
int cnt,n,m;
scanf("%lld%lld%lld",&cnt,&n,&m);
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
a[0]=1;
for(int i=1;i<=cnt;i++){
for(int j=1;j<=i&&j<=n;j++){
(a[i]+=a[i-j])%=mod;
}
}
b[0]=1;
for(int i=1;i<=cnt;i++){
for(int j=1;j<=i&&j<=m;j++){
(b[i]+=b[i-j]*a[j])%=mod;
}
}
printf("%lld\n",b[cnt]);
}
return 0;
}