告别

题目大意

一个n的排列，每次可以选择三个不同位置（这三个位置有序，(2,1,3)与(1,2,3)不同），然后对它们进行一次轮换。
只能操作m次，某一次操作过后当前排列变成了目标排列则停止操作。
每次均等概率选择三个有序位置进行操作。
在模意义下求m次操作内变成目标排列的概率。

设状态

首先可以扭转初始排列和目标排列，使得目标排列变成有序的（就是变成1，2，3……n）
我们知道可以用置换来表示一个排列，目标排列此时对应n个单位置换，同时n个单位置换也一定只对应目标排列。
对于两种置换，如果其中所有置换长度组成的多重集合完全相同，能拓展出的状态也完全相同，因此应该考虑缩到一起。
这样子状态只有135种。
暴力枚举每种状态，然后可以构造出与其对应的一个可行排列，再暴力枚举三个有序位置，然后再找会变成哪个状态（可以二分或map）
由于走到目标状态就停止操作，因此目标状态不能有连向其他状态的边，只能有一个自环（方便我们得到m次以内而不是恰好m次）
得到转移矩阵后可以做矩阵乘法。

#include<cstdio>#include<algorithm>#include<map>#include<vector>#define fo(i,a,b) for(i=a;i<=b;i++)using namespace std;const int maxn=150+10,mo=998244353;typedef long long ll;vector<int> now;map<vector<int>,int> ha;int c[maxn][20];int dis[maxn][maxn],ans[maxn][maxn],o[maxn][maxn],sta[60],a[20],b[20],d[20],ts[20],p[20];bool bz[20];int i,j,k,l,t,n,m,sid,tid,tot,top,sum,cnt,w1,w2,w3,num;int quicksortmi(int x,int y){    if (!y) return 1;    int t=quicksortmi(x,y/2);    t=(ll)t*t%mo;    if (y%2) t=(ll)t*x%mo;    return t;}void dfs(int x,int dow,int y){    if (!y){        int i;        top++;        fo(i,1,x-1) c[top][i]=p[i];        c[top][0]=x-1;        now.clear();        fo(i,1,x-1) now.push_back(p[i]);        ha[now]=top;        return;    }    int i;    fo(i,dow,y){        p[x]=i;        dfs(x+1,i,y-i);    }}int getid(){    int i;    now.clear();    fo(i,1,cnt) now.push_back(d[i]);    return ha[now];}int change(int x,int y,int z){    int i;    fo(i,1,n) b[i]=a[i];    int t=b[x];    b[x]=b[y];    b[y]=b[z];    b[z]=t;    cnt=0;    fo(i,1,n) bz[i]=0;    fo(i,1,n)        if (!bz[i]){            t=0;            j=i;            while (!bz[j]){                bz[j]=1;                t++;                j=b[j];            }            d[++cnt]=t;        }    sort(d+1,d+cnt+1);    return getid();}int main(){    freopen("goodbye.in","r",stdin);freopen("goodbye.out","w",stdout);    scanf("%d%d",&n,&m);    num=quicksortmi(n*(n-1)*(n-2),m);    num=quicksortmi(num,mo-2);    fo(i,1,n) scanf("%d",&a[i]);    fo(i,1,n) scanf("%d",&b[i]);    fo(i,1,n) ts[b[i]]=i;    fo(i,1,n) a[i]=ts[a[i]],b[i]=i;    dfs(1,1,n);    fo(i,1,n)        if (!bz[i]){            t=0;            j=i;            while (!bz[j]){                bz[j]=1;                t++;                j=a[j];            }            d[++cnt]=t;        }    sort(d+1,d+cnt+1);    sid=getid();    fo(i,1,n) d[i]=1;    cnt=n;    tid=getid();    fo(i,1,top){        j=1;k=1;        while (j<=c[i][0]){            fo(l,k,k+c[i][j]-2) a[l]=l+1;            a[k+c[i][j]-1]=k;            k+=c[i][j];            j++;        }        fo(w1,1,n)            fo(w2,1,n)                fo(w3,1,n)                      if (w1!=w2&&w1!=w3&&w2!=w3){                        j=change(w1,w2,w3);                        if (i==1) j=1;                        (dis[i][j]+=1)%=mo;                    }    }    //dis[tid][n+1]=dis[n+1][n+1]=1;    /*fo(i,1,n){        fo(j,1,n)            printf("%d ",dis[i][j]);        printf("\n");    }*/    //m++;    while (m){        sta[++sum]=m%2;        m/=2;    }    fo(i,1,top) ans[i][i]=1;    while (sum){        fo(i,1,top)            fo(j,1,top)                o[i][j]=0;        fo(k,1,top)            fo(i,1,top)                fo(j,1,top)                    o[i][j]=(o[i][j]+(ll)ans[i][k]*ans[k][j]%mo)%mo;        fo(i,1,top)            fo(j,1,top)                ans[i][j]=o[i][j];        if (sta[sum]){            fo(i,1,top)                fo(j,1,top)                    o[i][j]=0;            fo(k,1,top)                fo(i,1,top)                    fo(j,1,top)                        o[i][j]=(o[i][j]+(ll)ans[i][k]*dis[k][j]%mo)%mo;            fo(i,1,top)                fo(j,1,top)                    ans[i][j]=o[i][j];        }        sum--;    }    /*fo(i,1,n){        fo(j,1,n)            printf("%d ",ans[i][j]);        printf("\n");    }*/    fo(i,1,top)        fo(j,1,top)            dis[i][j]=0;    dis[1][sid]=1;    fo(i,1,top)        fo(j,1,top)            o[i][j]=0;    fo(k,1,top)        fo(i,1,top)            fo(j,1,top)                o[i][j]=(o[i][j]+(ll)dis[i][k]*ans[k][j]%mo)%mo;    fo(i,1,top)        fo(j,1,top)            ans[i][j]=o[i][j];    printf("%d\n",(ll)ans[1][1]*num%mo);}