P2615 [NOIP2015 提高组] 神奇的幻方

[NOIP2015 提高组] 神奇的幻方

题目背景

NOIp2015 提高组 Day1T1

题目描述

幻方是一种很神奇的 $N\times N$ 矩阵：它由数字 $1,2,3,\cdots \cdots ,N\times N$ 构成，且每行、每列及两条对角线上的数字之和都相同。

当 $N$ 为奇数时，我们可以通过下方法构建一个幻方：

首先将 $1$ 写在第一行的中间。

之后，按如下方式从小到大依次填写每个数 $K(K=2,3,\cdots ,N\times N)$ ：

1. 若 $(K-1)$ 在第一行但不在最后一列，则将 $K$ 填在最后一行， $(K-1)$ 所在列的右一列；
2. 若 $(K-1)$ 在最后一列但不在第一行，则将 $K$ 填在第一列， $(K-1)$ 所在行的上一行；
3. 若 $(K-1)$ 在第一行最后一列，则将 $K$ 填在 $(K-1)$ 的正下方；
4. 若 $(K-1)$ 既不在第一行，也不在最后一列，如果 $(K-1)$ 的右上方还未填数，则将 $K$ 填在 $(K-1)$ 的右上方，否则将 $K$ 填在 $(K-1)$ 的正下方。

现给定 $N$ ，请按上述方法构造 $N\times N$ 的幻方。

输入格式

一个正整数 $N$，即幻方的大小。

输出格式

共 $N$ 行，每行 $N$ 个整数，即按上述方法构造出的 $N\times N$ 的幻方，相邻两个整数之间用单空格隔开。

样例 #1

样例输入 #1

3

样例输出 #1

8 1 6
3 5 7
4 9 2

样例 #2

样例输入 #2

25

样例输出 #2

327 354 381 408 435 462 489 516 543 570 597 624 1 28 55 82 109 136 163 190 217 244 271 298 325
353 380 407 434 461 488 515 542 569 596 623 25 27 54 81 108 135 162 189 216 243 270 297 324 326
379 406 433 460 487 514 541 568 595 622 24 26 53 80 107 134 161 188 215 242 269 296 323 350 352
405 432 459 486 513 540 567 594 621 23 50 52 79 106 133 160 187 214 241 268 295 322 349 351 378
431 458 485 512 539 566 593 620 22 49 51 78 105 132 159 186 213 240 267 294 321 348 375 377 404
457 484 511 538 565 592 619 21 48 75 77 104 131 158 185 212 239 266 293 320 347 374 376 403 430
483 510 537 564 591 618 20 47 74 76 103 130 157 184 211 238 265 292 319 346 373 400 402 429 456
509 536 563 590 617 19 46 73 100 102 129 156 183 210 237 264 291 318 345 372 399 401 428 455 482
535 562 589 616 18 45 72 99 101 128 155 182 209 236 263 290 317 344 371 398 425 427 454 481 508
561 588 615 17 44 71 98 125 127 154 181 208 235 262 289 316 343 370 397 424 426 453 480 507 534
587 614 16 43 70 97 124 126 153 180 207 234 261 288 315 342 369 396 423 450 452 479 506 533 560
613 15 42 69 96 123 150 152 179 206 233 260 287 314 341 368 395 422 449 451 478 505 532 559 586
14 41 68 95 122 149 151 178 205 232 259 286 313 340 367 394 421 448 475 477 504 531 558 585 612
40 67 94 121 148 175 177 204 231 258 285 312 339 366 393 420 447 474 476 503 530 557 584 611 13
66 93 120 147 174 176 203 230 257 284 311 338 365 392 419 446 473 500 502 529 556 583 610 12 39
92 119 146 173 200 202 229 256 283 310 337 364 391 418 445 472 499 501 528 555 582 609 11 38 65
118 145 172 199 201 228 255 282 309 336 363 390 417 444 471 498 525 527 554 581 608 10 37 64 91
144 171 198 225 227 254 281 308 335 362 389 416 443 470 497 524 526 553 580 607 9 36 63 90 117
170 197 224 226 253 280 307 334 361 388 415 442 469 496 523 550 552 579 606 8 35 62 89 116 143
196 223 250 252 279 306 333 360 387 414 441 468 495 522 549 551 578 605 7 34 61 88 115 142 169
222 249 251 278 305 332 359 386 413 440 467 494 521 548 575 577 604 6 33 60 87 114 141 168 195
248 275 277 304 331 358 385 412 439 466 493 520 547 574 576 603 5 32 59 86 113 140 167 194 221
274 276 303 330 357 384 411 438 465 492 519 546 573 600 602 4 31 58 85 112 139 166 193 220 247
300 302 329 356 383 410 437 464 491 518 545 572 599 601 3 30 57 84 111 138 165 192 219 246 273
301 328 355 382 409 436 463 490 517 544 571 598 625 2 29 56 83 110 137 164 191 218 245 272 299

提示

对于 $100\%$ 的数据，对于全部数据， $1\le N\le 39$ 且 $N$ 为奇数。

---------------------------------------------------------分割线---------------------------------------------------------
作为一个蒟蒻,我第一时间想到的是模拟,以下分析全为个人理解:
引用Linux 的创始人 Linus Torvalds 的一句话:

CODE

//给个免费的赞吧,谢谢(●'◡'●)!
#include
using namespace std;
pair<int,int> vfind(int target,vector<vector<int>>arr) {
	pair<int,int>pos{1,1};
	for(vector<vector<int>>::iterator it1=arr.begin()+1;it1!=arr.end(); ++it1)//n行 
	{
		for(vector<int>::iterator it2=(*it1).begin()+1; it2!=(*it1).end(); ++it2)//n列 
		{
			if(*it2==target)
			{
				return pos;
			}
			++pos.second;
		}
		pos.second=1; 
		++pos.first;
	}
	return {-1,-1};
}
int main(){
	int n;cin>>n;
	vector<vector<int>>v(n+1);
	for(int i=0;i<v.size();++i)v[i].resize(n+1,0);//指定列数
	pair<int,int>p{0,0}; 
	
	v[1][n/2+1]=1;
//	p=vfind(1,v);cout<
	for(int i=2;i<n*n+1;++i){
		p=vfind(i-1,v);
		int x=p.first,y=p.second;
		if(x==-1)cout<<"ah";
		else{
			if(x==1&&y!=n)v[n][y+1]=i;			//在首行,但不在尾列 
			else if(x!=1&&y==n)v[x-1][1]=i;		//不在首行,但在尾列 
			else if(x==1&&y==n)v[x+1][y]=i;		//首行尾列 
			else if(x!=1&&y!=n){				//不在首行尾列 
//				cout<
				if(!v[x-1][y+1])v[x-1][y+1]=i;
				else v[x+1][y]=i;
			}
			else;
		}
	}
	for(int i=1;i<n+1;++i){
		for(int j=1;j<n+1;++j){
			cout<<v[i][j];
			if(j!=n)cout<<" ";
			if(j%n==0&&i!=n)cout<<endl;
		}
	}
	return 0;
}
/*
	还是很简单的一道题,刚好重温了一下STL的deque(double-ended queue),就用dq实现了,
	关于我的博客,如果你有任何不解的地方可私信我,
	我是学生党,技术有限,但是一定尽力解答,如果能帮助到你就是我的动力使然.
	*/

NOTE

关于几个比较重要的(个人做题过程中遇到的)一些一看就会一写就废的CodeBlock

//n维数组的从外到里的层次都是一个地址,分别为n-1层数组的首元素地址(至于这个元素是不是又是一个地址就看是不是最后一层),而最后一层才是数组元素值,不是地址.
//二维数组的行迭代器取值*仍然是一个地址(一维数组的首元素地址),且[这个地址.begin()/end()]表示首列,是一个数组元素值
//pair fx(int target,vector>arrname)定义了一个函数fx,返回类型是pair键值对,传参是目标值和二维数组名(vector name)
//pairpos{1,1};//初始化一个键值对,以二维数组遍历的起点(舍弃实际上的0行0列)开始
//vector>::iterator可以是auto
//auto it1=arr.begin()+1	遍历二维数组的行,并且main中我舍弃了首行首列,故逻辑上的首行(遍历起始行)应该是begin()+1
//vector::iterator也可以是auto
//auto it2=(*it1).begin()+1;首列
//it2!=(*it1).end();++it2因为步长是自增量1++,故!=end()不会越界,刚好遍历完(end()是最后元素的后一个元素)
//if(*it2==target)对列迭代器(最里层迭代器)取*则拿到我们需要的数组元素值了
//return pos返回键值对,并且这个键值对经过了偏移(++pos.second;++pos.first;),能正确反应在数组中的相对位置
//pos.second=1;二维数组的行遍历完,其列需要归位:因为行确实是遍历到几行找到了,这个行值就是坐标的x,而列是相对于行的,每一个行都有列,不能说我走完了一个五元素的行后,跳到下一行再走一步这个列就是元素的坐标y,(2,7)是错的,需要列归位,到下一行后,列值又从1开始,走一步到2,找到,返回(2,2)才是对的
//找不到的话就返回一个逻辑上不可能存在的坐标return {-1,-1};
}

最近事情挺多的,年底考研了,现在也得好好学一下数一英一了.
这个寒假看来也是加班加点的走进度了.
路好长啊,不过再坚持一年就能休息一下了.
如果你看到了这里,祝你天天开心ヽ(￣ω￣(￣ω￣〃)ゝ