codeforces round#524 D - Olya and magical square /// 大概算是数学规律题?
题目大意:
t 个测试用例 (1≤t≤103)
给定n k (1≤n≤10^9,1≤k≤10^18)
表示有一个边长为2^n的正方形格子 每次操作只能将一个格子切割为左上左下右上右下的四等分格子
问进行k次四等分切割后 能否使得 左下角的格子的边长 和 右上角的格子的边长 相等
并且存在一条左下角格子到右上角格子的路径上 经过的格子的边长 也和它们相等
若可以输出 “YES 切割后的log2(边长)” 若操作次数用不完输出“NO”
首先
对边长为2^1的格子四等分切割为1*1小格需要 1 次操作
对边长为2^2的格子四等分切割为1*1小格需要 1*4+1=5 次操作
对边长为2^3的格子四等分切割为1*1小格需要 5*4+1=21 次操作
......op[ i ] = op[ i-1 ] * 4 + 1
由此可预处理出边长为 2^i 的格子切割为1*1小格 需要 op[i] 次操作
然后我们可以发现当切割边长为2^31的格子时 op[31] 即操作次数超出了k的范围10^18
(k最大时不足以将一个2^31边长的格子切为1*1小格)
假设 边长为2^32时 先进行一次操作(分为四格2^31) 剩k-1次操作
之后只对右下的2^31的一格切割 那么k-1次操作绝对能用完
所以路径由左下经左上到右上 经过的三格的边长一样都是2^31 即输出log2(2^31)=31
以此类推 >31 的情况 只要这么处理 答案就是 n-1
n<=31时 考虑只切我们要走的路径的格子(假设我们走左边和上边的边缘圈的格子)
每次只对边缘圈的格子切割一次
第一次需要切割 1 格 (即切1次) (路径格子边长减为 2^(n-1) )
第二次需要切割 3 格 (即切3次) (路径格子边长减为 2^(n-2) )
第三次需要切割 7 格 (即切7次) (路径格子边长减为 2^(n-3) )
......
每次递推可得到下次需要切割的格子数 now(下次) = now(本次) * 2 + 1
累加得到边缘圈应切割次数 tot += now
但是仅仅只切割外围 k次操作很可能还是用不完的
那么此时我们考虑每次切割后不会成为外围圈的格子
因为它们不会影响到我们要走的路径 所以可以直接把它们切成1*1的小格
第一次红色格子可切割 共需切割次数 op[ n-1 ] * (3-2)
第二次绿色格子可切割 共需切割次数 op[ n-2 ] * (7-2)
第三次青色格子可切割 共需切割次数 op[ n-3 ] * (15-2)
.....(由于恰好对应下次切割要切割的外围圈格子往内的一圈 往内一圈会少两格 所以恰好是 now(下次)-2 格)
累加得到额外可切割次数 re
那么当只切外围圈的操作数 tot >= k 时 可得到答案
或者 当切外围也切内圈 tot+re>=k 时 也可得到答案
否则 k次操作 就不可能被用完
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll n,k,op[];
int main(){
for(int i=;i<=;i++) op[i]=op[i-]*4LL+1LL;
int t; scanf("%d",&t);
while(t--){
scanf("%I64d%I64d",&n,&k);
if(n>) {
printf("YES %I64d\n",n-);
continue;
}
ll tot=,now=,j=,re=;
while(now+tot<=k&&j<n){
tot+=now;
now=now*+; // now更新为下轮操作需要操作的边缘圈的格数
j++; // 对边缘圈的小格各操作一次 那么每格的边长又小了一半 即由2^(n-j)变为2^(n-(j+1))
re+=op[n-j]*(now-);
}
if(k>tot+re) printf("NO\n"); // 全部切到1*1小格的操作次数tot+re 仍然不够k次
else printf("YES %I64d\n",n-j); // n-j 即缩小到最后的 log2(边长)
}
return ;
}
codeforces round#524 D - Olya and magical square /// 大概算是数学规律题?的更多相关文章
- Codeforces Round #384 (Div. 2) B. Chloe and the sequence(规律题)
传送门 Description Chloe, the same as Vladik, is a competitive programmer. She didn't have any problems ...
- Codeforces Round #524 (Div. 2) D. Olya and magical square
D. Olya and magical square 题目链接:https://codeforces.com/contest/1080/problem/D 题意: 给出一个边长为2n的正方形,每次可以 ...
- 竞赛题解 - [CF 1080D]Olya and magical square
Olya and magical square - 竞赛题解 借鉴了一下神犇tly的博客QwQ(还是打一下广告) 终于弄懂了 Codeforces 传送门 『题目』(直接上翻译了) 给一个边长为 \( ...
- Codeforces Round #524 (Div. 2) Solution
A. Petya and Origami Water. #include <bits/stdc++.h> using namespace std; #define ll long long ...
- Codeforces Round #524 (Div. 2)(前三题题解)
这场比赛手速场+数学场,像我这样读题都读不大懂的蒟蒻表示呵呵呵. 第四题搞了半天,大概想出来了,但来不及(中途家里网炸了)查错,于是我交了两次丢了100分.幸亏这次没有掉rating. 比赛传送门:h ...
- Codeforces 715A & 716C Plus and Square Root【数学规律】 (Codeforces Round #372 (Div. 2))
C. Plus and Square Root time limit per test 2 seconds memory limit per test 256 megabytes input stan ...
- [Codeforces Round #247 (Div. 2)] A. Black Square
A. Black Square time limit per test 1 second memory limit per test 256 megabytes input standard inpu ...
- Codeforces Round #524 (Div. 2) F. Katya and Segments Sets(主席树)
https://codeforces.com/contest/1080/problem/F 题意 有k个区间,区间的种类有n种,有m个询问(n,m<=1e5,k<=3e5),每次询问a,b ...
- Codeforces Round #524 (Div. 2) E. Sonya and Matrix Beauty(字符串哈希,马拉车)
https://codeforces.com/contest/1080/problem/E 题意 有一个n*m(<=250)的字符矩阵,对于每个子矩阵的每一行可以任意交换字符的顺序,使得每一行每 ...
随机推荐
- A. Srdce and Triangle--“今日头条杯”首届湖北省大学程序设计竞赛(网络同步赛)
如下图这是“今日头条杯”首届湖北省大学程序设计竞赛的第一题,作为赛后补题 题目描述:链接点此 这套题的github地址(里面包含了数据,题解,现场排名):点此 Let be a regualr tr ...
- Codeforces 1176B - Merge it!
题目链接:http://codeforces.com/problemset/problem/1176/B 题意:给定序列,任意俩个元素可以相加成一个元素,求序列元素能被3整除的最大数量. 思路: 对于 ...
- 专题:『systemd』
[COMMAND] Systemd Daemons: systemd networkd journald logind "user session" Systemd units: ...
- zip压缩详细分析
该文章转自:http://www.cnblogs.com/esingchan/p/3958962.html (文章写得很详细,让我对zip压缩有了了解,感谢博主,贴在这是为了防止忘了有这么好的文章,侵 ...
- python中正则匹配之re模块
Python中正则表达式 re:re是提供正则表达式匹配操作的模块 一.什么是正则表达式 正则表达式是一个特殊的字符序列,它能帮助你方便的检查一个字符串是否与某个模式匹配,Python 自1.5版本起 ...
- Zabbix监控搭建
目录 Zabbix概述 zabbix是一个基于WEB界面的提供分布式系统监视以及网络监视功能的企业级的开源解决方案 ( 基于 GPL V2 )zabbix由 2 部分构成,zabbix ...
- 新项目UX设计0到1的正确开启方式
无论是在BAT还是创业小公司,都随时可能接到从0开始的新项目,那么作为负责新项目的主设OR独立设计师,我们应该从何开启工作呢?
- Peasy.NET学习之并发问题处理
Peasy.net之并发处理 BusinessServiceBase是ServiceBase的自定义实现,提供了额外的独特功能 首先,创建一个业务服务,该业务服务必须继承BusinessService ...
- Visio2016专业版永久激活码
Visio2016专业版永久激活码: [Key]:NKVJM-8MTT4-8YDFR-6738M-DPFJH [Key]:W9WC2-JN9W2-H4CBV-24QR7-M4HB8 [Key]:7K8 ...
- [转]Netty入门(最简单的Netty客户端/服务器程序)
Java中的NIO是一种解决阻塞式IO问题的基本技术,但是NIO的编写对java程序员是有比较高的要求的.那么Netty就是一种简化操作的一个成熟的网络IO编程框架.这里简单介绍一个程序,代码是< ...