难以言喻的满足感

#include
#include




int main(){
int N,i,j,k;
long double sum=0,x=0,y=0;


printf("N: ");
scanf("%d", &N);

for(i=1;i<=N;i++){
y=0;
for(j=1;j<=(i-1);j++){
x=0;
for(k=1;k<=(i-1);k++){
x=x-24*pow(-1,i+j+k)/sqrt(i*i+j*j+k*k);
}
y=y+x-24*pow(-1,i+j)/sqrt(i*i+j*j);
}


for(k=1;k<=(i-1);k++)
y=y-24*pow(-1,i+i+k)/sqrt(i*i+i*i+k*k);



sum=sum+y - 6*pow(-1,i)/i - 8*pow(-1,i)/(sqrt(3)*i) - 12/(sqrt(2)*i);
printf("N=%d\tM=%.20lf\n",i,sum);
}


printf("\nMadelung constant for NaCl at N=%d is\nM=%.20lf\n",N,sum);






printf("\n");






return 0;
}

都是自己爽的东西，

罢罢，

成功的满足感就是这样吧。