一秒记住【笔趣阁】biquge365.net，更新快，无弹窗！
第90章秒掉莫斯科
莫斯科的清晨来得早。
六点半，天已经大亮。
陆沉睁开眼，脑海中的倒计时刚好归零四个小时的睡眠，足够恢复精力了。
林枫还在睡，呼吸均匀。
昨晚他在终端前待到十一点，反覆熟悉BESM—6的指令集，最后是被值班老师赶回来的。
陆沉没有叫醒他。他安静地洗漱，然后坐在床边，闭上眼睛。
【写到这里我希望读者记一下我们域名台湾小说网超便捷，?????.???随时享】
意识沉入那片他称之为「构建空间」的领域。
这是他的超算大脑真正施展威力的地方。
在他的感知中，BESM—6的硬体结构像一张三维蓝图般展开。
寄存器组丶内存分页丶指令预取队列丶微码存储器—每一个模块都清晰可见，每一个信号传递的时序都精确到时钟周期。
他开始构建优化版数学库。
标准库的浮点乘法子程序需要147个时钟周期。
他的目标是：89个。
方法很简单—至少原理上很简单。标准库为了兼容所有数据类型和边界情况，加入了大量判断分支。而他只需要针对团体赛可能的计算场景，裁剪掉不需要的通用性，用最直接的寄存器操作完成运算。
但这需要精确到每一个微码步骤的调度。
如果调度错一个周期，乘法的部分积就会错位。如果寄存器冲突没处理好，流水线就会空转。
他需要同时考虑：
指令发射顺序寄存器依赖关系内存访问延迟微码跳转的预判惩罚这些变量互相耦合，构成一个极其复杂的优化问题。
但对于他而言—
这只是一道有趣的题。
七点十五分，林枫醒了。
他看到陆沉坐在床边，双眼微闭，以为他在打坐。
「修仙呢？」林枫揉着眼睛坐起来。
陆沉睁开眼。
「写完了。」
「啊？」
「数学库。」陆沉说，「乘除法丶浮点运算丶矩阵基本操作，大概四十几个函数，比赛的时候直接用就行。」
林枫沉默了五秒钟。
「————你一早上写的？」
「嗯。」
「在脑子里？」
陆沉想了想，觉得解释起来太复杂，于是只说：「差不多。
林枫张了张嘴，最终什么也没说，默默起床洗漱去了。
八点整，中国队集合完毕，前往计算中心。
团体赛的算法设计模块在主赛场进行。
十二支代表队各占据一片区域，每队一台BESM—6终端，配发题目和参考手册。
陆沉注意到，苏联队坐最前排，东德队和波兰队分列两侧。
中国队的位置在第三排左侧，旁边是匈牙利队。
落座时，一个戴眼镜的匈牙利男生朝他们点了点头，用口音浓重的英语说：「祝好运」
。
陆沉点头回礼。
八点三十分，主考官走到讲台前。
那是位五十多岁的苏联教授，灰白的头发梳得一丝不苟，他用俄语宣读比赛规则，旁边有翻译用英语复述。
「本次算法设计模块的题目，由苏联科学院计算中心提供。」
「限时四小时。」
「评分标准：正确性占50%，运行效率占40%，代码结构占10%。」
「可以使用任何程式语言，但评测将在BESM—6实机上运行。」
「现在，分发题目。」
题目是密封在牛皮纸袋里的。陆沉拆开，抽出里面的纸张。
只有一页纸。
上面是一段简短的描述，和一组数学公式。
陆沉的目光扫过题目。
然后，他的瞳孔微微收缩了一下。
不是因为难。
是因为他见过这道题。
准确地说，是前世见过。
这道题的核心，是求解大规模稀疏线性方程组的预条件共轭梯度法—一种在计算数学中极其重要的叠代算法。
预条件技术还处于萌芽阶段。
共轭梯度法虽然在五十年代就被提出，但因为数值稳定性和收敛速度的问题，一直没能在实际工程中广泛使用。
直到七十年代末，预条件技术的突破才让这个方法真正实用化。
而这道题，恰恰要求选手设计一个预条件子，并实现完整的求解流程。
陆沉放下题目。
他知道这道题的标准解法一不完整的Cholesky分解作为预条件子，结合共轭梯度叠代。这是他前世做数值计算时烂熟于心的东西。
但问题在于——
1988年，不完整的Cholesky分解的论文还没有发表。
原始论文是1990年才出现的。
如果他直接拿出这个解法，会引起什么反应？
不过这个念头只停留了一瞬。
他是来比赛的。
至于之后的事情，之后再说。
陆沉坐到终端前，开始敲击键盘。
中国队的分工是：陆沉负责核心算法设计，林枫负责输入输出和测试用例，另外两名队员负责文档和辅助验证。
但开赛十五分钟后，林枫就发现了一个问题。
「陆沉。」他压低声音，「你写的这个预条件子————是什么方法？」
「不完整的Cholesky分解。」
林枫翻了翻参考手册：「手册上没有。」
「嗯。
」
「那你怎么————」
「推导出来的。」陆沉没有停下手上的动作，「共轭梯度法的收敛速度取决于系数矩阵的条件数。如果能找到一个近似分解，使得预条件后的矩阵特徵值聚集，就能加速收敛。不完整分解的核心思想是，只保留原始矩阵稀疏结构内的元素，其他位置直接置零。」
林枫沉默了。
这些话他每一个字都听得懂，但连在一起就变成了一堵墙。
他决定不再问，专心做自己能做的事。
一个小时后，陆沉完成了核心代码。
他靠在椅背上，闭上眼睛，让代码在脑中的BESM—6模拟器上跑了一遍。
正确。
效率标准解法的约60%运行时间。
够用了。
他重新睁开眼睛，准备继续完善文档。
就在这时，他注意到赛场前方出现了一阵轻微的骚动。
苏联队的区域里，几个队员正在低声争论什么。
主考官走过去，听了一会儿，然后点了点头，对全场说：「苏联队请求更换终端。原终端出现故障。」
工作人员开始检查设备，陆沉的目光落在苏联队队长的脸上——那是一个十七八岁的金发少年，眉头紧锁，神情有些焦躁。
不是因为设备故障。
设备故障是藉口。
真正的原因是：他们被题目卡住了。
预条件共轭梯度法这道题，如果没有不完整分解的思路，就只能用最简单的对角预条件雅可比预条件。
但雅可比预条件对这道题的病态矩阵效果很差，收敛极慢，甚至可能不收敛。
苏联队大概试了雅可比，发现跑不通，以为是机器出了问题。
陆沉收回目光，继续写文档。
这是比赛。
赛场上，没有义务提醒对手。
十一点四十分，距离比赛结束还有五十分钟。
陆沉完成了全部代码和文档。林枫的测试用例也准备好了。
「提交吗？」林枫问。
「等一下。」陆沉说，「我再优化一个地方。」
他重新打开代码文件，定位到预条件子的构造部分。
其实现在提交，成绩已经足够好了。
不完整的Cholesky分解比雅可比预条件快至少一个数量级，中国队在这个模块上的优势是碾压级的。
但陆沉忽然想试试另一个思路。
如果不完整的Cholesky分解是1980年的解法，那么—
他能不能把时间线再往前推一点？
比如，1990年代才成熟的多重网格预条件？
不。时间不够。四十分钟写不完一个完整的多重网格框架。
但有一个折中的思路。
他可以在不完整分解的基础上，加入一个简单的块对角预处理，针对这道题矩阵的特殊分块结构。
这个改进不大，但能进一步压缩15%的运行时间。
而且，可以在现有代码的基础上用十分钟改完。
他动手了。
十一点五十五分。
陆沉按下保存键。
「提交。」
林枫长长地呼出一口气，像是终于卸下了什么重担。
评测过程是公开的。
各队的程序被依次送入BESM—6主机，在统一的测试集上运行，运行时间和结果实时显示在赛场前方的大屏幕上。
第一队，保加利亚。运行结果：部分正确。运行时间：47.3秒。
第二队，捷克斯洛伐克。运行结果：正确。运行时间：42.1秒。
第三队，东德。运行结果：正确。运行时间：38.7秒。
随着结果一个个公布，赛场里的气氛渐渐紧张起来。正确完成这道题的队伍并不多，大多数要么结果错误，要么运行超时。
第六队，匈牙利。运行结果：正确。运行时间：35.2秒。
匈牙利队的区域里发出一阵压低声音的欢呼。戴眼镜的男生朝陆沉这边看了一眼，比了个大拇指。
第七队，苏联。运行结果：正确。运行时间：
屏幕上跳动了一下。
31.4秒。
苏联队的区域里，队员们终于露出了笑容。金发队长靠在椅背上，神色明显放松下来0
31.4秒。目前全场最快。
第八队，波兰。运行结果：正确。运行时间：33.8秒。
第九队，罗马尼亚。运行结果：部分正确。运行时间：超时。
第十队，中国。
陆沉感觉到林枫的身体绷紧了。
屏幕上的数字开始跳动。
加载程序。
初始化数据。
开始计算。
计时器的数字飞速增长。
0.5秒。1.0秒。1.5秒。
林枫死死盯着屏幕。
陆沉倒是很平静。
他已经知道了结果。
2.1秒。
计时停止。
运行结果：正确。
.
——
运行时间：2.1秒。
赛场安静了一瞬。
然后「什么？」
声音是从苏联队的区域传出来的。
金发队长站了起来，盯着屏幕上的数字，脸上的表情从困惑变成难以置信。
2.1秒。
比第二名快了将近十五倍。
这不是「更快」。
这是不在同一个维度。
主考官也愣住了。他低头检查评测系统的日志，确认没有出错。
没有出错。
日志显示，中国队的程序确实在2.1秒内完成了全部计算任务。
而且，内存占用只有标准解法的三分之一。
主考官抬起头，看向中国队的方向。
「你们————」他用俄语说，然后切换成英语，「你们的预条件子，用的是什么方法？
」
翻译还没开口，陆沉已经用俄语回答了。
「不完整的Cholesky分解，结合块对角预处理。」
主考官的眼神变了。
「不完整的Cholesky分解？」他重复了一遍，「你是指————对系数矩阵做不完全的三角分解，舍弃填充元？」
「是的。」
「这个方法的收敛性」」
「有保证。」陆沉说，「对于M矩阵，不完整分解的存在性和稳定性都可以证明。具体来说，如果系数矩阵是M矩阵，则任意舍弃填充元的不完整Cholesky分解都存在，且预条件后的矩阵特徵值均为正实数。这道题的矩阵正好满足M矩阵条件。」
主考官沉默了好一会儿。
「你————证明过吗？」
「推导过。」陆沉说，「证明过程在我的文档附录里。」
主考官快步走向评测终端，调出中国队提交的文档。
屏幕上的内容一页页翻过。
算法描述。复杂度分析。数值实验。
然后，他看到了附录。
那里有一段简洁而严谨的证明，从M矩阵的定义出发，逐步推导出不完整Cholesky分解的存在性和稳定性条件。
推导过程乾净得像一场外科手术。
主考官看了很久。
然后他转过身，用一种很奇怪的语气说：「中国队。算法设计模块。满分。」
他顿了顿，又补充了一句：「附加分。10分。」
附加分。
这是本届比赛第一次出现附加分。
林枫猛地转过头看陆沉，眼睛瞪得像铜铃。
陆沉只是平静地收拾着桌上的纸笔。
但他的手指，微微收紧了。
不是因为紧张。
是因为他知道，这件事不会就这么结束。
一个十岁的中国孩子，在莫斯科的赛场上，拿出了一套比世界最前沿还要超前的数值算法。
这不是数学竞赛。
这是某种信号。
他抬头看了一眼赛场后方的观察席。
那里坐着几个没有佩戴参赛证件的人。
其中一个，是索科洛夫。
索科洛夫正低头在笔记本上写着什么，似乎感觉到陆沉的目光，抬起头来。
隔着整个赛场的距离，两人再次对视。
索科洛夫合上了笔记本。
他的表情，看不出任何情绪。
但他的眼神，像是在看一件他找了很久的东西。
