AGC027 D - Modulo Matrix 构造
目录
题目链接
题解
从第左上角第一个点开始染色,相邻不同色,染法唯一
那么一个点的四周与他不同色,我们另这个点比四周都大,那么这个点权值可以使lcm(四周的点权值)+1
于是我们就得到了一种构造方案,染色后对一种颜色的点进行赋值,然后另一种颜色的点取lcm
可是....直接这样瞎构造会爆掉1e15
对于一种染色点,可以按照i + j和i - j分为两类,每一类乘上一个相同的质数
对于当前格子的价值就是从左上角到右下角,和从右上角到左下角穿过他的第k素数乘积
这样构造的lcm最大就是四个素数(prime[n],prime[n - 1],prime[2 * n],prime[2 * n - 1])的乘积
不会超过1e15
代码
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define gc getchar()
#define pc putchar
#define LL long long
inline int read() {
int x = 0,f = 1;
char c = getchar();
while(c < '0' || c > '9') c = gc;
while(c <= '9' && c >= '0') x = x * 10 + c - '0',c = getchar();
return x * f;
}
void print(LL x) {
if(x < 0) {
pc('-');
x = -x;
}
if(x >= 10) print(x / 10);
pc(x % 10 + '0');
}
const int maxn = 20007;
int n;
int prime[maxn];
bool vis[maxn];
void pre(int lim = 10000) {
for(int cnt = 0,i = 2;i <= lim;++ i) {
if(!vis[i]) prime[++ cnt] = i;
for(int j = 1;j <= cnt && prime[j] * i <= lim;++ j) {
vis[i * prime[j]] = 1;
if(i % prime[j] == 0) break;
}
}
}
LL a[507][507];
LL gcd(LL x,LL y) {
return y == 0 ? x : gcd(y,x % y);
}
LL lcm(LL x,LL y) {
return x / gcd(x,y) * y;
}
int main() {
n = read();
pre();
if(n == 2) {
pc('4');pc(' ');pc('7');pc('\n');
print(23);pc(' '); print(10); pc('\n');
return 0;
}
for(int i = 0;i <= n + 1;++ i)
for(int j = 0;j <= n + 1;++ j)
a[i][j] = 1;
for(int i = 1;i <= n;++ i)
for(int j = 1;j <= n;++ j) {
if((j & 1) == (i & 1)) {
a[i][j] = prime[(i + j) / 2] * prime[(i + n + 1 - j) / 2 + n];
a[i + 1][j] = lcm(a[i + 1][j],a[i][j]);
a[i - 1][j] = lcm(a[i - 1][j],a[i][j]);
a[i][j + 1] = lcm(a[i][j + 1],a[i][j]);
a[i][j - 1] = lcm(a[i][j - 1],a[i][j]);
}
}
for(int i = 1;i <= n;++ i)
for(int j = 1;j <= n;++ j)
if((i & 1) == (j & 1));
else a[i][j] ++;
for(int i = 1;i <= n;++ i,pc('\n'))
for(int j = 1;j <= n;++ j)
print(a[i][j]),pc(' ');
return 0;
}
AGC027 D - Modulo Matrix 构造的更多相关文章
- AtCoder Grand Contest 027 (AGC017) D - Modulo Matrix 构造
原文链接https://www.cnblogs.com/zhouzhendong/p/AGC027C.html 题解 首先我们假装 max mod min = 1 然后对着这个构造. 将各自黑白染色, ...
- AGC 027D.Modulo Matrix(构造 黑白染色)
题目链接 \(Description\) 给定\(n\),要求构造一个\(n\times n\)的矩阵,矩阵内的元素两两不同,且任意相邻的两个元素\(x,y\),满足\(\max(x,y)\ \mat ...
- agc027D - Modulo Matrix(构造 黑白染色)
题意 题目链接 构造一个\(n * n\)的矩阵,要求任意相邻的两个数\(a,b\),使得\(max(a,b) \% min(a,b) \not = 0\) Sol 我的思路: 假设\(mod = 1 ...
- 「AGC027D」Modulo Matrix
「AGC027D」Modulo Matrix 传送门 神仙构造题. 首先考虑一个非常自然的思路,我们把棋盘黑白染色后会变成一个二分图,黑色棋子只会与白色棋子相邻. 也就是说,我们可以将二分图的一部随便 ...
- hdu 5015 233 Matrix(构造矩阵)
http://acm.hdu.edu.cn/showproblem.php?pid=5015 由于是个二维的递推式,当时没有想到能够这样构造矩阵.从列上看,当前这一列都是由前一列递推得到.依据这一点来 ...
- Educational Codeforces Round 48 (Rated for Div. 2) D 1016D Vasya And The Matrix (构造)
D. Vasya And The Matrix time limit per test 2 seconds memory limit per test 256 megabytes input stan ...
- UVA 11149 Power of Matrix 构造矩阵
题目大意:意思就是让求A(A是矩阵)+A2+A3+A4+A5+A6+······+AK,其中矩阵范围n<=40,k<=1000000. 解题思路:由于k的取值范围很大,所以很自然地想到了二 ...
- CF486B OR in Matrix(构造+思维)
CF486B 一道有趣的思维题 由于or的性质可知只要a[i][j]为1那么b中第i行,第j列将都变成1 相反的,如果b[i][j]是0那么a中第i行,第j列都必须是0 根据第二个性质我们可以构造出a ...
- [atAGC027D]Modulo Matrix
对网格图黑白染色,在黑色格中填不同的质数,白色格中填相邻黑色格的lcm+1,但这样会超过1e15的上限将网格图划分为两类对角线,每一条对角线选一个质数,然后每一个点就是两条对角线的质数相乘,而白格的值 ...
随机推荐
- SpringBoot整合Email(电子邮件服务)
(1).导入starter依赖 <dependency> <groupId>org.springframework.boot</groupId> <artif ...
- Ajax+json+jquery实现无限瀑布流布局
<!doctype html> <html> <head> <meta charset="utf-8"> <title> ...
- eclipse自动添加注释
自动添加注释 快捷键:alt shift jwindows-->preference Java-->Code Style-->Code Templates code-->new ...
- yum install oracle-validated
背景 当时心血来潮要在linux搞oracle,可一顿折腾,大约两个周时间,主要是各种环境的检测麻烦,在redhat上操作也不如centos有利. 命令 yum install oracle-vali ...
- javascript判断是用什么设备打开
var userAgentInfo = navigator.userAgent //查看浏览器用于 HTTP 请求的用户代理头的值 var agents = ["Android", ...
- CBAM: 卷积块注意模块
CBAM: Convolutional Block Attention Module 论文地址:https://arxiv.org/abs/1807.06521 简介:我们提出了卷积块注意模块 ( ...
- OneNET麒麟座应用开发之十:空气质量数据监测站项目总结
大气质量数据监测站用于测试空气质量监测及数据采集,实现野外或者室内空气质量的检测. 1.项目概述 本项目是一个定制项目,要求采集大气的压力.温度.湿度.PM25.位置等数据并上传到指定的后台服务器.但 ...
- STM32F412应用开发笔记之九:移植FreeRTOS到F412ZG平台
在开发实际应用系统时,我们经常需要考虑数据的实时性和多任务,嵌入式实时操作系统的出现为实现这一目的提供了很好的助力.FreeRTOS是近年来比较流行的嵌入式实时操作系统,而且是开源免费的,STM32C ...
- poj2185 kmp求最小覆盖矩阵,好题!
/* 特征值k=m-next[m]就是最小循环节的长度, m%k就是去末尾遗留长度 */ #include<iostream> #include<cstring> #inclu ...
- pytest十一:函数传参和 firture 传参数 request
为了提高代码的复用性,我们在写用例的时候,会用到函数,然后不同的用例去调用这个函数.比如登录操作,大部分的用例都会先登录,那就需要把登录单独抽出来写个函数,其它用例全部的调用这个登录函数就行.但是登录 ...