要用编程还原魔方，你可以选择层序法或CFOP法。以下是使用Python语言和pygame库实现层序法还原魔方的示例代码：
```python
import pygame
from pygame.math import Vector3
from pygame.locals import *
from OpenGL.GL import *
from OpenGL.GLUT import *
定义魔方的基本结构
class Cubie:
def __init__（self, pos, color）:
self.pos = pos
self.color = color
class Rubik:
def __init__（self）:
self.cubies = [
Cubie（Vector3（x, y, z）, （random（）, random（）, random（）））
for x in range（-1, 2）
for y in range（-1, 2）
for z in range（-1, 2）
]
def display（self）:
for cubie in self.cubies:
在这里添加代码来显示每个小立方体的颜色
pass
def rotate（self, axis, angle）:
在这里添加代码来实现魔方的旋转操作
pass
def solve（self）:
在这里添加代码来实现魔方的还原算法，例如层序法
pass
初始化pygame和OpenGL
pygame.init（）
display = （800, 600）
pygame.display.set_mode（display, DOUBLEBUF | OPENGL）
设置OpenGL的投影矩阵
glMatrixMode（GL_PROJECTION）
glLoadIdentity（）
glTranslatef（-0.5, 0.5, 0.0）
创建Rubik对象
rubik = Rubik（）
主循环
while True:
for event in pygame.event.get（）:
if event.type == pygame.QUIT:
pygame.quit（）
quit（）
在这里添加代码来更新魔方的状态和显示
pygame.display.flip（）
```
在这个示例中，我们首先定义了`Cubie`类来表示魔方的每个小立方体，以及`Rubik`类来组织整个魔方。`Rubik`类中包含了`display`方法来显示魔方的状态，`rotate`方法来实现魔方的旋转操作，以及`solve`方法来实现魔方的还原算法。
请注意，这个示例代码只是一个框架，你需要根据具体的算法和显示需求来填充`display`、`rotate`和`solve`方法中的代码。对于CFOP法，你需要实现交叉（Cross）、F2L（First Two Layers）、OLL（Orient Last Layer）和PLL（Permute the Last Layer）四个步骤，每个步骤都包含一系列旋转操作。
此外，如果你想要一个更复杂的3D可视化效果，你可以使用OpenGL库来渲染魔方的每个面，并实现更精细的旋转和移动操作。