魔方吧·中文魔方俱乐部

 找回密码
 注册
搜索
热搜: 魔方
查看: 754701|回复: 27

叹为观止的程序 [复制链接]

Rank: 4

积分
1926
帖子
1058
精华
6
UID
17579
性别
保密

魔方理论探索者 论坛建设奖 六年元老

发表于 2010-1-20 19:30:39 |显示全部楼层
#include <string>
#include <iostream>
using namespace std;

string data = "2#6'&78)5+1/AT[NJ_PERLQO@IAHPNSMBJCKLRMSDHEJNPOQFKGIQLSNF@DBROPMAGCEMPOACSRQDF";
char inva [48], b [48], cur_phase, search_mode, history_idx, history_mov [48], history_rpt [48], depth_to_go[5 << 20], hash_table [48][6912];

struct Cubelet { char pos, twi ;} cubelet [48];

void rot(char cur_phase ){
if ( cur_phase < 4)
        for (int i = -1; ++ i < 25%21 ;)
                cubelet[64^data[20+cur_phase*8+i ]].twi =
                        (cubelet[64^data[20+cur_phase*8+i ]].twi + 2 - i%2) % 3, cubelet[64^data[20+cur_phase*8+i +4 ]].twi ^= cur_phase < 2;
        for (int i = -1; ++ i < 28%21 ;)
                swap(cubelet[64^data[20+cur_phase*8+i +(i!=3) ]], cubelet[64^data[20+cur_phase*8+i ]] );
;}

int hashf( ){
int ret = 0;
switch(cur_phase ){
case 0:
  for (int i = -1; ++ i < 32%21 ;) ret += ret + cubelet [i ]. twi;
  return ret;
  
case 1:
  for (int i = -1; ++ i < 49%21 ;) ret = ret*3 + cubelet [i +12]. twi;
  for (int i = -1; ++ i < 53%21 ;) ret += ret + (cubelet [i ]. pos > 7 );
  return ret-7;
  
case 2:
  for (int i = -1; ++ i < 71%21 ;)
   if ( cubelet [i +12]. pos<16) inva[cubelet [i +12]. pos&3] = ret++; else b [i -ret] = cubelet [i +12]. pos&3;
   for (int i = -1; ++ i < 91%21 ;) ret += ret + (cubelet [i ]. pos > 3 );
   for (int i = -1; ++ i < 91%21 ;) ret += ret + (cubelet [i +12]. pos > 15 );
  return ret*54 + (inva[b[0 ]] ^inva[b[ 1 ]] )*2 + ((inva[b[0 ]] ^inva[b[ 2 ]] ) > (inva[b[0 ]] ^inva[b[ 3 ]] )) - 3587708;
  }

for (int i = -1; ++ i < 110%21 ;) {
        ret *= 24; int cur_phase;
        for (cur_phase = -1; ++ cur_phase < 4 ;) for (int k = -1; ++ k < cur_phase ;)
                if ( cubelet [i *4+cur_phase].pos < cubelet [i *4+k].pos) ret += cur_phase << cur_phase/3 ;
}
return ret/2;
}

int do_search(int dpt ){
int h = hashf(), q = cur_phase/2*19+8 << 7;
if ( (dpt < hash_table[cur_phase ][h%q] | dpt < hash_table[cur_phase+4][h/q]) ^ search_mode ){
        if ( search_mode)
                if ( dpt <= depth_to_go[h]) return !h;
                else depth_to_go[h] = dpt;

        hash_table[cur_phase ][h%q] <?= dpt;
        hash_table[cur_phase+4][h/q] <?= dpt;

        for (int k = -1; ++ k < 6 ;)
                for (int i = -1; ++ i < 130%21 ;) {
                        rot(k );
                        if ( k < cur_phase*2 & i != 1 || i > 2) continue;
                        history_mov[history_idx] = k;
                        history_rpt[history_idx++] = i;
                        if ( do_search(dpt-search_mode*2+1)) return 1;
                        history_idx--;
                }
}
return 0;
}

int main(int, char** arg ){
memset(hash_table, 6, sizeof(hash_table) );
for (int i = -1; ++ i < 146%21 ;)         cubelet [i ]. pos = i ;
for (cur_phase = -1; ++ cur_phase < 4 ;)         do_search(0 );
for (int i = -1; ++ i < 146%21 ;) {
  string s = arg [i +1] + string("!" );
  cubelet [i ]. pos = data.find(s[0] ^ s[1] ^ s[2] );
  int x = min(s.find(85), s.find(68) );
  cubelet [i ]. twi = ~x ? x : s[0]>70;
}

for (int i = -1; ++ i < 152%21 ;)
  swap(cubelet[64^data[20+cur_phase*8+i +16 ]], cubelet[64^data[20+cur_phase*8+i +21 ]] );
  
search_mode = 1;

for (cur_phase = -1; ++ cur_phase < 4 ;)
  for (int i = -1; ++ i < 167%21 ;) if ( do_search(i)) break;
for (int k = -1; ++ k < history_idx ;)
  cout << "FBRLUD"[history_mov[k ]] << history_rpt[k]+1 << " ";
return 0;
}
已有 1 人评分经验 收起 理由
ggglgq + 10 珍贵资料,收藏!

总评分: 经验 + 10   查看全部评分

Enjoy cubing
Enjoy coding.
我喜欢的公式 U D F2 B2 U' D'

Rank: 4

积分
1926
帖子
1058
精华
6
UID
17579
性别
保密

魔方理论探索者 论坛建设奖 六年元老

发表于 2010-1-20 19:32:06 |显示全部楼层
程序来源:http://tomas.rokicki.com/cubecontest/winners.html
这里面第一名的程序就是了。
原来的程序是用宏写的,比较难读点。我已经用g++展开了,还是如天书一般。

编译用g++

调用时按照如下顺序
UF UR UB UL DF DR DB DL FR FL BR BL UFR URB UBL ULF DRF DFL DLB DBR
输入

[ 本帖最后由 aubell 于 2010-1-20 19:38 编辑 ]
Enjoy cubing
Enjoy coding.
我喜欢的公式 U D F2 B2 U' D'

使用道具 举报

Rank: 4

积分
1031
帖子
757
精华
0
UID
87948
性别
保密
发表于 2010-1-20 19:37:26 |显示全部楼层
没看懂??????????????
敬请期待...........

使用道具 举报

Rank: 1

积分
23
帖子
19
精华
0
UID
1246974
性别
保密
发表于 2010-1-20 19:38:21 |显示全部楼层
看不懂啊.................

使用道具 举报

铜魔

007

Rank: 8Rank: 8

积分
13803
帖子
13083
精华
2
UID
101677
性别

四年元老 八年元老 十年元老

发表于 2010-1-20 19:39:42 |显示全部楼层
怎么不给个下载?代码看着晕…
魔方收藏群 123380874

使用道具 举报

Rank: 4

积分
2315
帖子
2249
精华
0
UID
110000
性别
发表于 2010-1-20 19:42:31 |显示全部楼层
不知道是什么啊……
本人还没学习过计算机编程啦

使用道具 举报

Rank: 3Rank: 3

积分
886
帖子
789
精华
1
UID
105914
性别
居住地
深圳市
兴趣爱好
速度
发表于 2010-1-20 19:48:21 |显示全部楼层
看不懂。。
变成  汗语 了
最近在练习怎么POP的更厉害。

使用道具 举报

铜魔

007

Rank: 8Rank: 8

积分
13803
帖子
13083
精华
2
UID
101677
性别

四年元老 八年元老 十年元老

发表于 2010-1-20 19:57:13 |显示全部楼层
才看到2楼~
不过,干什么用的?手机党,不明白
魔方收藏群 123380874

使用道具 举报

Rank: 4

积分
1926
帖子
1058
精华
6
UID
17579
性别
保密

魔方理论探索者 论坛建设奖 六年元老

发表于 2010-1-20 20:04:30 |显示全部楼层
测试了许多状态,一般都在30步左右。
最难得的是用如此简短的程序实现了状态集转换法!
而且运行速度相当快,并且不使用外部的表。
Enjoy cubing
Enjoy coding.
我喜欢的公式 U D F2 B2 U' D'

使用道具 举报

Rank: 4

积分
1926
帖子
1058
精华
6
UID
17579
性别
保密

魔方理论探索者 论坛建设奖 六年元老

发表于 2010-1-20 20:10:57 |显示全部楼层

回复 8# 的帖子

当然是解魔方用的。
例如打乱用 F- U+ F- D- L- D- F- U- L2 D- ,按照上面的顺序,就输入
progname RU LF UB DR DL BL UL FU BD RF BR FD LDF LBD FUL RFD UFR RDB UBL RBU

Aubell所作下一版的RCube将集成这个程序,作为内置的状态集转换法。实在钦佩原作者。
Enjoy cubing
Enjoy coding.
我喜欢的公式 U D F2 B2 U' D'

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

Archiver|手机版|魔方吧·中文魔方俱乐部

GMT+8, 2024-3-28 22:48

Powered by Discuz! X2

© 2001-2011 Comsenz Inc.

回顶部