HDU6038 - Function

---

做法：
展開后有，\(f(i) = b_{f(a_i)} = b_{b_{f(a_{a_i})}} = ... = b_{b_{..b_{f(i)}}}\)，可以發現當 \(a_i\) 所在的循環節中確定一個時，整個循環都確定了，根據這個式子還可以發現，對于一個\(a_i\)所在的大小為\(x\)的循環節，一定映射到一個\(b_i\)中的一個大小為\(y\)循環節，且\(y|x\) 。因此答案就是：
\[ \prod_{i=1}^k \sum_{d|l_i} d·B_d \]
（其中\(a\)有\(k\)個大小為\(l_i\)的循環節，扎著\(b\)中大小為\(i\)的循環節有 \(B_i\) 個）

#include <bits/stdc++.h> typedef long long ll; const int N = 1e5 + 7; const int mod = 1e9 + 7; using namespace std; int n, m, a[N], b[N], in[N], CC = 0; ll B[N]; ll cal(int st, int a[]) {int now = st; ll ans = 0;while(!in[now]) {in[now] = 1; now = a[now];++ans;}return ans; } int main() {while(scanf("%d%d",&n,&m)!=EOF) {for(int i = 0; i < n; ++i) scanf("%d", &a[i]);for(int i = 0; i < m; ++i) scanf("%d", &b[i]);memset(B,0,sizeof(B));memset(in,0,sizeof(in));for(int i = 0; i < m; ++i) if(!in[i]) ++B[cal(i,b)];memset(in,0,sizeof(in));ll ans = 1LL;for(int i = 0; i < n; ++i) if(!in[i]) {ll x = cal(i,a), tmp = 0;for(ll j = 1; j*j <= x; ++j) if(x%j == 0) {(tmp += j*B[j]%mod)%=mod;if(j*j==x) continue;(tmp += ((x/j)*B[x/j])%mod)%=mod;}(ans *= tmp)%=mod;}printf("Case #%d: %lld\n",++CC,ans);} }

轉載于:https://www.cnblogs.com/RRRR-wys/p/9699781.html

總結

以上是生活随笔為你收集整理的HDU6038 - Function的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。