HDU 6038 Function (数学)
来源:互联网 发布:教师培训课程 大数据 编辑:程序博客网 时间:2024/06/07 13:31
Description
You are given a permutation a from 0 to n−1 and a permutation b from 0 to m−1.
Define that the domain of function f is the set of integers from 0 to n−1, and the range of it is the set of integers from 0 to m−1.
Please calculate the quantity of different functions f satisfying that f(i)=bf(ai) for each i from 0 to n−1.
Two functions are different if and only if there exists at least one integer from 0 to n−1 mapped into different integers in these two functions.
The answer may be too large, so please output it in modulo 10^9+7.
Input
The input contains multiple test cases.
For each case:
The first line contains two numbers n, m. (1≤n≤100000,1≤m≤100000)
The second line contains n numbers, ranged from 0 to n−1, the i-th number of which represents ai−1.
The third line contains m numbers, ranged from 0 to m−1, the i-th number of which represents bi−1.
It is guaranteed that ∑n≤106, ∑m≤106.
Output
For each test case, output “Case #x: y” in one line (without quotes), where x indicates the case number starting from 1 and y denotes the answer of corresponding case.
Sample Input
3 21 0 20 13 42 0 10 2 3 1
Sample Output
Case #1: 4Case #2: 4
题意
给出两个序列 a 和 b ,求满足
思路
根据样例 1 我们可以得到:
代换以后有
考虑置换
于是
而如果
因此满足条件的函数个数就是
AC 代码
#include<iostream>#include<vector>#include<cstring>#include<algorithm>#include<cmath>#include<cstdio>using namespace std;typedef __int64 LL;const int maxn=110000;const int mod = 1e9+7;int a[maxn],b[maxn];vector<int> na,nb;bool vis[maxn];void findnum(int *a,int n,vector<int> &res) //寻找 res 中的循环节个数以及长度{ memset(vis,false,sizeof(vis)); for(int i=0; i<n; i++) { if(!vis[i]) { int now=a[i],len=0; while(!vis[now]) { ++len; vis[now]=true; now=a[now]; } res.push_back(len); } }}int main(){ ios::sync_with_stdio(false); int n,m,ti=0; while(cin>>n>>m) { na.clear(); nb.clear(); for(int i=0; i<n; i++) cin>>a[i]; for(int i=0; i<m; i++) cin>>b[i]; findnum(a,n,na); findnum(b,m,nb); LL ans=1; int lena=na.size(); int lenb=nb.size(); for(int i=0; i<lena; i++) { LL res=0; for(int j=0; j<lenb; j++) { if(na[i]%nb[j]==0) res=(res+nb[j])%mod; } ans=(ans*res)%mod; } cout<<"Case #"<<++ti<<": "<<ans<<endl; } return 0;}
- HDU 6038 Function (数学)
- 【HDU 6038】Function(图+组合数学)
- HDU 6038 Function(组合数学)
- HDU 3903 Trigonometric Function(数学定理)
- HDU 4498 Function Curve(数学)
- 【HDU 6050 Funny Function】 + 数学
- Function (HDU 6038)
- HDU-6038 Function(思维)
- hdoj function 5608 (数学)
- HDU 6038 Function (强连通+计数)
- hdu 6038 Function (思维题搞一搞)
- HDU 6038(Function) 思维 Java
- HDU 6038 Function(强连通分量)
- HDU 6038 Function(找循环节)
- HDU 6038 Function(置换->循环节)
- HDU 6038 Function
- HDU-6038 Function
- hdu 6038 Function(思维)
- 取数游戏
- DOM基础一
- Java结合postgreSql的增删改查遇到问题2
- 《android手势密码》
- 如何证明我妈是我妈?如今终于有了答案!
- HDU 6038 Function (数学)
- Spring AOP : AspectJ Pointcut 切点
- 剑指offer解题报告(Java版)——二叉搜索树转换为双向链表 27
- WPF C# 多屏情况下,实现窗体显示到指定的屏幕内
- linux命令的理解
- C# 之 类复制 MemberwiseClone与Clone(深 浅 Clone)
- 手把手教你写基于C++ Winsock的图片下载的网络爬虫
- 剑指offer题解C++【3】
- 戴尔DELL 笔记本 改WIN7 BIOS 设置详解