挑战：以任何语言实现Delacorte数的计算。最短的代码获胜。

对于不同整数的一个给定的方阵1..n²（可能边长Ñ在3和27之间至少），其德拉科特数是产品的总和满足gcd（A，B）×distance²（A，B）为每个不同的对整数{a，b}。

以下示例显示了Delacorte数为160的3×3正方形。

3 2 9
4 1 8
5 6 7

在这个正方形中，我们有36个不同的对来计算，例如对4和6：gcd（4，6）× 距离 ²（4，6）= 4

另一个用于测试的示例正方形-的Delacorte数为5957：

10  8 11 14 12
21  4 19  7  9
 5 13 23  1 16
18  3 17  2 15
24 22 25  6 20

Delacorte Numbers摘自该编程竞赛 -有关详细信息，请参阅此处。竞赛于2015年1月结束。这真是太好了！

规则：

必要的换行符计为1个字符。您可以在打高尔夫球的解决方案中添加换行符，但是只有在必要时才使用该语言进行计数。

您可以选择处理输入和输出的方式，而不必计算必要的语言框架，例如标准包含或主函数标头。仅实际的代码计数（包括快捷方式/别名定义），例如此C＃示例：

namespace System
{
    using Collections.Generic;
    using I=Int32; //this complete line counts
    class Delacorte
    {
        static I l(I[]a){return a.Length;} //of course this complete line counts

        static void CalculateSquare(int[] a, out int r)
        {
            r=0;for(I i=l(a);i-->0;)r+=a[i]; //here only this line counts
        }

        static void Main()
        {
            int result;
            CalculateSquare(new int[] { 1, 2, 3, 4, 5, 6, 7, 8, 9 }, out result);
            Console.Write(result); //should output 140 for the example
            Console.ReadKey();
        }
    }
}

您也可以将正方形输入为二维数组，也可以从提示输入，也可以输入为字符串或某些标准收集类型。二维数组是不必自己计算正方形的边长的唯一方法。
不需要用于实际工作的子函数，您也可以将代码直接放在Main（）中。

甚至可以免费进行更多准备，例如：

using System;
unsafe class Delacorte
{
    static void CalculateSquare(int* a, out int r)
    {
        r=0;while(*a>0)r+=*a++; //only this line counts
    }

    static void Main()
    {
        var input = new int[] { 1, 2, 3, 4, 5, 6, 7, 8, 9, 0 }; //adding a terminator
        int result;
        fixed (int* a = &input[0]) //necessary in C#
            CalculateSquare(a, out result);
        Console.Write(result);
        Console.ReadKey();
    }
}

如果您不确定自己的冗长准备是否符合这些规则的精​​神，或者可以称为作弊，请询问:)

APL（38）

{.5×+/∊∘.{(∨/Z[⍺⍵])×+/⊃×⍨⍺-⍵}⍨⊂¨⍳⍴Z←⍵}

这是一个将矩阵作为正确参数的函数，如下所示：

      sq5←↑(10 8 11 14 12)(21 4 19 7 9)(5 13 23 1 16)(18 3 17 2 15)(24 22 25 6 20)
      sq5
10  8 11 14 12
21  4 19  7  9
 5 13 23  1 16
18  3 17  2 15
24 22 25  6 20
      {.5×+/∊∘.{(∨/Z[⍺⍵])×+/⊃×⍨⍺-⍵}⍨⊂¨⍳⍴Z←⍵}sq5
5957

说明：

• ⊂¨⍳⍴Z←⍵：将矩阵存储在中Z。在中列出每个可能的坐标对Z。
• ∘.{... }⍨：对于每对坐标，并结合每对坐标：
  • +/⊃×⍨⍺-⍵：计算distance^2：从第二对中减去第一对坐标，将它们自己相乘并求和
  • ∨/Z[⍺⍵]：获取Z两个坐标对的数字，并找到GCD
  • ×：彼此相乘
• +/∊：求和结果的元素
• .5×：乘以0.5（因为我们之前对每个非零对进行了两次计数）

Mathematica（83 82 79 69 67 66）

制备

a={{10,8,11,14,12},{21,4,19,7,9},{5,13,23,1,16},{18,3,17,2,15},{24,22,25,6,20}}

码

#/2&@@Tr[ArrayRules@a~Tuples~2/.{t_->u_,v_->w_}->u~GCD~w#.#&[t-v]]

如果我们使用Unicode字符计数：62：

Tr[ArrayRules@a~Tuples~2/.{t_u_,v_w_}u~GCD~w#.#&[t-v]]〚1〛/2

蟒蛇-128112 90 89 88

制备：

import pylab as pl
from fractions import gcd
from numpy.linalg import norm
from itertools import product

A = pl.array([
    [10,  8, 11, 14, 12],
    [21,  4, 19,  7,  9],
    [ 5, 13, 23,  1, 16],
    [18,  3, 17,  2, 15],
    [24, 22, 25,  6, 20]])

计算Delacorte数（计数的行）：

D=sum(gcd(A[i,j],A[m,n])*norm([m-i,n-j])**2for j,n,i,m in product(*[range(len(A))]*4))/2

输出：

print D

结果：

5957

腐霉病 43

这个答案几乎可以肯定会进一步发展。我特别不喜欢距离计算。

K^lJ.5VJFdUN~Z*i@JN@Jd+^-/dK/NK2^-%dK%NK2;Z

要进行设置，请将线性化的数组存储在变量J中。您可以这样编写：

J[3 2 9 4 1 8 5 6 7)

在线尝试。

输出浮点数。我认为这是合法的，如果我违反规则，请告诉我：)

说明：

                                             : Z=0 (implicit)
K^lJ.5                                       : K=sqrt(len(J))
      VJ                                     : for N in range(len(J))
        FdUN                                 : for d in range(N)
            ~Z*                              : Z+= the product of
               i@JN@Jd                       : GCD(J[N],J[d])
                      +^-/dK/NK2^-%dK%NK2    : (d/K-N/K)^2 + (d%K-N%K)^2 (distance)
                                         ;Z  : end all loops, and print Z

果酱，55

q~:Q__,mqi:L;m*{_~{_@\%}h;\[Qf#_Lf/\Lf%]{~-_*}/+*}%:+2/

采用以下格式将矩阵作为STDIN：

[10  8 11 14 12
 21  4 19  7  9
  5 13 23  1 16
 18  3 17  2 15
 24 22 25  6 20]

在这里在线尝试