排列 (搜索)

排序

内存限制:128 MiB 时间限制:1000 ms 标准输入输出
 
 
题目描述

image

输入格式

image

输出格式

一行,一个整数,表示可以将数组A从小到大排序的不同的操作序列的个数。

样例
样例输入
3
7 8 5 6 1 2 4 3
样例输出
6
数据范围与提示

对于30%的数据,1<=N<=4; 对于全部的数据,1<=N<=12。

一群数论题中出了一个搜索。

甚至我考试时一点也不会。

自己算复杂度比较大然后交上去只跑了90ms。

这也太

不多说了。

引理1

当前可以将原数列变成排序的操作,其每一种其他排列都可以将原序列变成排序。

证明 不太好证,举几个例子

例如 3 4 2 1这个序列可以经过操作1 操作2 变成 1 2 3 4

那么显然 先进行操作2 再进行操作1 也可以变成 1 2 3 4

于是推而广之 (数学归纳)用cnt次操作可以达到目的,那么所有符合的方案数为cnt!

这样通过引理1 可以剪很多枝

我们从小到大枚举每一种方案,每次验证符合,若符合ans+cnt!就完了

引理2

假设枚举到第i个方案 设划分序列为区间长度为1<<(i+1)的序列,如果有三个及以上不符合递增且相连的区间,那么这个序列没救了。

这个我们可以验证

首先我们保证如果我们进行到第i个方案那么所有长度1<<i的序列都符合递增且相连

如果有一个长度1<<(i+1)不是递增且相连的区间,我们把它划分成两个1<<i,i<<i的序列,我们直接用第i种操作(可能)就行了

然后如果有两个长度那么把它划分成四个,然后也(可能)行

但如果有三个及以上我们最多保证经过i操作后,两个序列符合条件

所以得证

然后搜就完了

以下是本人丑陋的代码

 1 #include<bits/stdc++.h>
 2 #define ll long long
 3 #define A 1000000
 4 using namespace std;
 5 ll xx[A],a[A],b[A],ans,n;
 6 ll read()
 7 {
 8         ll x=0,f=1;char ch=getchar();
 9         while(!isdigit(ch)){
10                 if(ch=='-') f=-1;
11                 ch=getchar();
12         }
13         while(isdigit(ch)){
14                 x=x*10+ch-'0';
15                 ch=getchar();
16         }
17         return f*x;
18 }
19 void SWAP(ll x,ll k,ll len)
20 {
21         for(ll i=x,j=k;len;len--,i++,j++) 
22                 swap(a[i],a[j]);
23 }
24 bool check(ll l,ll r)
25 {
26         for(ll i=l;i<r;i++)
27                 if(a[i]+1!=a[i+1])
28                         return 0;
29         return 1;
30 }
31 void shuchu()
32 {
33         for(ll i=1;i<=(1<<n);i++)
34                 cout<<a[i]<<"";
35         cout<<endl;
36 }
37 
38 void dfs(ll k,ll cnt)
39 {
40 //        shuchu();
41 //        cout<<k<<endl;
42         if(k==n)
43         {
44                 ans+=xx[cnt];
45                 return ;
46         }
47         ll now=(1<<(k+1)),x=(1<<k),b[5],tot=0;
48         for(ll i=1;i<=(1<<n);i+=now)
49         {       
50 //                printf("检查i=%lld j=%lld\n",i,i+now-1); 
51                 if(!check(i,i+now-1))
52                 {
53                         b[++tot]=i;
54                         if(tot==3) return ;
55                 }
56         }
57 //        cout<<"tot="<<tot<<"x="<<x<<"now="<<now<<endl;
58 //        for(ll i=1;i<=n;i++)
59 //        {
60 //                
61 //        }
62         if(tot==0) {dfs(k+1,cnt);return ;}
63         if(tot==1)
64         {
65                 SWAP(b[1],b[1]+x,x);
66                 if(check(b[1],b[1]+x)) dfs(k+1,cnt+1);
67                 SWAP(b[1],b[1]+x,x);
68         }
69         if(tot==2)
70         {
71                 SWAP(b[1],b[2],x);
72                 if(check(b[1],b[1]+x)&&check(b[2],b[2]+x)) dfs(k+1,cnt+1);
73                 SWAP(b[1],b[2],x);
74                 SWAP(b[1]+x,b[2],x);
75                 if(check(b[1],b[1]+x)&&check(b[2],b[2]+x)) dfs(k+1,cnt+1);
76                 SWAP(b[1]+x,b[2],x);
77                 SWAP(b[1],b[2]+x,x);
78                 if(check(b[1],b[1]+x)&&check(b[2],b[2]+x)) dfs(k+1,cnt+1);
79                 SWAP(b[1],b[2]+x,x);
80                 SWAP(b[1]+x,b[2]+x,x);
81                 if(check(b[1],b[1]+x)&&check(b[2],b[2]+x)) dfs(k+1,cnt+1);
82                 SWAP(b[1]+x,b[2]+x,x);
83         }
84 //        shuchu();
85 }
86 int main()
87 {
88         n=read();xx[0]=1;
89         for(ll i=1;i<=(1<<n);i++) a[i]=read();
90         for(ll i=1;i<=n;i++) xx[i]=xx[i-1]*i;
91         dfs(0,0);
92         cout<<ans<<endl;
93 }
View Code

 

转载于:https://www.cnblogs.com/znsbc-13/p/11157192.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值