题目来源: mostleg
基准时间限制:1 秒 空间限制:131072 KB 分值: 40 难度:4级算法题
 收藏
 关注
如果一个由正整数组成的n*n的方阵,满足以下条件:

1,每个数字各不相同
2,每行以及每列的和,都是互不相同的完全平方数

我们称这种方阵为超级完全平方数方阵。
输入n,输出一个n*n的超级完全平方数方阵。如果存在多个方阵满足条件,输出将所有元素按行、列顺序排列后字典序最小的一个答案。例如
n=3时,下面两个方阵都符合条件

 1   2   6
 3   4   9
21 30 49

21 30 49
 3   4   9
 1   2   6

按行、列顺序排列后,第一个方阵表示为[1, 2, 6, 3, 4, 9, 21, 30, 49],第二个方阵表示为[21, 30, 49, 3, 4, 9, 1, 2, 6]。第一个方阵字典序更小一些。

如果不存在这样的方阵,输出No Solution。

Input
仅一行,为一个正整数n。(1 <= n <= 64)
Output
输出n行,每行为n个整数,之间用空格隔开,表示所求的n*n方阵。或者,输出No Solution。
Input示例
3
Output示例
1 2 6
3 4 9
21 30 49

真真正正地被虐了一下午。。。其实本质上就是一个dfs,但是做起来是真的麻烦啊,各种错误百出的。

官方题解:

首先,n=1时无解。


接下来处理n>=2的情况。由于题目要求字典序最小的方阵,使用贪心算法的思想,不难发现,每一行每一列其实只需要依靠最后一个数字(最右边和最下边的数字)就足够使得该行该列的和达到一个没有使用过的完全平方数。因此,按照题目中对方阵序列化的次序,对无关紧要的位置都尽力使用最小的数字;每当到达一行的最后一个位置,或者最后一行的时候,再去寻找符合题目要求的最后一个数字。这样做直到右下角的最后一个位置。

此时,最后一行和最后一列都需要满足和为完全平方数的条件。搜索最小的符合条件的数字。如果找不到解,就加大倒数第二个位置的数字(因为这样做对字典序的影响最小),再重新搜索最后一个位置。

怎样快速发现最后一个位置找不到解呢?不难发现,最后一列的和必定小于最后一行的和,设它们的差为d。我们可以枚举较小的一个完全平方数x,如果发现x的下一个完全平方数与x的差已经大于d,则在最后一个位置无解。

代码:

#include <iostream>
#include <algorithm>
#include <cmath>
#include <vector>
#include <string>
#include <cstring>
#include <map>
#pragma warning(disable:4996)
using namespace std; int n,wang=0;
int square_flag[10000];
int val_flag[64005];
long long val[70][70]; int sear(int su)
{
long long i;
for (i = 2; ; i++)
{
if (i*i >= su && ((i<=9999&&square_flag[i]==0)||(i>9999)))
return i;
}
} void dfs(int x,int y,long long value)
{
if (wang==1)
{
return;
}
if (x == n&&y == n)
{
long long i, j, h, k, sum2 = 0, sum3 = 0;
for (i = 1; i <= n - 1; i++)
{
sum2 += val[i][y];
}
for (i = 1; i <= n - 1; i++)
{
sum3 += val[x][i];
}
for (i = 2;; i++)
{
if ((i <= 9999 && square_flag[i] == 1))continue;
double g = sqrt((double)(i*i - sum2 + sum3)); if (i*i - sum2 > 0 && g == (long long)g && ((g<=9999)&&(square_flag[(long long)g] == 0)||g>9999))
{
val[x][y] = i*i - sum2;
for (h = 1; h <= n; h++)
{
for (k = 1; k <= n; k++)
{
cout << val[h][k]<< " ";
}
cout << endl;
}
wang = 1;
return;
}
long long wa = sum3 - sum2;
if ((i + 1)*(i + 1) - (i*i) > wa)
{
long long op, sum_op = 0;
for (op = 1; op <= n; op++)
{
sum_op += val[op][y - 1];
}
square_flag[(long long)sqrt((double)sum_op)] = 0;
dfs(x,y-1,value+1);
return;
}
}
}
else if (x == n)
{
long long i, sum2 = 0;
for (i = 1; i <= n - 1; i++)
{
sum2 += val[i][y];
}
i = sear(sum2 + value);
while (val_flag[i*i - sum2] == 1||square_flag[i]==1)
{
i++;
}
val[x][y] = i*i - sum2;
val_flag[i*i - sum2] = 1;
square_flag[i] = 1; dfs(x, y+1, value);
}
else if (y == n)
{
long long i,sum2=0;
for (i = 1; i <= n - 1; i++)
{
sum2 += val[x][i];
}
i = sear(sum2 + value); while (val_flag[i*i - sum2] == 1)
{
i++;
}
val[x][y] = i*i - sum2;
val_flag[i*i - sum2] = 1;
square_flag[i] = 1; dfs(x + 1, 1, value);
}
else
{
val[x][y] = value;
val_flag[value] = 1;
if (val_flag[value+1] == 0)
{
dfs(x, y + 1, value+1);
}
else
{
while (val_flag[value+1] == 1)
{
value++;
}
dfs(x, y + 1, value+1);
}
}
} int main()
{
scanf("%d", &n);
if (n == 1)
{
cout << "No Solution" << endl;
}
else
{
memset(square_flag,0,sizeof(square_flag));
memset(val_flag, 0, sizeof(val_flag)); dfs(1, 1, 1);
}
return 0;
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

51nod 1438:方阵与完全平方数的更多相关文章

  1. 51Nod 欢乐手速场1 C 开心的小Q[莫比乌斯函数]

    开心的小Q tangjz (命题人) quailty (测试)   基准时间限制:1 秒 空间限制:131072 KB 分值: 80 如果一个数字存在一个约数是完全平方数,那么小Q就认为这个数是有趣的 ...

  2. 51nod1787最大子方阵

    51nod1787最大子方阵 我在51nod上面切的第一道题 我在51nod上面切的第一道8级题 我在51nod上面切的第一道8级题的一血 题目大意 有一个n*m的矩阵,矩阵中的每一个元素是'X'或者 ...

  3. BZOJ 2440: [中山市选2011]完全平方数 [容斥原理 莫比乌斯函数]

    2440: [中山市选2011]完全平方数 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 3028  Solved: 1460[Submit][Sta ...

  4. 【51Nod 1244】莫比乌斯函数之和

    http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1244 模板题... 杜教筛和基于质因子分解的筛法都写了一下模板. 杜教筛 ...

  5. [LeetCode] Valid Perfect Square 检验完全平方数

    Given a positive integer num, write a function which returns True if num is a perfect square else Fa ...

  6. [LeetCode] Perfect Squares 完全平方数

    Given a positive integer n, find the least number of perfect square numbers (for example, 1, 4, 9, 1 ...

  7. 51Nod 1268 和为K的组合

    51Nod  1268  和为K的组合 1268 和为K的组合 基准时间限制:1 秒 空间限制:131072 KB 分值: 20 难度:3级算法题 给出N个正整数组成的数组A,求能否从中选出若干个,使 ...

  8. 51Nod 1428 活动安排问题

    51Nod   1428  活动安排问题 Link: http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1428 1428 活 ...

  9. 51Nod 1278 相离的圆

    51Nod 1278 相离的圆 Link: http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1278 1278 相离的圆 基 ...

随机推荐

  1. 新闻网大数据实时分析可视化系统项目——19、Spark Streaming实时数据分析

    1.Spark Streaming功能介绍 1)定义 Spark Streaming is an extension of the core Spark API that enables scalab ...

  2. uboot源码分析1-启动第一阶段

    1.不简单的头文件包含 #include <config.h>:这个文件的内容其实是包含了一个头文件:#include <configs/x210_sd.h>". # ...

  3. 006.Delphi插件之QPlugins,多服务演示

    演示效果如下 演示工程,全部就一个文件,代码如下 unit Frm_Main; interface uses Winapi.Windows, Winapi.Messages, System.SysUt ...

  4. CentOS7 防火墙设置

    CentOS7 防火墙命令 最近在公司服务器上安装了oracle12c数据库,在用数据库客户端连接的时候,连接不了.最后查找资料的原因是因为oracle的服务端口未开放. 首先还是还是输入以往的开启某 ...

  5. AS3.0判断数组中最大值

    function getMax(Arr) { if (typeof Arr !="object") {     return null;  }  for (var i=0,max= ...

  6. C语言备忘录——qsort

    写了这么久的排序感觉还是用现成的最舒服.其实C语言其实自己带了一个快速排序在stdlib 库里,但是其函数调用的接口过于复杂,所以让人望而却步.为了加深自己的记忆,所以写下这篇博客 先来看一下函数原型 ...

  7. js实现二叉查找树

    二叉树的特点:   像一颗树一样,从顶端往下延伸,最顶端的为根节点,每个节点下面子节点的数不超过两个,没有任何子节点的节点被称为叶子节点, 除了根节点和叶子节点的被称为中间节点. 二叉查找树: 每个节 ...

  8. [LeetCode] 929. Unique Email Addresses 独特的邮件地址

    Every email consists of a local name and a domain name, separated by the @ sign. For example, in ali ...

  9. 五、JavaScript之点击按钮调用相关JavaScript函数

    一.代码如下 二.执行之后 三.点击按钮之后,内容被改变

  10. 六十、SAP中的加减乘除等算数运算符

    一.代码和显示如下,很奇怪的事情是,负数的负号居然在数字后面