CodeForces 272 D.Dima and Two Sequences（组合数学）

Description

给出两个长度为$n$的点对$(a_1,1),(a_2,2),...,(a_n,n)$和$(b_1,1),(b_2,2),...,(b_n,n)$，要求用这两个序列构成一个长度为$2n$的点对序列，使得$x$坐标非减，问方案数

Input

第一行一整数$n$，之后输入两个长度为$n$的序列$a_i$和$b_i$，最后输入模数$m(1\le n\le 10^5,1\le a_i,b_i\le 10^9,2\le m\le 10^9+7)$

Output

输出方案数，结果模$m$

Sample Input

1
1
2
7

Sample Output

1

Solution

假设$a,b$序列中值为$x$的数字有$cnt_x$个，那么这$cnt_x$个数字对应的点对必然在一段连续区间，其所有排列均可，但是有重复的情况，即当$a_i=b_i$时，相当于有两个相同的点对$(a_i,i)$，这会被算两次，假设所有相同点对有$num$个，那么答案即为$\frac{\prod\limits_{x}cnt_x!}{2^{num}}$ ，由于模数不一定是素数，故要在算阶乘的时候把这$num$个 2 <script type="math/tex" id="MathJax-Element-105">2</script>除掉

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=100001;
map<int,int>M,vis;
map<int,int>::iterator it;
int n,a[maxn],b[maxn],m;
int main()
{
    while(~scanf("%d",&n))
    {
        M.clear();
        for(int i=1;i<=n;i++)scanf("%d",&a[i]),M[a[i]]++;
        for(int i=1;i<=n;i++)scanf("%d",&b[i]),M[b[i]]++;
        scanf("%d",&m);
        int num=0;
        for(int i=1;i<=n;i++)if(a[i]==b[i])num++;
        int ans=1;
        for(it=M.begin();it!=M.end();it++)
        {
            for(int i=2;i<=it->second;i++)
            {
                int j=i;
                while(num&&j%2==0)num--,j/=2;   
                ans=(ll)ans*j%m;
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}