码迷,mamicode.com
首页 > 其他好文 > 详细

counting swaps

时间:2018-06-16 21:09:21      阅读:178      评论:0      收藏:0      [点我收藏+]

标签:white   namespace   格式   sizeof   layer   span   vol   fir   ref   

3602 Counting Swaps 0x30「数学知识」例题

背景

https://ipsc.ksp.sk/2016/real/problems/c.html

Just like yesterday (in problem U of the practice session), Bob is busy, so Alice keeps on playing some single-player games and puzzles. In her newest puzzle she has a permutation of numbers from 1 to n. The goal of the puzzle is to sort the permutation using the smallest possible number of swaps.

Instead of simply solving the puzzle, Alice is wondering about the probability of winning it just by playing at random. In order to answer this question, she needs to know the number of optimal solutions to her puzzle.

描述

给定一个 1~n 的排列 p_1,p_2,…,p_n,可进行若干次操作,每次选择两个整数 x,y,交换 p_x,p_y。设把 p_1,p_2,…,p_n 变成单调递增的排列 1,2,…,n 至少需要 m 次交换。求有多少种操作方法可以只用 m 次交换达到上述目标。因为结果可能很大,你只需要输出对 10^9+9 取模之后的值。1≤n≤10^5。
例如排列 2,3,1 至少需要2次交换才能变为 1,2,3。操作方法共有3种,分别是:
先交换数字2,3,变成 3,2,1,再交换数字3,1,变成 1,2,3。
先交换数字2,1,变成 1,3,2,再交换数字3,2,变成 1,2,3。
先交换数字3,1,变成 2,1,3,再交换数字2,1,变成 1,2,3。

You are given a permutation p1,?…,?pn of the numbers 1 through n. In each step you can choose two numbers x?<?y and swap px with py.

Let m be the minimum number of such swaps needed to sort the given permutation. Compute the number of different sequences of exactly m swaps that sort the given permutation. Since this number may be large, compute it modulo 109?+?9.

输入格式

The first line of the input file contains an integer t specifying the number of test cases. Each test case is preceded by a blank line.

Each test case consists of two lines. The first line contains the integer n. The second line contains the sequence p1,?…,?pn: a permutation of 1,?…,?n.

In the easy subproblem C11?≤?n?≤?10.

In the hard subproblem C21?≤?n?≤?105.

输出格式

For each test case, output a single line with a single integer: x mod(10^9+9), where x is the number of ways to sort the given sequence using as few swaps as possible.

样例输入

3

3
2 3 1

4
2 1 4 3

2
1 2

样例输出

3
2
1

样例解释

In the first test case, we can sort the permutation in two swaps. We can make the first swap arbitrarily; for each of them, there’s exactly one optimal second swap. For example, one of the three shortest solutions is “swap p1 with p2 and then swap p1 with p3”.

In the second test case, the optimal solution involves swapping p1 with p2 and swapping p3 with p4. We can do these two swaps in either order.

The third sequence is already sorted. The optimal number of swaps is 0, and thus the only optimal solution is an empty sequence of swaps.

 

求长度为n的环变成自环的公式

f[n]=n^(n-2)

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int maxn=100000+10;
const long long mod=1000000009;

long long a[maxn];
long long num[maxn];
long long l[maxn];
bool v[maxn];

long long pow(long long x,long long k){
    long long ans=1%mod;
    for (;k;k>>=1){
        if(k&1) ans=ans*x%mod;
        x=(long long)x*x%mod;
    }
    return ans%mod;
}

int main(){
    //freopen("data.in","r",stdin);
    //freopen("data.out","w",stdout);
    a[0]=1;
    for(int i=1;i<=100000;i++) a[i]=a[i-1]*i%mod;
    int t;
    scanf("%d",&t);
    while(t--){
        memset(v,0,sizeof(v));
        int n;
        scanf("%d",&n);
        for (int i=1;i<=n;i++) scanf("%lld",&num[i]);
        int top=0;
        for (int i=1;i<=n;i++){
            if(!v[i]){
                int pp=1;
                v[i]=true;
                for (int j=num[i];j!=i;j=num[j]){
                    ++pp;
                    v[j]=true;
                }
                l[++top]=pp;
            }
        }
        long long ans=1;
        for (int i=1;i<=top;i++){
            ans=ans*(l[i]==1 ? 1 : pow(l[i],l[i]-2))%mod;
            ans=ans*pow(a[l[i]-1],mod-2)%mod;
        }
        ans=ans*a[n-top]%mod;
        printf("%lld\n",ans%mod);
    }
return 0;
}

 

counting swaps

标签:white   namespace   格式   sizeof   layer   span   vol   fir   ref   

原文地址:https://www.cnblogs.com/lmjer/p/9191341.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!