loj.ac:#10024. 「一本通 1.3 练习 3」质数方阵
题目描述
质数方阵是一个\(5×5\)的方阵,每行、每列、两条对角线上的数字可以看作是五位的素数。方格中的行按照从左到右的顺序组成一个素数,而列按照从上到下的顺序。两条对角线也是按照从左到右的顺序来组成。这些素数每一位上的数之和必须相等。 左上角的数字是预先定好的。 一个素数可能在方阵中重复多次。不计含有前导 \(0\) 的五位素数,如\(00003\) 不是五位素数。
给出每一位上的数之和,以及左上角的数字,请输出方阵所有可能的填数方案。
如果不只有一个解,将它们全部输出(按照这 \(25\) 个数字组成的 \(25\) 位数的大小排序)。
输入格式
一行,包括两个被空格分开的整数:每一位上的数之和,以及左上角的数字。
输出格式
对于每一个找到的方案输出 \(5\) 行,每行 \(5\) 个字符,每行可以转化为一个 \(5\) 位的质数。在两组方案中间输出一个空行。如果没有解就单独输出一行 \(NONE\)。
样例
样例输入
11 1
样例输出
11351
14033
30323
53201
13313
11351
33203
30323
14033
33311
13313
13043
32303
50231
13331
不怎么华丽的分割线
这道题是真的水难,暴力中还带着技巧。
总的来说,我的填写方法有点奇特
填写顺序騛常奇特
这样填是我随便想的,结果过了 (除先填左上到右下)
代码很丑,照着填写顺序思路,代码写的很通俗,只要你不被判断绕晕
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
struct node
{
int map[10][10];
}ans[110000];
bool cmp(node a,node b)
{
for(int i=0;i<=4;i++)for(int j=0;j<=4;j++)if (a.map[i][j]!=b.map[i][j])return a.map[i][j]<b.map[i][j];
}
int a[10][10],n,k;
bool b[110000];
void bt()
{
k++;for(int i=0;i<=4;i++)for(int j=0;j<=4;j++)ans[k].map[i][j]=a[i][j];
}
void dfs6()
{
for(int i=1;i<=9;i++)
if(a[0][0]+a[0][4]+i<=n&&a[1][1]+a[2][1]+a[3][1]+i<n&&n-i-a[1][1]-a[2][1]-a[3][1]<=9&&!b[i*10000+a[1][1]*1000+a[2][1]*100+a[3][1]*10+n-i-a[1][1]-a[2][1]-a[3][1]])
{
a[0][1]=i,a[4][1]=n-i-a[1][1]-a[2][1]-a[3][1];
for(int j=1;j<=9;j++)
if (a[0][0]+a[0][1]+j+a[0][4]<n&&n-j-a[1][2]-a[2][2]-a[3][2]<=9&&j+a[1][2]+a[2][2]+a[3][2]<n&&n-a[0][0]-a[0][1]-j-a[0][4]<=9&&!b[j*10000+a[1][2]*1000+a[2][2]*100+a[3][2]*10+n-j-a[1][2]-a[2][2]-a[3][2]]&&!b[a[0][0]*10000+a[0][1]*1000+j*100+(n-a[0][0]-a[0][1]-j-a[0][4])*10+a[0][4]])
{
int n42=n-j-a[1][2]-a[2][2]-a[3][2],n03=n-a[0][0]-a[0][1]-j-a[0][4];
if(a[4][0]+a[4][1]+n42+a[4][4]<=n&&n-a[4][0]-a[4][1]-n42-a[4][4]<=9&&n03+a[1][3]+a[2][3]+a[3][3]<=n&&n-n03-a[1][3]-a[2][3]-a[3][3]<=9&&!b[a[4][0]*10000+a[4][1]*1000+n42*100+(n-a[4][0]-a[4][1]-n42-a[4][4])*10+a[4][4]]&&!b[n03*10000+a[1][3]*1000+a[2][3]*100+a[3][3]*10+n-n03-a[1][3]-a[2][3]-a[3][3]])
{
a[0][2]=j,a[4][2]=n42,a[0][3]=n03,a[4][3]=n-a[4][0]-a[4][1]-n42-a[4][4];
bt();
a[0][2]=a[4][2]=a[0][3]=a[4][3]=0;
}
}
a[0][1]=a[4][1]=0;
}
}
void dfs5()
{
for(int i=1;i<=9;i++)
if(a[3][3]+i<=n&&a[0][0]+a[1][0]+a[2][0]+i<n&&n-a[0][0]-a[1][0]-a[2][0]-i<=9&&!b[a[0][0]*10000+a[1][0]*1000+a[2][0]*100+i*10+n-a[0][0]-a[1][0]-a[2][0]-i])
{
a[3][0]=i,a[4][0]=n-a[0][0]-a[1][0]-a[2][0]-i;
for(int j=0;j<=9;j++)
if(a[1][1]+a[2][1]+j<=n&&a[3][0]+a[3][3]+j<=n&&a[4][0]+a[2][2]+a[1][3]+j<n&&n-a[4][0]-j-a[2][2]-a[1][3]<=9&&!b[a[4][0]*10000+j*1000+a[2][2]*100+a[1][3]*10+n-a[4][0]-j-a[2][2]-a[1][3]])
{
int n04=n-a[4][0]-j-a[2][2]-a[1][3];
if(n04+a[1][4]+a[2][4]+a[4][4]<=n&&n-n04-a[1][4]-a[2][4]-a[4][4]<=9&&!b[n04*10000+a[1][4]*1000+a[2][4]*100+(n-n04-a[1][4]-a[2][4]-a[4][4])*10+a[4][4]])
{
int n34=n-n04-a[1][4]-a[2][4]-a[4][4];
if(a[3][0]+j+a[3][3]+n34<=n&&n-a[3][0]-j-a[3][3]-n34<=9&&!b[a[3][0]*10000+j*1000+(n-a[3][0]-j-a[3][3]-n34)*100+a[3][3]*10+n34])
{
a[3][1]=j,a[0][4]=n04,a[3][4]=n34,a[3][2]=n-a[3][0]-j-a[3][3]-n34;
dfs6();
a[3][1]=a[0][4]=a[3][4]=a[3][2]=0;
}
}
}
a[3][0]=a[4][0]=0;
}
}
void dfs4()
{
for(int i=0;i<=9;i++)
if(a[1][1]+i<=n&&a[2][0]+a[2][2]+i<=n)
{
a[2][1]=i;
for(int j=0;j<=9;j++)
if(a[1][3]+j<=n&&a[2][0]+a[2][1]+a[2][2]+j<n&&n-a[2][0]-a[2][1]-a[2][2]-j<=9&&!b[a[2][0]*10000+a[2][1]*1000+a[2][2]*100+j*10+n-a[2][0]-a[2][1]-a[2][2]-j])
{
a[2][3]=j,a[2][4]=n-a[2][0]-a[2][1]-a[2][2]-j;
dfs5();
a[2][3]=a[2][4]=0;
}
a[2][1]=0;
}
}
void dfs3()
{
for(int i=0;i<=9;i++)
if(a[1][0]+a[1][1]+a[1][2]+i<n&&n-a[1][0]-a[1][1]-a[1][2]-i<=9&&a[3][3]+i<=n&&!b[a[1][0]*10000+a[1][1]*1000+a[1][2]*100+i*10+n-a[1][0]-a[1][1]-a[1][2]-i])
{
a[1][3]=i,a[1][4]=n-a[1][0]-a[1][1]-a[1][2]-i;
for(int j=1;j<=9;j++)
if(a[0][0]+a[1][0]+j<=n&&a[2][2]+j<=n)
{
a[2][0]=j;
dfs4();
a[2][0]=0;
}
a[1][3]=a[1][4]=0;
}
}
void dfs2()
{
for(int i=1;i<=9;i++)
if(a[0][0]+i<=n&&a[1][1]+i<=n)
{
a[1][0]=i;
for(int j=0;j<=9;j++)
if(a[1][0]+a[1][1]+j<=n&&a[2][2]+j<=n)
{
a[1][2]=j;
dfs3();
a[1][2]=0;
}
a[1][0]=0;
}
}
void dfs1()
{
for(int i=0;i<=9;i++)
if(a[0][0]+i<=n)
{
a[1][1]=i;
for(int j=0;j<=9;j++)
if(a[0][0]+a[1][1]+j<=n)
{
a[2][2]=j;
for(int g=0;g<=9;g++)
if(a[0][0]+a[1][1]+a[2][2]+g<=n&&n-a[0][0]-a[1][1]-a[2][2]-g<=9&&!b[a[0][0]*10000+a[1][1]*1000+a[2][2]*100+g*10+n-a[0][0]-a[1][1]-a[2][2]-g])
{
a[3][3]=g,a[4][4]=n-a[0][0]-a[1][1]-a[2][2]-g;
dfs2();
a[3][3]=a[4][4]=0;
}
a[2][2]=0;
}
a[1][1]=0;
}
}
int prime[110000],pr=0;
int main()
{
scanf("%d%d",&n,&a[0][0]);//
memset(prime,0,sizeof(prime));
memset(b,false,sizeof(b));
for(int i=2;i<=100000;i++)
{
if(b[i]==false)
{
prime[++pr]=i;
}
for(int j=1;(j<=pr)&& (i*prime[j]<=100000);j++)
{
b[i*prime[j]]=true;
if(i%prime[j]==0) break;
}
}
dfs1();
if(!k)printf("NONE\n");
else
{
sort(ans+1,ans+k+1,cmp);
for(int l=1;l<=k;l++)
{
for(int i=0;i<=4;i++)
{
for(int j=0;j<=4;j++)printf("%d",ans[l].map[i][j]);
printf("\n");
}
printf("\n");
}
}
return 0;
}
愿各位能活着看到这
(阿门)
loj.ac:#10024. 「一本通 1.3 练习 3」质数方阵的更多相关文章
- 「LOJ#10051」「一本通 2.3 例 3」Nikitosh 和异或(Trie
题目描述 原题来自:CODECHEF September Challenge 2015 REBXOR 1≤r1<l2≤r2≤N,x⨁yx\bigoplus yx⨁y 表示 ...
- LOJ#10117. 「一本通 4.1 练习 2」简单题
LOJ#10117. 「一本通 4.1 练习 2」简单题 题目描述 题目来源:$CQOI 2006$ 有一个$n$个元素的数组,每个元素初始均为$0$.有$m$条指令,要么让其中一段连续序列数字反转— ...
- LOJ#10064. 「一本通 3.1 例 1」黑暗城堡
LOJ#10064. 「一本通 3.1 例 1」黑暗城堡 题目描述 你知道黑暗城堡有$N$个房间,$M$条可以制造的双向通道,以及每条通道的长度. 城堡是树形的并且满足下面的条件: 设$D_i$为如果 ...
- 「LOJ#10056」「一本通 2.3 练习 5」The XOR-longest Path (Trie
#10056. 「一本通 2.3 练习 5」The XOR-longest Path 题目描述 原题来自:POJ 3764 给定一棵 nnn 个点的带权树,求树上最长的异或和路径. 输入格式 第一行一 ...
- LOJ #10131 「一本通 4.4 例 2」暗的连锁
LOJ #10131 「一本通 4.4 例 2」暗的连锁 给一棵 \(n\) 个点的树加上 \(m\) 条非树边 , 现在需要断开一条树边和一条非树边使得图不连通 , 求方案数 . $n \le 10 ...
- 「LOJ#10042」「一本通 2.1 练习 8」收集雪花 (map
题目描述 不同的雪花往往有不同的形状.在北方的同学想将雪花收集起来,作为礼物送给在南方的同学们.一共有 n 个时刻,给出每个时刻下落雪花的形状,用不同的整数表示不同的形状.在收集的过程中,同学们不希望 ...
- 「LOJ#10043」「一本通 2.2 例 1」剪花布条 (KMP
题目描述 原题来自:HDU 2087 一块花布条,里面有些图案,另有一块直接可用的小饰条,里面也有一些图案.对于给定的花布条和小饰条,计算一下能从花布条中尽可能剪出几块小饰条来呢? 输入格式 输入数据 ...
- 「LOJ#10015」「一本通 1.2 练习 2」扩散(并查集
题目描述 一个点每过一个单位时间就会向 444 个方向扩散一个距离,如图所示:两个点 a .b 连通,记作 e(a,b),当且仅当 a .b的扩散区域有公共部分.连通块的定义是块内的任意两个点 u.v ...
- #10042. 「一本通 2.1 练习 8」收集雪花 || 离散化 || 双指针法 || C++ || LOJ
题目:#10042. 「一本通 2.1 练习 8」收集雪花 看到网上没有这道题的题解,所以写一下. 要标记数字是否存在,看到x<=1e9,所以考虑用离散化,然后开一个last数组,last[i] ...
随机推荐
- Centos7.2(linux)minimal install之后需要的操作
minimal install之后,很多命令都不存在,例如ifconfig, wget等等 首先,需要先配置网络,保证机器可以连上互联纲 ip addr可以查看网卡的基本信息 一般默认就只有两个,一个 ...
- Windows环境下最新OpenCV和Contribute代码的联合编译【20180926更新红字】
解决这个问题,目的在于获得并使用最新的完全版本的代码,主要方法是对CMake能够熟练使用,并且对编译等基础支持有所了解. 因为这篇博客经过多次修改,所以里面的内容和配图可能有不是完全比对的地方,但是只 ...
- 关于Excel报错
今天,用POI做一个Excel的导出时,总是报空指针错误.后来经过我仔细的查找翻阅日志才发现问题的真正所在: 报错如下: HTTP Status 500 - Request processing fa ...
- Price Channel Breakout 交易系统简价及源码
Price Channel Breakout 交易系统简价及源码 既然这个版有交易系统的模型报告,小弟先在这裡野人献曝一下,把目前正在用的系统拿来请大家批评指教一下. ================ ...
- 当前 .NET SDK 不支持将 .NET Core 2.2 设置为目标。请将 .NET Core 2.1 或更低版本设置
场景重现 Windows 10 操作系统 (64-bit) 已经安装 .NET Core SDKs 如下: C:\Users\taadis>dotnet --list-sdks 2.1.202 ...
- [math]本博客已经支持书写数学公式
本博客已经支持mathjax格式公式 使用方法 使用方法单美元符号加单行公式. 使用方法双美元符号加多行公式. 展示 单行公式:\(x^2+2x+1=0\) 多行公式:\[x=\frac{{-b}\p ...
- TCGA一些数据库
最出名,http://www.cbioportal.org/ 特色:最基本的简单分析基因突变.共表达/共突变的基因,下载数据也可以,最常看的应该还是oncoPrint那个. 详细用法:TCGA数据库的 ...
- 分布式拒绝服务攻击(DDoS:Distributed Denial of Service)
DDoS攻击通过大量合法的请求占用大量网络资源,以达到瘫痪网络的目的. 指借助于客户/服务器技术,将多个计算机联合起来作为攻击平台,对一个或多个目标发动DDoS攻击,从而成倍地提高拒绝服务攻击的威力. ...
- 常见adb指令
1. adb –-help 查看帮助文档 2. adb start-server 当adb没有启动或被手动杀掉时,可以使用该命令启动服务 3. adb kill-server 杀死adb服务 4. a ...
- js中关于数组处理的一些小技巧
1 reduce方法同时实现map和filter 假设现在有一个数组,然后遍历它的每一项(map的功能)然后筛选出其中的一部分(filter的功能).如果使用map和filter的话,我们需要遍历这个 ...