bzoj3990[SDOI2015]排序

http://www.lydsy.com/JudgeOnline/problem.php?id=3990

DFS

好吧,表示不会做。

发现对于这些搜索的题我比较弱,看来需要加强一下。

回归正题。

我们发现对于一个操作方案(不妨记操作数为$cnt$),我们任意改变操作的顺序,总可以满足条件。

根据最小表示法的原理,我们规定按照编号从小到大进行操作,如果可行,那么$ans+=cnt!$

假设我们做到第$i$个操作,此时对于第$1$到第$i-1$个操作,我们已经知道了各个操作是否用到;并且我们保证:如果序列划分成$2^{N-i+1}$个长度为$2^{i-1}$的段,那么每一段的数字都是递增且连续的。

那么我们在做第做到第$i$个操作时,我们要保证:如果序列划分成$2^{N-i}$个长度为$2^{i}$的段,那么每一段的数字都是递增且连续的。

怎么保证呢?

我们将序列划分成$2^{N-i}$个长度为$2^{i}$的段,看看每一段是否是递增且连续的。

如果有多于2个长度为$2^{i}$的段不是递增且连续的,那么肯定是没救的。

如果有1个长度为$2^{i}$的段不是递增且连续的,我们一定能将这个长度为$2^{i}$的段平均分成2个长度为$2^{i-1}$的段,并且这2个长度为$2^{i-1}$的段一定是递增且连续的(为什么?因为我们保证如果序列划分成$2^{N-i+1}$个长度为$2^{i-1}$的段,那么每一段的数字都是递增且连续的

那么交换这两段并判断是否可行,继续搜索。

如果有2个长度为$2^{i}$的段不是递增且连续的,类似得,我们分别将2个长度为$2^{i}$的段平均分成2个长度为$2^{i-1}$的段,总共有4个长度为$2^{i-1}$的段,有4种交换方法,交换并判断是否可行后继续搜索即可。

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<fstream>
#include<algorithm>
#include<cstring>
#include<string>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<utility>
#include<set>
#include<bitset>
#include<vector>
#include<functional>
#include<deque>
#include<cctype>
#include<climits>
#include<complex>
//#include<bits/stdc++.h>适用于CF,UOJ,但不适用于poj

using namespace std;

typedef long long LL;
typedef double DB;
typedef pair<int,int> PII;
typedef complex<DB> CP;

#define mmst(a,v) memset(a,v,sizeof(a))
#define mmcy(a,b) memcpy(a,b,sizeof(a))
#define fill(a,l,r,v) fill(a+l,a+r+1,v)
#define re(i,a,b)  for(i=(a);i<=(b);i++)
#define red(i,a,b) for(i=(a);i>=(b);i--)
#define ire(i,x) for(typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define fi first
#define se second
#define m_p(a,b) make_pair(a,b)
#define p_b(a) push_back(a)
#define SF scanf
#define PF printf
#define two(k) (1<<(k))

template<class T>inline T sqr(T x){return x*x;}
template<class T>inline void upmin(T &t,T tmp){if(t>tmp)t=tmp;}
template<class T>inline void upmax(T &t,T tmp){if(t<tmp)t=tmp;}

const DB EPS=1e-9;
inline int sgn(DB x){if(abs(x)<EPS)return 0;return(x>0)?1:-1;}
const DB Pi=acos(-1.0);

inline int gint()
  {
        int res=0;bool neg=0;char z;
        for(z=getchar();z!=EOF && z!=‘-‘ && !isdigit(z);z=getchar());
        if(z==EOF)return 0;
        if(z==‘-‘){neg=1;z=getchar();}
        for(;z!=EOF && isdigit(z);res=res*10+z-‘0‘,z=getchar());
        return (neg)?-res:res;
    }
inline LL gll()
  {
      LL res=0;bool neg=0;char z;
        for(z=getchar();z!=EOF && z!=‘-‘ && !isdigit(z);z=getchar());
        if(z==EOF)return 0;
        if(z==‘-‘){neg=1;z=getchar();}
        for(;z!=EOF && isdigit(z);res=res*10+z-‘0‘,z=getchar());
        return (neg)?-res:res;
    }

const int maxN=12;

int N;
int fac[maxN+10];
vector<int>A;
int ans;

inline int check(vector<int> &A,int l,int r)
  {
      int i;
      re(i,l+1,r)if(A[i]!=A[i-1]+1)return 0;
      return 1;
  }
inline void SWAP(vector<int> &A,int x,int y,int len)
  {
      int i;
      re(i,1,len)swap(A[x+i-1],A[y+i-1]);
  }
inline void DFS(vector<int> A,int k,int cnt)
  {
      if(k==N){ans+=fac[cnt];return;}
      int i,b[5],tot=0;
      for(i=1;i<=two(N);i+=two(k+1))
        if(!check(A,i,i+two(k+1)-1))
          {
              if(tot==4)return;
              b[++tot]=i,b[++tot]=i+two(k);
          }
        if(tot==0)
          {
              DFS(A,k+1,cnt);
              return;
          }
        vector<int>B;
        if(tot==2)
          {
              if(A[b[2]]+two(k)==A[b[1]])
                {
                    B=A;
                    SWAP(B,b[1],b[2],two(k));
                    DFS(B,k+1,cnt+1);
                }
              return;
          }
        if(tot==4)
          {
              if(A[b[3]]+two(k)==A[b[2]] && A[b[1]]+two(k)==A[b[4]])
                {
                      B=A;
                        SWAP(B,b[1],b[3],two(k));
                        DFS(B,k+1,cnt+1);
                }
              if(A[b[4]]+two(k)==A[b[2]] && A[b[3]]+two(k)==A[b[1]])
                {
                    B=A;
                    SWAP(B,b[1],b[4],two(k));
                    DFS(B,k+1,cnt+1);
                }
              if(A[b[1]]+two(k)==A[b[3]] && A[b[2]]+two(k)==A[b[4]])
                {
                    B=A;
                    SWAP(B,b[2],b[3],two(k));
                    DFS(B,k+1,cnt+1);
                }
              if(A[b[1]]+two(k)==A[b[4]] && A[b[3]]+two(k)==A[b[2]])
                {
                    B=A;
                    SWAP(B,b[2],b[4],two(k));
                    DFS(B,k+1,cnt+1);
                }
              return;
          }
    }

int main()
  {
      freopen("sort.in","r",stdin);
      freopen("sort.out","w",stdout);
      int i;
      N=gint();
      fac[0]=1;re(i,1,N)fac[i]=fac[i-1]*i;
      A.resize(two(N)+1);re(i,1,two(N))A[i]=gint();
      DFS(A,0,0);
      cout<<ans<<endl;
      return 0;
  }

时间: 2024-10-05 04:58:59

bzoj3990[SDOI2015]排序的相关文章

BZOJ3990 SDOI2015 排序 DFS

题意:给定一个长度为2^N的序列和N个操作,每个操作i为将2^N分为2^(N-i+1)段,然后任意交换其中两段,求有多少种不同的交换方案使得序列升序 题解: 由于一个合法的方案中,交换操作的先后顺序,方案依然合法,所以我们只需要确定使用哪些操作. 按i的大小从小到大枚举每一个操作i,然后将序列分为2^(N-i)段,看有多少段不是递增的.如果有两段以上显然无解,如果有一段那么直接折半交换,如果有两段则枚举四种情况看是否合法. #include <cstdio> #include <cstr

【BZOJ 3990】 [SDOI2015]排序

3990: [SDOI2015]排序 Time Limit: 20 Sec Memory Limit: 128 MB Submit: 152 Solved: 83 [Submit][Status][Discuss] Description 小A有一个1-2^N的排列A[1..2^N],他希望将A数组从小到大排序,小A可以执行的操作有N种,每种操作最多可以执行一次,对于所有的i(1<=i<=N),第i中操作为将序列从左到右划分为2^{N-i+1}段,每段恰好包括2^{i-1}个数,然后整体交换其

[SDOI2015]排序 题解 (搜索)

Description 小A有一个1-2^N的排列A[1..2^N],他希望将A数组从小到大排序,小A可以执行的操作有N种,每种操作最多可以执行一次,对于所有的i(1<=i<=N),第i中操作为将序列从左到右划分为2^{N-i+1}段,每段恰好包括2^{i-1}个数,然后整体交换其中两段.小A想知道可以将数组A从小到大排序的不同的操作序列有多少个,小A认为两个操作序列不同,当且仅当操作个数不同,或者至少一个操作不同(种类不同或者操作位置不同). 下面是一个操作事例: N=3,A[1..8]=[

BZOJ 3990 Sdoi2015 排序 DFS

题目大意:给定一个长度为2^n的排列,有n个操作,第i个操作为[将序列分成2^(n-i+1)段,每段长2^(i-1),然后任选两段交换],每个操作最多用一次,求有多少操作序列能把序列排出来 Orz dzy 首先我们很容易发现一个操作序列是否合法与序列的顺序是无关的 因此我们只需要确定某个操作序列中每个操作选不选就行了 那么这类操作序列对答案的贡献就是选择的操作数的阶乘 我们从小到大DFS,对于第i次操作我们将序列分成2^(n-i)段,每段长度2^i 我们找到序列中不是连续递增的段,如果这样的段超

luoguP3322 [SDOI2015]排序

首先我们可以容易地知道任意区间交换的顺序对答案没有影响. 所以我们可以按照区间的长度进行搜索. 又因为每一种长度的区间只能交换一次,所以我们可以进行剪枝. 对于当前搜索区间的长度\(2^x\),我们可以对于每一个长度为\(2^{x+1}\)的区间,判断它是不是单调递增且相邻两数之间差\(1\),如果不是,则打上标记. 如果被标记了的区间超过了\(2\),那么无论如何都不能排好序了,直接回溯(想一想,为什么). 下面我们分标记区间的个数情况讨论: \(0\)个:直接进行下一层的搜索. \(1\)个

【搜索】BZOJ 3990: 【Sdoi 2015】排序

3990: [SDOI2015]排序 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 336  Solved: 164[Submit][Status][Discuss] Description 小A有一个1-2^N的排列A[1..2^N],他希望将A数组从小到大排序,小A可以执行的操作有N种,每种操作最多可以执行一次,对于所有的i(1<=i<=N),第i中操作为将序列从左到右划分为2^{N-i+1}段,每段恰好包括2^{i-1}个数,然后整体交换

【SDOI2015】bzoj3990 排序

A. 排序 题目描述 输入格式 输出格式 一行,一个整数,表示可以将数组A从小到大排序的不同的操作序列的个数. 样例 样例输入 3 7 8 5 6 1 2 4 3 样例输出 6 数据范围与提示 对于30%的数据,1<=N<=4: 对于全部的数据,1<=N<=12. 一个伪装成计数题的大搜索… 首先要知道一个结论:对于一个操作序列,如果他是合法的,那么他的全排列都是合法的,所以从小到大搜索,同时记录操作数k. 对于第x次操作,把序列分成2x-1段,暴力搜索有几段不满足a[i]=a[i

BZOJ3990 排序

题目:www.lydsy.com/JudgeOnline/problem.php?id=3990 这题很不错. 刚开始时无从下手,想了好多$O((2^n)log(2^n))$ 的idea,但是都不行. 后来去看题解发现操作序列是满足交换率的,然后竟然是搜索. 因为swap是swap的逆运算(歪歪的) 然后只要从小到大枚举操作序列就可以了. 这样类似分治下去,当你在计算长度为$2^i$的序列时已经保证了所有长度为$2^{i-1}$的序列的「连续且递增」. 注意是「连续且递增」,开始W了好多发,然后

【SDOI2015第1轮第1试】排序

#include<cstdio> #include<iostream> using namespace std; typedef long long LL; const int N = 1 << 13; int n , a[N]; LL ans , fac[13]; inline bool check(int len) { for(register int i = 1; i <= 1 << n - len; i++) if (a[(i - 1) * (