uva11212_Editing a Book(A*算法)


///
作者:tt2767
声明:本文遵循以下协议自由转载-非商用-非衍生-保持署名|Creative Commons BY-NC-ND 3.0
查看本文更新与讨论请点击:http://blog.csdn.net/tt2767
链接被删请百度: CSDN tt2767
///


这题学到了很多啊,也碰到一些原来没有遇见过的情况。

题解: A*算法,移动一块最多可以导致3个后缀改变,如果改变都是我们所期望的,那么我们可以认为最少进行h/3次搜索就可以搜到了,如果 当前搜索步数+最少修改步数 > 最大步数,就不可能完成了,剪枝即可。

这题有意思的一点是h(),求相邻逆序数函数,我刚开始的时候是这样写的:

int h()//计算相邻逆序数对的个数
{
    int cnt = 0;
    for(int i = 0; i < n-1 ; i++)
        if(a[i]> a[i+1])
            cnt++;
    return cnt;
}

然后 数组最大长度N设为210,结果跑了2900ms才过,后来我又读了一边题,把N值改为了10, 也要1300ms,我百思不得其解,发现代码仓库是这样写的:

int h()
{
    int cnt = 0;
    for(int i = 1; i < n; i++)
        if(a[i] -1 != a[i-1])
            cnt++;
    return cnt;
}

理论上来说,都是扫了一遍,但为什么用 > 比用!= 慢了这么多呢?上网上查了之后也是不太明白,这真是让我十分不解。如果各位大神有明白这个问题的无比指点我一下。。。。

AC正解:


#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
//
#include<iostream>
#include<algorithm>
#include<string>
#include <iterator>
#include<sstream>
#include<functional>
#include<numeric>
///
#include<vector>
#include<map>
#include <stack>
#include<queue>
#include<set>
#include <bitset>
#include <list>
using namespace std;
#define lch(x)  ((x) << 1)
#define rch(x) ((x)<<1|1)
#define dad(x) ((x)>>1)
#define lowbit(x) ((x)&(-x))
typedef  long long int LL;
const int INF = ~0U>>1;
const double eps = 1e-6;
const long double PI = acos(0.0) * 2.0;
const int N = 10 ;

#define lx (j-i+1)    //片段x,y的长度
#define ly (k-j)

int a[N];
int ans,maxd;
int n;

int solve();
bool dfs(int d,int maxd);
int h(); //计算相邻逆序数对的个数
bool ordered(); //判断当前是否有序

int main()
{
    //ios::sync_with_stdio(false);
#ifdef ONLINE_JUDGE
#else
    freopen("in.txt", "r", stdin);
    //freopen("out.txt", "w", stdout);
#endif
    int tot = 1;
    while(scanf("%d",&n)==1&& n)
    {
        for(int  i =  0 ; i < n ; i++)    scanf("%d",&a[i]);
        printf("Case %d: %d\n",tot++,solve());
    }
    return 0;
}

int solve()
{
    if(ordered()) return 0;  // n<=9 时最多 5 次
    for(int maxd = 1 ; maxd < 5 ; maxd++)
        if(dfs(0,maxd))
            return maxd;
    return 5;
}

bool dfs(int d,int maxd)
{
    if(3*d + h() > 3*maxd )     return false;
    if(ordered())    return true;
    int cpy[N],tmp[N];              //tmp为缓存
    memcpy(cpy , a ,sizeof(a));//制作a的副本

    for(int i = 0 ; i < n ; i++)
    {
        if(i == 0 || cpy[i]-1 != cpy[i-1])  //尽量对连续的片段进行操作
        for(int  j = i ; j < n ; j++)
        {
            while(j+1 < n && cpy[j+1] == cpy[j]+1)  j++;
            memcpy(tmp,cpy+i , sizeof(int)*lx);        //将片段x放入缓存区

            for(int k = j+1 ; k < n ;k++)
            {
                while(k+1<n && cpy[k+1]==cpy[k]+1)  k++;
                memcpy(a+i,cpy+j+1,sizeof(int)*ly);    //将片段y放入片段x的位置
                memcpy(a+i+ly,tmp,sizeof(int)*lx);         //把x插入y的后面
                if(dfs(d+1,maxd))   return 1;
                memcpy(a,cpy,sizeof(a));
            }
        }
    }
    return 0;
}

int h()
{
    int cnt = 0;
    for(int i = 1; i < n; i++)
        if(a[i] -1 != a[i-1])
            cnt++;
    return cnt;
}

bool ordered()
{
    for(int i = 0 ; i < n ; i++)
        if(a[i] != i+1)
            return 0;
    return 1;
}

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

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值