[AGC025D]Choosing Points

题目大意:

输⼊\(n(n\le300),d_1,d_2\),你要找到\(n^2\)个整点\((x,y)\)满⾜\(0\le x,y<2n\)。并且找到的任意两个点距离,既不是\(\sqrt{d_1}\),也不是\(\sqrt{d_2}\)。

思路:

所有距离为\(\sqrt{d_1}\)的点连边,可以得到一个⼆分图。\(d_2\)同理。注意到满足\(a^2+b^2=d\)的\((a,b)\)只有\(\mathcal O(n)\)个,可以得到\(\mathcal O(n^3)\)的二分图染色做法。

源代码:

#include<cmath>
#include<cstdio>
#include<cctype>
#include<vector>
inline int getint() {
register char ch;
while(!isdigit(ch=getchar()));
register int x=ch^'0';
while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
return x;
}
const int N=600;
int n,d1,d2,col[2][N*N];
std::vector<std::pair<int,int>> v[2];
inline int id(const int &x,const int &y) {
return x*N+y;
}
inline bool check(const int &x) {
return 0<=x&&x<n*2;
}
void dfs(const int &x,const int &t) {
const int i=x/N,j=x%N;
for(auto &d:v[t]) {
const int nx=i+d.first,ny=j+d.second;
if(!check(nx)||!check(ny)) continue;
const int y=id(i+d.first,j+d.second);
if(!col[t][y]) {
col[t][y]=col[t][x]==1?2:1;
dfs(y,t);
}
}
}
int main() {
n=getint(),d1=getint(),d2=getint();
for(register int i=0;i<=d1;i++) {
const int j=sqrt(d1-i*i);
if(i*i+j*j!=d1) continue;
v[0].emplace_back(std::make_pair(i,j));
v[0].emplace_back(std::make_pair(-i,j));
v[0].emplace_back(std::make_pair(i,-j));
v[0].emplace_back(std::make_pair(-i,-j));
}
for(register int i=0;i<=d2;i++) {
const int j=sqrt(d2-i*i);
if(i*i+j*j!=d2) continue;
v[1].emplace_back(std::make_pair(i,j));
v[1].emplace_back(std::make_pair(-i,j));
v[1].emplace_back(std::make_pair(i,-j));
v[1].emplace_back(std::make_pair(-i,-j));
}
for(register int i=0;i<n*2;i++) {
for(register int j=0;j<n*2;j++) {
if(!col[0][id(i,j)]) {
col[0][id(i,j)]=1;
dfs(id(i,j),0);
}
if(!col[1][id(i,j)]) {
col[1][id(i,j)]=1;
dfs(id(i,j),1);
}
}
}
for(register int i=0,cnt=0;i<n*2;i++) {
for(register int j=0;j<n*2;j++) {
if(col[0][id(i,j)]==1&&col[1][id(i,j)]==1) {
printf("%d %d\n",i,j);
if(++cnt==n*n) return 0;
}
}
}
}

最新文章

  1. DIV+CSS设计IE6浮动产生双倍距离
  2. centos6 下安装xfce+vnc
  3. Js中把JSON字符串转换为JSON对象(eval()、new Function())
  4. 浅谈 Python 的 with 语句
  5. Code First 更新数据库结构
  6. C# Winform窗口之间传值的多种方法浅析(转)
  7. 用phpcms如何将静态页面制作成企业网站,头部加尾部
  8. Oracle中用户(User)和模式(Schema)的概念
  9. python/匿名函数和内置函数
  10. jQuery的版本兼容问题
  11. ASP 基础一 网站开发 初步认识
  12. kaggle PredictingRedHatBusinessValue 简单的xgboost的交叉验证
  13. Java 面向切面 AOP
  14. 使用iTextSharp修改PDF文件(一)
  15. Linux 客户端bind函数的使用
  16. Monkey学习笔记&lt;五&gt;:检查内存泄露
  17. linux中配置JDK环境变量
  18. IOS8 Playground介绍
  19. oracle 分区表详解
  20. 在Web API 2 中实现带JSON的Patch请求

热门文章

  1. CentOS7安装MySQL5.7以及修改密码
  2. netlink socket编程
  3. tex src
  4. rhel-server srpms iso
  5. mysql连接池优化笔记
  6. 1000: 恶意IP 课程作业
  7. 获取windows 网卡GUID和ip信息
  8. [How to]如何通过xib来自定义UIViewController
  9. Vim常用命令(转)—默写版
  10. Log4Net中配置文件的解释