Ruotian's I/O

2023年终总结

Sun, 31 Dec 2023 08:11:06 +0800

2023 结束了，简单记录一下。

🎓 研究生毕业：度过了有惊无险的毕业季，从写论文、预答辩，到最终答辩，最后参加毕业典礼。离开了熟悉的校园，离开了熟悉的课题组老师、同学。

九寨沟之行

Mon, 01 May 2023 16:29:34 +0800

五一和超超去了九寨沟玩，风景非常好。

Figure 1

Figure 2

Figure 3

2022年终总结

Sat, 31 Dec 2022 23:59:59 +0800

2022有点累，不想总结了。2023打算躺了。

Elden Ring

Sun, 18 Dec 2022 20:55:25 +0800

玩了一会，有点难。

Elden Ring

打过了第一个 Boss。

恶兆妖鬼玛尔基特

谷歌实习初体验

Sun, 24 Jul 2022 16:34:35 +0800

很棒！

归档

Fri, 27 May 2022 11:27:36 +0800

💬 Chat

Wed, 25 May 2022 16:50:17 +0800

End-to-end encrypted chat.
基于 Firebase Realtime Database.

Peers

(loading...) online

Message

Speak

(your id):

(ctrl or shift + enter)

Action

board id: connecting...

danger

Source code

https://github.com/zrt/chat

绕过清华校园网的二层隔离

Mon, 16 May 2022 15:45:45 +0800

写了个脚本绕过校园网二层隔离，自动化这个教程里面提到的加路由的方法。

依赖 ipcalc 计算子网地址。apt install ipcalc。

#!/bin/bash

IPA=$(ip address list dev eth0 | awk ' /inet/ {print $2}' | head -1)
GATE=$(ip route list dev eth0 | awk ' /^default/ {print $3}')

IP=$(echo $IPA | awk -F'/' '{print $1}')
MASK=$(echo $IPA | awk -F'/' '{print $2}')
NETWORK=$(ipcalc $IPA -n -b | grep Network | awk '{print $2}')

echo =========================
echo Your IP: $IP
echo Netmask: $MASK
echo Gateway: $GATE

echo adding route: ip r a $NETWORK via $GATE
sudo ip r a $NETWORK via $GATE
sudo ip r a $GATE dev eth0
echo ok

PyScript 尝试

Wed, 11 May 2022 21:04:03 +0800

- numpy - matplotlib

PyScript让你可以在 html 里写 Python ，感觉很有趣。

比如下面这个图（大概要等5分钟😂）：

  <div id="plot"></div>
  <py-script output="plot">
  import matplotlib.pyplot as plt
  import numpy as np

  x = np.random.randn(1000)
  y = np.random.randn(1000)

  fig, ax = plt.subplots()
  ax.scatter(x, y)
  fig
  </py-script>

REPL

Python 协程 (async/await)

Tue, 10 May 2022 13:47:48 +0800

用Telegram 机器人库 Telethon的时候发现他家的 API 接口走的都是 async。

这就来学一下 Python 的 async 和 await (官方文档)，需要 Python 版本3.7以上。

import asyncio

async def f():
  print('start')
  await asyncio.sleep(1) # sleep 1s
  print('end')

asyncio.run(f())

f() 返回了一个协程对象，asyncio.run调用会执行这个对象，输出start和end之间会等一秒钟。 async和await的语法类似js。

大端序与小端序的C++存取

Sun, 08 May 2022 17:52:53 +0800

把int等数字转换成字节流，需要考虑字节的顺序，有大端序或小端序两种。这两种顺序很多平台、算法、协议里都有用到。

大端序

int数字 0xAABBCCDD按照大端序存储，顺序应该是0xAA, 0xBB, 0xCC, 0xDD。

关于

Sat, 07 May 2022 00:47:27 +0800

Telegram channel: @chanstring
GitHub: @zrt
Mail: admin@这个网站的域名
投资专栏：investing.ruotian.io

Python 使用子进程调用函数并得到 exit code

Tue, 03 May 2022 18:16:31 +0800

Python 的 subprocess 并不支持直接调用函数。可以用 multiprocessing 库来搞。

def twosum(a, b):
    exit(a+b)

from multiprocessing import Process

if __name__ == '__main__':
    p = Process(target=twosum, args=(2, 33))
    p.start()
    p.join()
    print('exit code:', p.exitcode)

如果需要传输其他类型的变量，可以用 multiprocessing 里的 Queue 。

Scrooge

Sat, 26 Mar 2022 13:41:03 +0800

买了一本漫画书，《史高治的光辉岁月》。

🧶 Unravel Two

Tue, 22 Mar 2022 13:38:36 +0800

《Unravel Two》

创建一个 Service/Systemctl 服务

Fri, 04 Mar 2022 18:12:45 +0800

写入{taskname}.service文件到/etc/systemd/system/文件夹。
sudo systemctl daemon-reload
sudo systemctl enable {taskname}.service
sudo systemctl start {taskname}.service
systemctl status {taskname}.service查看状态。

{taskname}.service 例子:

使用有线网络连接树莓派

Sat, 19 Feb 2022 13:55:13 +0800

使用此方法需要电脑上有网线接口，或者自备USB网卡。
使用 USB-C 电源线连接树莓派，树莓派会自动开机。
拿一根网线连接树莓派和电脑，可以看到网口的状态灯亮起。
开启网络共享，使得树莓派被分配 IP 地址并且可以通过电脑上网。
1. Windows 参考资料：Link1, Link2
2. Linux 参考资料：Link
3. macOS 参考资料：Link1, Link2
4. 需要注意的是，如果共享的网络使用 802.1X 认证（例如 Tsinghua-Secure 和 eduroam），可能会共享失败，请切换到其他网络（例如Tsinghua/Tsinghua-5G/手机热点）。
找到树莓派的 IP 地址
1. 【推荐方法】在命令行终端内尝试 ping sysintro-{ID}.local得到树莓派的 IP 地址({ID}替换为SD卡上的编号，如sysintro-23.local)。
2. 使用arp -a命令列出各个网口上通过 ARP 协议发现的设备的 IP 地址，找到其中树莓派的地址。
  1. 可以通过网络接口 (interface) 的名称/MAC地址，并结合ipconfig, ifconfig的信息找到有线网卡或USB网卡对应的 ARP 记录。
使用 ssh 命令/或者 VS Code 尝试连接树莓派，例如 ssh pi@192.168.2.2，其中192.168.2.2需替换为上一步找到的树莓派 IP 地址。密码为 SysIntro2022!。

参考：https://lab.cs.tsinghua.edu.cn/router/doc/howto/raspberry_pi/

树莓派SD卡空间扩容

Sat, 19 Feb 2022 13:50:03 +0800

使用 ssh 连接到 RPi 上的终端。
使用 df -h 查看当前空间情况，Mounted on / 为 SD 卡空间，扩容前为 4GB左右。
使用 sudo fdisk /dev/mmcblk0 进入fdisk 的操作界面，按p回车，大概如图：

Apache 反向代理使用 Basic Auth

Sat, 05 Feb 2022 18:10:11 +0800

生成密码 sudo htpasswd -c /etc/apache2/.htpasswd {username}
修改 /etc/apache2/sites-available/{site}.conf 文件，添加:

<Proxy *>
        Order deny,allow
        Allow from all
        Authtype Basic
        Authname "Password Required"
        AuthUserFile /etc/apache2/.htpasswd
        Require valid-user
</Proxy>

3. 重启 Apache 服务 sudo systemctl restart apache2

购入 Xbox Seires X

Thu, 27 Jan 2022 14:03:45 +0800

溢价购入 Xbox Series X，体验很棒！

2021年终总结

Fri, 31 Dec 2021 08:00:00 +0800

2021 结束了，记录一下一年来发生的有趣的事情。

奇思妙想

Wed, 15 Dec 2021 13:08:47 +0800

这里记录一些奇思妙想。

投掷计算器

Tue, 14 Dec 2021 21:33:42 +0800

皮卡丘可以把一个篮球扔出去多远?

Sublime Text 开启 GPU 渲染

Tue, 14 Dec 2021 18:05:33 +0800

Sublime Text 4 中加入了 GPU 渲染，不过默认没有开启。可以通过设置开启一下，渲染速度飞快。

Steve Jobs的经典演讲

Thu, 09 Dec 2021 00:48:29 +0800

Squid Game

Mon, 06 Dec 2021 17:37:54 +0800

在Netflix上看完了《鱿鱼游戏》(Squid Game)。

这个剧的音效很棒，可以很好的感受人物的情绪变化。
一些设计很棒，或者说剧中把“一些人很用心的设计了整个游戏”这个理念传达的很到位。：

Sony电视 Setup

Mon, 29 Nov 2021 10:41:43 +0000

收货
安装
- 打孔
- 挂起来
校准声音
校准图像
用u盘安装Clash小猫咪，配置网络
安装Aptoide应用商店
安装Netflix, Youtube, Apple TV, 旧版Bilibili, Kodi等应用
- Youtube用FireTV版本因为国行没google框架
- Bilibili 1.6.4目前可用
连接ADB更换launcher, 卸载默认瀑布流启动器，新启动器用的“leanback launcher fire tv”
连接Switch

感谢亚马逊的Fire tv，给没有google框架的国行机器创造了无限可能。

BWT

Fri, 08 Oct 2021 05:42:21 +0000

BWT: Burrows-Wheeler transform (aka. 块排序压缩)。

一篇 Genome Biology 2009 的论文

《Ultrafast and memory-efficient alignment of short DNA sequences to the human genome》(Bowtie)

这篇文章使用BWT对人类基因组建立索引，并用贪心+回溯支持模糊匹配。

滑板信用卡🛹

Sat, 28 Aug 2021 14:00:17 +0000

交通银行出品
美国运通卡组织，（滑板带闪付，但没地方刷）
无工本费，年费刷免
GSC制卡
还赠一个台灯盒子

macOS使用Secoclient VPN时ssh, scp卡死的解决办法

Tue, 10 Aug 2021 04:01:41 +0000

这个issue指出这个问题主要是因为mtu太大了。

运行此命令可以缓解：

#!/bin/bash
ip link set dev $(ip r get 192.168.1.238 | awk '{print $5}') mtu 750

Hello, macOS Big Sur!

Thu, 05 Aug 2021 08:43:53 +0000

MacBook Pro 16′ 2019更新了Big Sur。

抹盘安装 + 重置SMC。目前无卡顿。

使用体验

这个UI画风很喜欢。
字体变小了，任务栏变窄了。
风扇控制比Catalina要好一些。
TimeMachine支持APFS的快照，安全感++。
没有发现软件不兼容的情况。

Issues

联想ThinkPad E480 CPU降频修复

Mon, 02 Aug 2021 14:14:02 +0000

i7-8550U CPU锁频0.9GHZ，调整电源选项无效。BIOS中关闭节能降频后可以到1.8GHZ，但无法睿频，CPU-Z中查看CPU电压较低。

CPU 0.79GHz运行

解决方法：刷BIOS降到1.19解决。

重庆之行

Sun, 01 Aug 2021 07:21:07 +0000

和超超去了重庆玩～

使用Sony A6000+Sigma 30mm f/1.4拍摄。

修Macbook Pro 16' (2)

Sun, 01 Aug 2021 06:22:50 +0000

7.19 晚，打开微信，一个学姐让我推荐笔记本电脑，我想了一下说了“m1 Macbook”。

接着我从我的包里拿出我的Macbook，发现按开机键没有反应，插电没有反应。

清华大学学费代交服务(2021)

Mon, 12 Jul 2021 15:29:44 +0000

今年活动继续，代交学费立减1%！

（¥5000学费立省¥50，一顿香锅钱就出来了~）

名额有限，欢迎留言、微信咨询预约~

温馨提示：学校8.12号开启交费通道，如果想参与代交学费活动注意学校绑定的银行卡里不要放钱，避免自动扣款。

☁️

Thu, 08 Jul 2021 09:21:13 +0000

清华大学美术学院毕业作品展

Mon, 05 Jul 2021 17:47:33 +0000

向多样的世界提问

去看了美术学院在美院和艺术博物馆展出的研究生&本科生毕设展。整体很不错~

图片点击不能放大。

研究生

本科生

官方还有在线展，有兴趣的可以去找一找，这里就不放链接了。

王兴已经两个月没发过饭否了

Mon, 05 Jul 2021 22:15:01 +0800

如果我一整天都没看到、想到、或做过什么值得在饭否上说的事，那这一天就太浑浑噩噩了。

自述

前两天一首唐诗让我颇为感叹：秦国警惕书生，但后来灭秦的刘项却不读书。这给我提了个醒，最危险的对手往往不是预料中的那些。这些年阿里一直盯着京东，最后却是拼多多斜刺里杀出来，用户数一举超过淘宝。同理，美团外卖最大的对手看起来是饿了么，但更可能颠覆外卖的却是我们还没关注到的公司和模式。

清洗HHKB键盘

Wed, 19 May 2021 21:59:32 +0800

把键盘拆下来洗了洗，因为是无刻的，所以需要记住每个键在哪一行。

小程序「THU洗衣」开通“流量主”

Thu, 13 May 2021 13:38:03 +0000

小程序够累积1000uv可以开启“流量主”广告服务。

前段时间洗衣提醒小程序恰好够了，就加上了广告，主要是想体验一下微信广告的API。

微信广告“流量主”的说明文档

扫码直达广告页面：

利用Service Worker实现CDN

Fri, 09 Apr 2021 22:02:51 +0800

Service worker(详细介绍)是一个浏览器 API，它可以拦截并修改网页所有发出的网络请求。

利用这一点可以将所有需要加速的网络请求的 URL 修改为 CDN 地址实现加速。

以下service-worker.js代码使用监听fetch事件的方法拦截并修改所有的网络请求。

无锡之行

Sun, 04 Apr 2021 01:02:49 +0800

和超超去了无锡玩~

小米对讲机1s

Sat, 27 Feb 2021 18:21:33 +0000

最近入手了小米对讲机1s，终于有了不依赖电信基础设施的局域通信信道(~~黑电台~~)。

长这个样子

点阵LED显示。
支持3.5mm耳机、蓝牙耳机。
支持409.7500MHz~409.9875MHz的20个公共频道(L)，0.5W功率发射。
支持430.000MHz~440.000MHz的20个远距离频道(H)或自定义频道，3W发射。（需要业余无线电操作证才能合法发射）
支持FM收音机。

实际测试能够清晰通信的极限距离：

如何在清华吃汤圆？

Thu, 25 Feb 2021 17:28:46 +0800

清华吃汤圆，三步走策略

照澜院超市购买汤圆
清青小火锅点一锅火锅
偷偷在火锅中煮汤圆

小程序：THU洗衣

Mon, 22 Feb 2021 08:42:29 +0000

开学了，发个广告~

【洗衣机空闲提醒小程序】

源代码：https://github.com/zrt/thu-xiyi

2021.2.22：用户超过100。
2021.2.23：用户超过200。
2021.2.26：用户超过300。

工商银行Clear卡

Sat, 23 Jan 2021 16:48:22 +0000

透明卡片
美国运通
官方链接
到2021年底，每月送1杯奶茶

早安喵

Mon, 18 Jan 2021 11:00:17 +0800

想记录一下和卿卿在一起的心动瞬间~

上个月买了一个蓝色相片册，用来收集一起玩耍的证据。卿卿好像猜到了我在收集，每次都会主动把票给我。
卿卿突然送了我一个bilibili的兔子围巾，好可爱啊。想为这个围巾买一件可爱的衣服来配，但是为什么找不到我想要的款式呢，要是衣服可以3d打印就好了。

宜家鲨鱼洗澡记

Thu, 07 Jan 2021 10:46:19 +0800

今天的主角是这只鲨鲨~

但是它最近脏了，怎么洗呢？

从官网可以了解到，它是可以机洗的。

于是把鲨鲨扔进洗衣机里hiahiahia

可能是受到了旋转的离心力影响，出浴时鲨鲨内部填充物都集中到了肚子，变成了生无可恋的球鲨。

2020年终总结

Thu, 31 Dec 2020 08:00:00 +0800

TL;DR 拥抱变化

经常从阿里相关新闻听到“拥抱变化”，2020发生的变化很多，我也在学着适应许多变化。
今年年初开端就是疫情，过年没有串门拜年，买了很多口罩。
紧接着，疫情改变了一切。
本科的最后一学期在家度过，完成了两个毕设，毕业成功。
假期去了躺海南，肉眼看了一次火箭发射。
8月底研究生线下入学成功，开启了全新的研究生生涯。
因为新生舞会认识了YC同学。和YC同学在一起啦。和YC同学到处玩玩玩玩玩，然后就封校了。
然后就到了2020的结尾。

本科阶段的结束

寒假放假时大家谁也没有想到之后就很难再见到了（或许mingming想到了，因为他把电脑搬回家了）。
接着疫情越来越严重，记着当时微博40个热搜都是疫情的时候，我和爸爸去社区买了一些口罩，紧接着就变的严峻起来。

《光遇》里那些好看的截图

Thu, 17 Dec 2020 10:55:31 +0800

脱单啦

Sat, 21 Nov 2020 14:00:19 +0000

🥰🎉❤️🥳

今天下了今年的第一场雪。

下午在雪后的校园中一起走了很久。

晚上吃完夜宵在桃李门口表白成功啦！❤️

然后深夜一起看了电影《海兽之子》。

2020初雪

国庆中秋，双节同庆

Thu, 01 Oct 2020 04:15:23 +0000

今年中秋赶上了国庆，虽然只放4天假，但还是要好好庆祝一下~

彩虹

Thu, 03 Sep 2020 14:12:45 +0000

第一次在清华观测到双层彩虹。

洗衣机状态

Sun, 30 Aug 2020 15:19:08 +0000

由于洗衣机厂商升级，本文已经失效

获取中…

基于Cloudflare workers。

Zenithal同学对这个接口很感兴趣，做成了一个开放式的查询接口。
接口的形式大概是这个样子的`https://washer.zenithal.workers.dev/?s=二号楼二层`，可以返回json结果，文档请参考这里。

修Macbook Pro 16'

Fri, 28 Aug 2020 12:24:43 +0000

今天去Apple 王府井修电脑(Macbook Pro 16′)。

主要有两个问题，一个是屏幕底下的位置的玻璃碎了，另外一个问题是方向键被上次拆坏了（上次可乐撒到了键盘上，拆下来清洗来着）。

开学典礼

Thu, 27 Aug 2020 07:00:46 +0000

8.26是开学典礼的日子。

上午成功举办了开学典礼，开学典礼后开了第一次班会。

成功选出了班委。

下午搬家。

另， https://mp.weixin.qq.com/s/8rVFL3tzXEU1y-o5XYm6wQ

新生报到

Tue, 25 Aug 2020 15:28:15 +0000

今天8.25是新生报到的日子。

早上坐火车去首都北京。
打车从北京西到清华东南门。
东南门检查录取通知书、核酸证明、绿码批准进校。
FIT 1区2层报到，详见 https://mp.weixin.qq.com/s/EerbExNXX2kCVP5xVWSfqQ
带着行李坐校车去宿舍楼报到。
抢救自行车。
圈存校园卡，水卡。
去超市购物。
晚上和同学出去去五道口吃饭，七夕。
逛卜蜂莲花超市。
刷学生卡进校。
收拾新宿舍。

清华大学学费代交服务（2020）

Sat, 15 Aug 2020 08:09:01 +0000

限时99折特惠~ (RMB 5000学费立省¥50，一顿外卖、奶茶钱就出来了~)

学校确认后再转钱，无任何风险。

心动不如行动，立即打开 info-综合-财务 -网上交费 – 交费 – 微信/支付宝，把二维码拍照发给我来交学费吧~

观看了一次火箭发射

Tue, 11 Aug 2020 01:02:59 +0800

UOJ多评测机一键配置脚本

Sat, 08 Aug 2020 00:40:00 +0000

OI赛制比赛的终测，一个评测机显然是不够的。200人的比赛，要10台评测机，7分钟才能测完。（成本大概10块钱？）

但是这些评测机平时闲置着又很浪费，于是我参考uoj开源版文档，写了一个一键配置脚本。配合阿里云的机器hostname自动顺序命名来使用。

RSA来了，又走了

Fri, 07 Aug 2020 15:14:03 +0000

没选择PGP、GPG，选择了裸的RSA。

#303 -> #36896

建立起了信任体系，又毁掉了。

美国运通人民币信用卡

Thu, 06 Aug 2020 14:20:00 +0000

美国运通今天0点在官网上发布了《携手启程一起通往精彩无数》一文，紧接着国内四家商业银行浦发、平安、民生、广发发布了美国运通单标人民币信用卡。

在看到广发的广告后随手申请了一张“广发美国运通生活+卡”，然后瞬间就通过了（卡号还挺吉利）。这张卡似乎定位挺普通的，没有机场、高铁权益，没有积分、里程兑换，线下刷卡完全没有好处（不过线下应该也没地方刷）。线上网购宣称返现10%，不过其规则极其复杂，简单疏理了一下，不考虑分期的话，大概一个月消费20000以内是无脑1%返现，还算可以吧，可以用来交学费(雾)。

Apple Watch 使用初体验

Wed, 05 Aug 2020 11:33:04 +0000

Apple Watch Series 5 使用了两个星期，来谈谈感受。

能看时间。
能看手机通知。
能刷地铁（手机也行，而且比左手戴表要方便）。
能用Siri。
能测心电图。
watch上的众多第三方app都没什么卵用。
那三个圆圈确实能刺激去运动。
一天一充还是有一点点麻烦的。

冲动消费不可取。

大学四年回顾

Thu, 11 Jun 2020 12:46:35 +0000

长图(13MB)

Arduino连接两个MPU6050

Thu, 09 Apr 2020 12:45:00 +0000

关于 MPU 6050

MPU 6050是一个陀螺仪+加速度仪传感器芯片，通过I2C总线的方式与Arduino通信。

传感器还包含一个DMP(Digital Motion Processor)模块，能够对传感器的数据进行一些处理。

Codeforces April Fools Day Contest 2020

Wed, 01 Apr 2020 09:41:00 +0000

愚人节专场，打一乐。官方tutorial

A. Is it rated?

puts("no");

B. Limericks

看样例35->57, 57->319, 391->1723，敏锐的觉察到可能是质因数分解。

for(int i=2;i<=x;i++){
        while(x%i==0){
            printf("%d",i);
            x/=i;
        }
    }

拒绝推荐系统

Thu, 12 Mar 2020 12:13:44 +0800

最近效率有点低，需要改进一下时间管理系统。

加一个限制，每天在推荐系统上沉迷的时间不能超过1h。

先从2020.3.12开始坚持两周吧。

M1 Finance

Tue, 10 Mar 2020 12:25:42 +0800

昨天美股触发了熔断，不过我重仓的是TMF 3x美国长期国债，涨了不少，赶紧撤出了。

撤出了之后，发现波动这么大的时候感觉好累，打算疫情结束之前不碰美股了。

Chrome中使用新版QQ邮箱的提醒插件

Tue, 10 Mar 2020 12:25:30 +0800

首先是打开微信登录的QQ邮箱。

“QQ邮箱-设置-邮件通知-通知插件下载”有一个按钮下载，可以尝试按一下，会提示包无效“CRX_REQUIRED_PROOF_MISSING”，这是因为Chrome/Edge出于安全考虑禁止了.crx的下载、安装。

数字娱乐设计二学位

Sat, 07 Mar 2020 20:00:13 +0800

“数字媒体艺术专业(数字娱乐设计方向)”是清华大学美术学院面向非美院学生开设的第二学士学位，也简称美双。这个简称在2018年之前是没有歧义的，因为这是美术学院开设的唯一一个二学位项目，班号也类似美双16、美双15。在2018年美术学院开设了另一个二学位项目，于是班号也就变成了美双71、美双72 / 美双81、美双82。

童话镇

Sun, 01 Mar 2020 20:05:59 +0800

Youtube给推荐了陈一发儿的《童话镇》，听了听，还是很好听。

只有睿智的河水知道
睡美人逃避了生活的煎熬
小人鱼把阳光抹成眼影
投入泡沫的怀抱

想起来，这首歌根据网易云的统计应该是2017年有200+天都听过这首歌，所以现在听到有一种回到过去的感觉。

2.29

Sat, 29 Feb 2020 20:03:35 +0800

今天是4年一遇的2.29。

上一次是2016年，应该还在保送后的长假里。

祝世界越来越美好吧。

不同种类冠状病毒生成的进化树

Thu, 27 Feb 2020 17:43:28 +0800

计算生物学课程的第一次作业。

用不同种类冠状病毒生成的进化树

红色的是2019-nCov，最接近的那个是个穿山甲体内的病毒。

进军美股

Sun, 23 Feb 2020 10:25:18 +0800

准备去美股送点钱，体验一下资本主义国家资本市场的险恶。

在和雪球比较紧密的“雪盈证券”开的美股帐户，选的低佣金方案。

找同学收了点美元，加上之前境外没花完的美元，准备一共撒$3900。

在线上课的一周

Sun, 23 Feb 2020 00:01:40 +0800

这周是开学的第一周，也是线上上课的第一周。记录我的感受。

计算生物学

老板的课，使用了Zoom会议软件，只有老师开摄像头、语音、屏幕共享，授课状况良好。

雨课堂出现了一些故障，没有使用。

Namebase Airdrop

Fri, 21 Feb 2020 23:18:24 +0800

因为这个让我薅到了真金白银，还是记录一下。

那天tuna的群pin了条消息，点进去看了眼那个官网，结果没看懂，就忽略了…

昨天又看到某中文领钱教程才下手的。

🔒 Emoji AES

Tue, 11 Feb 2020 10:31:35 +0800

GitHub上发现了一个库Emoji-AES，把AES加密后的串Base64，然后mapping到了emoji。

Message

OI(信息学竞赛)课件

Mon, 10 Feb 2020 10:15:22 +0800

最近有一阵风，大家纷纷开始公开自己做过的课件。

这里记录一下方便日后引用。

自己的是没精力去整理了，毕竟每次都是讲课前一晚上拼拼凑凑。

hzwer: https://github.com/hzwer/shareOI
Trinkle23897: https://github.com/Trinkle23897/oi_slides
huzecong: https://github.com/huzecong/oi-slides
haoxizhong: https://github.com/haoxizhong/problem

几个远程会议软件

Sun, 09 Feb 2020 15:34:06 +0800

作为下学期某门课程的助教，响应学校远程教学的政策，帮老师测试了一下这些远程会议软件。

雨课堂

ppt没有动画
交流不方便，学生没法说话
因为在微信里，所以和微信群聊冲突
课程结束下发ppt时，没有播放的ppt也会被下发

华为WeLink

不支持mac
拉人困难，需要一个一个添加手机号
其他音视频效果良好

Zoom

音视频效果良好
支持共享屏幕时标注
支持白板
软件内聊天、举手不会有明显提示
全平台支持(除linux)

腾讯会议

音视频效果良好
没有标注功能
没有白板功能
注册简便

Skype

有些卡顿
官方宣传的多方视频通话建议5个人

最终可能使用Zoom配套雨课堂，腾讯会议备选。

现金流与杠杆：我的信用卡们

Sun, 09 Feb 2020 12:44:14 +0800

第一张信用卡

大学刚入学的时候，需要给自行车喷号码。在登记自行车的摊位旁边就是建行信用卡摆摊的摊位。业务员说可以办信用卡，但是想要额度的话需要家长填一个还款保证协议。但是我的家长不在，也不想麻烦他们，就直接提交了申请。

本周主题：防疫、英语和科研

Tue, 21 Jan 2020 19:49:00 +0000

昨天晚上把本学期的最后一份作业交了，希望老师不要挂我。

今天和久违的switch亲密了一下，终于把四神兽都打死了(子弹时间+炸弹箭)，然后去玩了会dlc的必杀之剑任务。

2019年终总结

Tue, 31 Dec 2019 08:00:00 +0800

又到年底了，写一下年终总结。就只讲开心的事情吧。

1月

开了一个telegram频道@chanstring，已经破产。

大三上的期末季。

寒假和班里同学一起去新加坡~~参加TCS的一个冬令营~~玩。

2月

春节过年。

假期在家学了波怎么做菜，最后也没学会几个。

2019年的最后一周

Mon, 23 Dec 2019 21:20:42 +0800

这周大概是今年的最后一周，要做的事情有很多。有一些是最近新产生的锅，有一些是之前该做的事情拖到了现在。

Date	Weekday
12.24	周二	网络原理最后一次书面作业
12.24	周二	选好二学位毕设的导师
12.26	周四	算法分析与设计最后一次作业
12.26	周四	网络原理路由器组网验收
12.27	周五	二学位毕设预开题
12.31	周二	小毕设终期答辩
1.2	周四	网络科学考试
1.4	周六	网络原理考试
1.5	周日	小毕设结题报告提交

网络科学project
算法分析与设计论文阅读报告

努力把这些事情一个一个一个一个一个做完。

[视频] 早起去上课

Sat, 14 Dec 2019 21:20:25 +0800

Youtube 嵌入：

youtu.be/P8ogEiHo_jY

Bilibili 嵌入：

av79142325

2019年冬天的第一场雪

Sat, 30 Nov 2019 13:57:00 +0000

今天是2019.11.29，北京下了今年冬天的第一场雪。

这场雪比前几年来的更早一些。

出去走了走，心情不错。

Macbook Pro 16寸

Sat, 30 Nov 2019 13:55:00 +0000

现在用的电脑是Macbook Pro 13′ Early 2015，已经4年多了，决定换一下。

11.13 晚在微博上看到了新Macbook Pro在官网上架，实体Esc，键盘换回去了，16’屏幕，感觉还不错。

成都

Fri, 20 Sep 2019 01:02:45 +0800

和我在成都的街头走一走，直到所有的灯都熄灭了也不停留。

这首歌还是蛮好听的。之前因为参加竞赛的原因去过两次绵阳，但一直没有去过成都。于是在回国之后、开学之前，决定去成都看一看。

9.15

Sun, 15 Sep 2019 21:19:13 +0800

上午完成了保研面试。

中午去20公里外吃了一顿饭。

晚上开了班会，听各位大佬讲暑研成果。

New York City

Thu, 29 Aug 2019 01:03:13 +0800

前往纽约&新泽西旅游。

枯燥的一天

Thu, 22 Aug 2019 21:19:00 +0800

早上8:20，手机闹铃响了， Siri聪明的建议“要不要启用9:25的闹钟”。

于是我接受了Siri的建议。

早上9:25，闹钟又响了，匆匆起来，倒了一杯 1% fat的牛奶，拿出了几块点心作为早餐。

用信息买醉

Thu, 15 Aug 2019 14:48:00 +0000

当对某件事、某个话题、某个东西突然特别感兴趣时，
在Google翻几十页关于这个关键词的搜索结果，
在知乎、微博、豆瓣等平台上寻找关于这个话题的一切信息，
津津有味的阅读，越看越上瘾，时间默默流走。

浮生若梦 Chicago

Tue, 16 Jul 2019 07:20:15 +0000

这次来芝加哥已经有4个半月了，预期还要再待一个半月。
在一个地方一个人待久了，就会进入另一种状态。
进入一种旁观者的状态，感觉每天自己做的事情和这个世界很独立，世界发生什么好像都与自己无关。

“拥有对时间控制的绝对自由”，这绝对是一件有挑战的事情。之前的浮躁和拖延将会被无限放大，可能只有在产生了某些不好的后果之后，才能痛定思痛做出改变吧。

Sublime Text 弹窗 Patch

Mon, 08 Jul 2019 21:13:13 +0800

Sublime 是一款优秀的文本编辑器。

什么，你是 vs code 用户? 那我只好推荐你用这款插件了。

正版售价$80，但是可以无限试用，不过在一定时间后会弹窗提示。

弹窗

下面记录一下使用 Hopper 软件将该弹窗去掉的过程。

首先用 Hopper 打开 /Applications/Sublime\ Text.app/Contents/MacOS/Sublime\ Text。

用搜索功能找到提示的文字，并找到对该地址的引用。

又是一年高考季

Fri, 07 Jun 2019 21:18:18 +0800

murmur	date
食堂的自助一顿要18usd，不过真的好好吃…	2019/6/7
自动驾驶又左转撞直行了..	2018/6/7
这微积分真TM难，这都啥，抄抄抄。。	2017/6/7
“为了艾泽拉斯！《魔兽》首映还不错。”	2016/6/7
“越努力、越幸运！”	2015/6/7
“希望不用高考啊&怎么每年高考都下雨。”	2014/6/7

Margin

Sun, 26 May 2019 14:47:00 +0000

记着当时看完《阿丽塔》有一些关于“在规则内反抗规则制定者”的感叹，可能是旁观视角才让我得到这种感受吧，现实的规则制定者是不会让你有这种感觉的。

飞聊flipchat

Mon, 20 May 2019 21:17:45 +0800

头条旗下的飞聊上线了，跟之前推出的"多闪"比基本去掉了所有短视频的元素，但整个应用框架、UI相似。应该比之前短暂火过的子弹短信能打得多。

收藏夹分享-2

Tue, 30 Apr 2019 21:17:31 +0800

一些有趣的链接，顺序随机。

http://jxgl.cic.tsinghua.edu.cn/jxpg/f/xssq/exwfx/xssqb/cksqs 通过二学位的接口查询 GPA。(先登录info)

https://www.behance.net/ Behance，一个创意设计社区。

https://flarum.org/ Flarum，一个php写的轻论坛。

收藏夹分享-1

Sat, 27 Apr 2019 21:17:25 +0800

一些有趣的链接，顺序随机。

https://github.com/hzwer/LearningToPaint A painting AI that can reproduce paintings stroke by stroke using deep reinforcement learning.

http://bin.mailgun.net Mailgun 提供的postbin服务，会生成一个随机的url，向这个url POST信息会记录下来。

《告白の夜》小提琴谱

Fri, 26 Apr 2019 21:17:01 +0800

很好听，收藏一下谱子。

来源：http://www.qupu123.com/puyou/shangchuan/p340100.html

upd

https://www.youtube.com/watch?v=qDvKVxFe7VA

发现了这个视频，把位和揉弦都能看得很清楚，适合学习。

不会受伤的机器人

Sat, 20 Apr 2019 21:16:09 +0800

当你和我在一起的时候，无论你对我做什么，我都感到开心。我不在身边的时候，我知道你去做对你而言更重要的事了，我一样也感到开心。

我不知道我的开心是发自内心的，还是因为我只有这一种感情。但这也无所谓，毕竟，我只是一台机器而已。

Fancybox如何直接兼容img标签

Wed, 17 Apr 2019 21:15:48 +0800

TL;DR

把下面这段代码加到html里就可以直接点开<img>标签了。

function createFancyboxEvent() {
    $('img').each(function(){
        this.onclick = function(){
            $.fancybox.open({
                src: this.src,
                type: 'image',
            })
        }
    })
}
$(function(){
    createFancyboxEvent();
})

详细

点击上面这个图片是可以放大的。

这个功能是靠 Fancybox 来实现的，我看 hexo 插件里有这个，但是我设置了fancybox=true后并不能点开图片。

低垂的果实

Sat, 13 Apr 2019 12:50:00 +0000

摘完了低垂的果实后，不要变得懒惰。

说的还是简单，做到就太难了。

有一个老故事，叫“上帝已经给了你三次机会”。

现实中，上帝可能比这宽容的多，但如果上帝已经给了你十次机会了，而你还在测试上帝的底线，是不是就有点作死了。
本文不代表任何宗教观点

窗边日落

Thu, 04 Apr 2019 21:14:27 +0800

日落

还记着来看房时的那次日落让我下决心住在这里。

Numba是个好东西

Tue, 02 Apr 2019 21:14:22 +0800

logo里的这个python快如闪电⚡️。

要用python做计算密集型任务的可以考虑一下这个库，非常香。

http://numba.pydata.org/

饭否关键词分析

Mon, 18 Mar 2019 17:52:13 +0800

饭否是中国大陆地区首批提供微型博客服务的类 Twitter 网站之一，用户可以通过网页、WAP 页面、手机短信/彩信（现已不可用）、IM 软件（包括 MSN、Google Talk）和上百种 API 应用发布消息（限 140 字以内）或上传图片。用户间通过互相关注、私信、或@对话等方式互动。(wikipedia)

其 API 完善，搜索消息接口可以完整遍历含有某个关键词的所有消息。

为了好玩，写了个 ipynb 来对某个关键词进行自动分析。

芝加哥初见闻

Thu, 07 Mar 2019 21:14:08 +0800

到芝加哥大概一周了，验证了地球是圆的这个事实。

来的时候做的是美联航的飞机。到的前两天吃饭很狼狈，后来去了趟chinatown，买了点材料，稍微好了一些。

前往芝加哥

Sat, 02 Mar 2019 01:03:20 +0800

去TTIC research intern 半年，生物信息学方向。

重新开始

Mon, 28 Jan 2019 21:13:50 +0800

重新安装了macos，感觉速度快了很多。

应该是硬盘清空了，之前硬盘可用空间一直<5G，难为ssd主控了。

在新加坡的时候，趁着有十几MB的上传速度，把所有有用的数据上传到了dropbox，大概有15GB左右，新的数据也有美区icloud同步，不自己弄owncloud了（上次差点把数据弄丢）…

前往新加坡

Tue, 22 Jan 2019 01:03:25 +0800

来到了新加坡，来参加Tsinghua-NTU Theoretical Computer Science Winter School 2019。

这好像是第二次来新加坡，上一次好像是初二的时候去参观法嘉中学。

MD5碰撞

Thu, 17 Jan 2019 21:13:30 +0800

这两张图的md5相同，都是02d465adff54099619f704f7b1960b42。

ok.jpg

gg.jpg

Reference: https://github.com/cr-marcstevens/hashclash

2018年终总结

Mon, 31 Dec 2018 08:00:00 +0800

总觉得2018跨到2019其实没什么变化，但是周围的环境又会让你在这个时刻有些反思，确定一些新的目标。

下面是一些随机的，现在能想到的一些感想。

感觉2018有很长一段时间是在焦虑中度过的，手头的事情很多，有些事无从下手，有些事无能为力，但是还是需要一直想着..现在倒是没那么焦虑了，列个ddl的list一件一件去做就好了..

Codeforces可以改handle了

Sat, 29 Dec 2018 21:11:37 +0800

Codeforces可以改handle了。

今年还可以改成已有的非活跃的账号的名字。

然后随便试了几个词，发现hello这个号last visit 8 yrs ago，然后就改成了hello了 233333。

Zeitgeist 2018

Thu, 27 Dec 2018 19:43:00 +0000

2019加油。

R君的弯道

Thu, 23 Aug 2018 21:11:01 +0800

上个周末，见到了R君，R君刚过完20岁生日，相谈甚欢，聊着聊着就聊到了R君的人生经历。

R君说，他的经历，就像两次“弯道超车”，于是慢慢听R君讲起了人生路上的每一个节点。

20岁生日

Mon, 13 Aug 2018 21:10:12 +0800

今天是20岁的生日。

10:00起床。

步行十分钟去上班。

改之前代码的bug。

午饭在万利达二楼食堂。

写新的bug。

晚饭在万利达三楼食堂。

21:00下班。

买了回去上小学期的机票。

其实都一样

Sat, 28 Jul 2018 21:09:56 +0800

某兴恢复了经营，有了这次的政府背书，以后股价会更高吧…

其实这个明面上的事件，和一些阴暗面的事件是一样的。

但是阴暗面又不能说。

经历过越大的风波，之后就越有可吹嘘的，接下来的泡沫就能吹的更大。

腾讯实习两周感悟

Tue, 17 Jul 2018 21:09:49 +0800

原先觉得很娱乐、很随意的东西，原来有这么一群人在认真的讨论、决策。

游戏直播头部真赚钱。

微信小程序野心很大。

实习生项目得了个奖。

去讲课的课件还没做。

外面开始下暴雨了。

腾讯实习一周感悟

Tue, 10 Jul 2018 21:09:43 +0800

🙊

实习生入职，这一批工卡的带子和其他的不一样。

毕竟大厂大厂，内部的网络等管理还是都有一套的，服务质量也很高很高。

漫步在深圳街头，慢慢就不觉得热了。

然后导师交代的第一个任务做了三天才完成。

射击课

Sat, 14 Apr 2018 00:56:34 +0800

这学期选了射击课，还蛮有趣的。上半学期是气手枪，下半学期是气步枪。在这里记录一下成绩。

Week1

2017年终总结

Sun, 31 Dec 2017 08:00:00 +0800

再见2017

2017结束了。

2017是完整在大学生活的第一年。

一切都在改变，希望能不忘初心。

希望在2018，找到努力的方向，多做自己感兴趣的事，多认识有趣的人。

屠龙故事

Wed, 22 Nov 2017 21:08:05 +0800

1.

有一条恶龙，每年要求村庄献祭一个处女，每年这个村庄都会有一个少年英雄去与恶龙搏斗，但无人生还。又一个英雄出发时，有人悄悄尾随。龙穴铺满金银财宝，英雄用剑刺死恶龙，然后坐在尸身上，看着闪烁的珠宝，慢慢地长出鳞片、尾巴和触角，最终变成恶龙。

小Tip: 使用其他设备连接Tsinghua WiFi

Fri, 13 Oct 2017 21:18:06 +0800

有些时候，需要使用非智能设备去连接Tsinghua WiFi。

可能没法打开net.tsinghua.edu.cn/ usereg.tsinghua.edu.cn去认证，没法获取到设备的IP地址，甚至一些设备只有在通过这个WiFi能够连接网络的时候才会使用这个WiFi(比如Switch)。

消除网络学堂红色的0

Wed, 20 Sep 2017 21:07:43 +0800

消除网络学堂红色的0。

示例效果：

网络学堂

Chrome下载插件tampermonkey，这是一个可以在匹配页面跑js脚本的插件。

然后新建插件，写下面内容，启用插件：

OAuth2试玩

Tue, 29 Aug 2017 21:07:30 +0800

温馨提示：这篇文章已经过时

把Isso评论魔改成了github登录发表only 233333。

用requests和flask手撸了个OAuth2验证。

具体魔改的内容可以看https://github.com/zrt/isso-OAuth2 2333。

加载速度过度优化

Thu, 24 Aug 2017 21:07:08 +0800

** 温馨提示：这篇文章已经过时 **

这篇文章讨论如何过度优化静态博客的加载速度。

闲的无聊，于是思考，怎样才能让zrt.io/blog加载的更快呢？

这个blog是由hexo生成，所以可以看作纯静态的一堆html。

修复Apache2 mpm prefork与HTTP/2不兼容问题

Wed, 23 Aug 2017 21:06:53 +0800

配置好了http2，通过curl -v –http2发现并没有走h2，在/var/log/大概报个这样的warning：

[http2:warn] AH10034: The mpm module (prefork.c) is not supported by mod_http2. The mpm determines how things are processed in your server. HTTP/2 has more demands in this regard and the currently selected mpm will just not do. This is an advisory warning. Your server will continue to work, but the HTTP/2 protocol will be inactive.

大概是mpm的prefork模式和http2不兼容。Google到修补方法：修改mpm模式为event。参考https://http2.pro/doc/Apache#prefork-http2

六教的包

Wed, 24 May 2017 21:06:40 +0800

wireshark挂着，抓抓空气中的数据。

写了个脚本利用抓到的cookie。

大概长这样：

def cookie_to_js(s):
    res = ''
    a=s.split(';')
    for x in a:
        if '=' in x:
            res+='document.cookie="%s";'%(x)
            
    return res

print(cookie_to_js('a=b;c=d;'))

有其他好用的工具的话就见笑了，求推荐。

http://apilive.zhihu.com/apilive 这个websocket好像会被客户端访问，没有用https，可以抓到cookie。

Memobird咕咕机Python开发

Sat, 04 Mar 2017 21:06:33 +0800

Memobird咕咕机是一个热敏纸打印机。官网。

官方提供了开发的web API接口: 开放平台，API文档。

通过web API可以做许多奇怪的事了。

不过官方并没有提供示例程序，于是我按照接口说明写了一下Python3的代码实现。

小米路由器mini安装 OpenVPN client

Mon, 13 Feb 2017 07:51:00 +0000

温馨提示：此篇内容可能过时，请谨慎参考

因为刚到的memobird需要用无需认证的wifi环境才能用，然而发现寝室并没有这样的wifi，于是打算改造一下~~已经废了的~~201A_5G。

大概思路是先把路由器刷成PandoraBox，然后用opkg安装openvpn，然后接入vpn9，此处感谢贵系科协，最后路由器用iptables转发一下连接路由器的流量。

注意我的路由器版本是小米路由器mini，其他版本应该类似。

2016年终总结

Sat, 31 Dec 2016 08:00:00 +0800

再见2016

2016结束了。

2016: 中学毕业，来到大学。

2016: 前3/4年喊着要学微积分，后1/4年被微积分虐傻。

还不知道未来的方向，还是只单身汪。

希望在2017，多做自己感兴趣的事，多认识有趣的人。

启用了BBR

Tue, 20 Dec 2016 18:02:56 +0800

温馨提醒：这篇文章已经过时，请不要参考

Ubuntu如何启用BBR

关于BBR –知乎

在把DigitalOcean SF1的Ubuntu16.04启用BBR后网速x20，油管4k无压力，效果非常好，记录一下步骤。

Keybase.io 介绍

Mon, 21 Nov 2016 18:02:45 +0800

Keybase是什么？

Keybase.io是一个用于交换公钥的社交网站，用于证明你是你，以及加密解密信息。

Keybase.io提供了一套PGP加密、解密、签名、验证签名的工具，还提供了其他社交网络账号、github、比特币地址、域名所有权的验证服务。

word.ruotian.io

Wed, 21 Sep 2016 18:02:33 +0800

温馨提示：这篇文章已经过时

决定背新东方的TOEFL绿宝书，一周一个LIST，一年背完。

随手写了个网页随机出现这周要背的词，方便背和复习。

欢迎其他同学按我的进度一块背。

毕业蹭饭图制作

Tue, 20 Sep 2016 18:02:16 +0800

图多预警！

昨天，做了个图：

蹭饭图

…很有ingress的风格，发出来后收到许多同学的好评。

下面说说是怎么做出来的。

首先想想要做成什么样的呢？

用Python3写了个OJ

Mon, 19 Sep 2016 18:02:03 +0800

这个OJ以后可能搭在学校吧..(论造轮子的修养)

https://github.com/zrt/XOJ

主要为了练习下python和HTML。

分为web和judger两个部分。

web主要处理OJ非评测的部分，judger主要处理评测部分。

手写公式，生成LaTex代码!

Sun, 18 Sep 2016 07:44:00 +0000

经常需要写LaTeX数学公式。
但是并不是LaTeX语言大师。
碰到不知道怎么写的情况怎么办呢？

有一个好用的网站(Chrome打开): http://webdemo.myscript.com/views/math/index.html，可以手写公式,告诉你对应的LaTeX语句。

Hack ghost: http://detexify.kirelabs.org/classify.html

我的初中高中

Fri, 10 Jun 2016 18:01:22 +0800

记录留在脑海里的一切。

阶段1 初三开始-最后一次分班前

初中在石家庄41中。

初中所在的班级不是特别厉害的那种班，但至少还是学习为主流的那种，有个非常负责的班主任。
感觉压力不大的度过了初一初二的生活。我当时在班里的成绩top5左右，年级成绩大概100-200名（一共年级1300人左右）。好像是数学课代表（曾经有过几次没做作业上去讲题的经历233），然后语文非常渣。
也有几个比较要好的同学。感觉一切都很正常，到了初三也很努力，按部就班的在学习，成绩也一点一点的进步。
后来什么中考体育、实验、微机也都轻轻松松。
在记事本中还发现一些莫名其妙的话。

我的集训队作业

Fri, 13 Nov 2015 17:18:01 +0800

第一次作业

我的题目分别是Codechef Nov 2014的三道题：

Chef and Churu https://www.codechef.com/problems/FNCS
Sereja and Order https://www.codechef.com/problems/SEAORD
The Spelling Problem https://www.codechef.com/problems/SPELL

福利内容？

自认为写的题解还不如官方题解。QAQ。
T1简要题解: 分块+树状数组。对函数分块，预处理每个位置在每一块里出现次数。
T2简要题解: 贪心+排序。证明答案一定等于max(max(ai+bi),max(sumA,sumB))，然后贪心或随机化。卡掉了几个不正确的贪心。
T3简要题解: hash或trie树。给字典中的单词赋一个权值。对于每个不在字典中的词，尝试所有的出错方法，取权值最高的那个词。

旅游美景

Fri, 07 Aug 2015 17:17:46 +0800

NOI考完了心情不错。
既然来了南方，就好好转一转吧。

路线：杭州->上海->苏州->南京。

[图片挂了]

NOI2015 行记

Tue, 21 Jul 2015 17:17:20 +0800

NOI2015这么快就到了呢。

作为弱省高一开始学OI的蒟蒻。

一辈子也就参加这一次NOI吧。

写写日记记录下吧。

感觉时间过的好快。

好像昨天才听说了OI，今天却已经到了国家赛的现场。

THUSC2015 行记

Tue, 02 Jun 2015 21:26:46 +0800

THUSC: Tsinghua University Summer Camp.

我参加的当然是OI营咯。

之前一直纠结是去P大夏令营还是T大夏令营。

经过与学长的交流还有自己的思考终于决定去THU的啦。

sqy去了P大的夏令营。

HEOI2015 行记

Sun, 26 Apr 2015 17:17:05 +0800

陪你开始OI的人未必能陪你走到散场。一次擦肩而过之后可能就再也见不到了。

这几天过得好乱，也不知该说些什么。

Day 0

今年省选在河北师大，和去年noip一样。

题解 CodeChef March Lunchtime 2015 No Unpaired Chefs

Sun, 29 Mar 2015 17:16:55 +0800

中文题面
 题目链接
题意:
给你一个矩阵$M$,满足$M_{i,j}=M_{j,i}$，且$M_{i,j}>=0$。
让你构造一个矩阵$P$,满足：

$P$是对称矩阵。
$P_{i,j}$要么等于零，要么等于$M_{i,j}$。
存在一个$S>0$,使得$P+P^2+…+P^S$中每个元素都是正数。
最小化$P$中元素和。

$M$是一个邻接矩阵，可以看成一个完全图。我们知道$(M^k)_{i,j}$对应所有走$k$条边从$i$到$j$边权乘积的和。
满足第3个条件就需要所有点两两可达。
于是就是最小生成树了。

题解 BZOJ 1089 [SCOI2003]严格n元树

Tue, 24 Mar 2015 17:16:37 +0800

设$f[i]$表示深度不超过i的严格n元树的数量。
有$f[i]=f[i-1]^n+1$。
$Ans=f[d]-f[d-1]$
递推式很显然。

代码:

n,d=map(int,raw_input().split())
f=[0]*17
f[0]=1
for i in range(1,d+1):
	f[i]=f[i-1]**n+1
print f[d]-f[d-1]

这个题目答案可能很大很大.. 所以最坏情况下是跑不出来的…

题解 BZOJ 3586 字符串生成器

Sat, 21 Mar 2015 17:16:27 +0800

题意:

　　有一个字符串生成器，初始时生成的字符串为空串，它每次按照给定概率随机生成一个小写字母，加在当前已生成字符串的后面。
　　给定N个长度为L的字符串，每个字符串由小写字母组成。
　　如果在某个时候，发现每个给定字符串都在当前已生成的字符串中作为子串出现过，生成器就会停下来，将当前生成的字符串作为输出。
　　求输出字符串长度的期望值。

题解 BZOJ 3489 A simple rmq problem

Fri, 20 Mar 2015 17:15:12 +0800

给出M个询问：在[l,r]之间找到一个在这个区间里只出现过一次的数强制在线。

设$p[i]$为$a[i]$上一次出现位置,$q[i]$为$a[i]$下一次出现位置。

每次查询的就是
$ max _{l<=i<=r} a[i] $ 且$q[i] >r$ , $p[i]< l $。

可以排序然后树套树做。

偷懒写了KD-TREE。。复杂度$O(n^{5/3})$，不强制在线的话可以做到$O(n^{3/2})$ 。

题解 BZOJ 2280 [Poi2011]Plot

Fri, 20 Mar 2015 17:15:07 +0800

给出一系列点p_1, p_2, … , p_n，将其分成不多余m个连续的段，第i段内求一个点q_i，使得q_i到这段内点的距离的最大值的最大值最小

最大值最小，先二分答案。
再从每个开始点二分，用随机增量法的最小圆覆盖判断最长能扩展长度，段数是不是$<=m$。
但是有个问题。
如果m=n的话，每段长度都是1。
这样做是$O(n^2lgn)$的。

题解 BZOJ 3611 [Heoi2014]大工程

Wed, 18 Mar 2015 17:14:53 +0800

去年省选的题。
突然发现还没有A这个题。
就抽时间写了写。

询问的那个树形dp很好搞。
我们把原树中有用的点(被选的点，和他们之间的lca)拿出来，建一棵虚树。
在虚树上dp就行了。
类似"消耗战"(为啥听着像卡组名)。

题解 BZOJ 2212 [Poi2011]Tree Rotations

Wed, 18 Mar 2015 17:14:43 +0800

学习 http://www.docin.com/p-598829500.html 线段树合并中的一道例题。
详细见那个课件。
复杂度$O(nlgn)$.
需要注意下空间

code:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
int n;
int LS[400005],RS[400005];
int cc;
int Root;
int w[400005];
int siz[400005];
int ls[15000500],rs[15000500],sum[15000500];
int num;
int stk[15000500],top;
int getnew(){
	if(top>0){
		ls[stk[top-1]]=rs[stk[top-1]]=sum[stk[top-1]]=0;
		return stk[--top];
	}
	return ++num;
}
void add(int &o,int l,int r,int pos){
	if(!o) o=getnew();
	if(l==r){
		sum[o]=1;
	}else{
		int mid=(l+r)>>1;
		if(pos<=mid) add(ls[o],l,mid,pos);else add(rs[o],mid+1,r,pos);
		sum[o]=1;
	}
}
int read(){
	int x;
	scanf("%d",&x);
	int now=++cc;
	if(!x){
		LS[now]=read();
		RS[now]=read();
		siz[now]=siz[LS[now]]+siz[RS[now]]+1;
	}else{
		add(w[now],1,n,x);
		siz[now]=1;
	}
	return now;
}
LL ans;
LL ans0,ans1;
int merge(int a,int b,int l,int r){
	if(!a) return b;
	if(!b) return a;
	int now=getnew();
	int mid=(l+r)>>1;
	int lft=sum[ls[a]]+sum[ls[b]];
	int rgt=sum[rs[a]]+sum[rs[b]];
	ans0+=sum[rs[a]]*1LL*sum[ls[b]];
	ans1+=sum[ls[a]]*1LL*sum[rs[b]];
	if(lft<rgt){
		rs[now]=merge(rs[a],rs[b],mid+1,r);
		ls[now]=merge(ls[a],ls[b],l,mid);
	}else{
		ls[now]=merge(ls[a],ls[b],l,mid);
		rs[now]=merge(rs[a],rs[b],mid+1,r);
	}
	stk[top++]=a;
	stk[top++]=b;
	sum[now]=sum[ls[now]]+sum[rs[now]];
	return now;
}
void dfs(int x){
	if(!LS[x]) return;
	else{
		if(siz[LS[x]]<siz[RS[x]]) dfs(RS[x]),dfs(LS[x]);
		else dfs(LS[x]),dfs(RS[x]);
		ans0=ans1=0;
		w[x]=merge(w[LS[x]],w[RS[x]],1,n);
		ans+=min(ans0,ans1);
	}
}
int main(){
	http://www.docin.com/p-598829500.html
	scanf("%d",&n);
	Root=read();
	dfs(1);
	printf("%lld\n",ans);
	return 0;
}

题解 BZOJ 2758 Blinker的噩梦

Tue, 17 Mar 2015 17:14:35 +0800

题意:

给你n个不相交的凸多边形或圆，每个图形有一个权值。
询问m次，从一个点到另一个点，穿越过的图形的权值异或和。

题解:

扫描线+树链剖分。
这些图形的互相包含关系形成了一些树。
每次相当于求从树上一个点走到另一个点边上的权值异或和。
扫描线算法$O(nlgn)$求出这颗树，然后树链剖分维护下就好了。

bzoj-problem-chooser

Thu, 12 Mar 2015 17:14:24 +0800

今天上午无聊写了这个。

BZOJ抽题助手。

估价函数综合 AC率，AC人数，提交人数，随机因素。

还能指定难度。

记录已A的题，跳过设置的黑名单。

计算出这道题与你输入的难度的差距难度系数。

题解 BZOJ 3460 Jc的宿舍

Thu, 12 Mar 2015 17:14:13 +0800

贪心一下。
接水的时候肯定是从小到大接。
好像没什么性质。莫队好了。
观察到奇怪的强制在线方法。。。
类似[糖果公园/苹果树]做法。
块大小没仔细算，直接$n^{0.5}$的，卡了70s评测姬 233。
复杂度不明。。

题解 BZOJ 3227 [Sdoi2008]红黑树(tree)

Thu, 12 Mar 2015 17:14:04 +0800

求红黑树红色节点的最大最小值。
参考百度百科

性质1. 节点是红色或黑色。性质2. 根节点是黑色。性质3 每个叶节点（NIL节点，空节点）是黑色的。性质4 每个红色节点的两个子节点都是黑色。(从每个叶子到根的所有路径上不能有两个连续的红色节点) 性质5. 从任一节点到其每个叶子的所有路径都包含相同数目的黑色节点。

题解 BZOJ 2565 最长双回文串

Thu, 12 Mar 2015 17:13:58 +0800

求最长双回文串.
首先Manacher求回文半径p[i]。
然后求每个点最左是可以由哪一个点扩展来的left[i]。
具体就是拿个数组弄一个 i+p[i] 的后缀min。
刚开始想单调队列好像也行。
然后对于一个点以它为右回文串中心的最大双回文串长度就是i-left[i-p[i]+1].
取个max就行了。

题解 BZOJ 1970 [Ahoi2005]Code 矿藏编码

Thu, 12 Mar 2015 17:13:50 +0800

写个高精度。
然后扫过去就行了。
抽题器竟然抽了这么水的题233.

code:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int k;
struct N{
	int a[100];
	friend N operator + (const N&a,const N&b){
		N c;
		for(int i=0;i<100;i++) c.a[i]=a.a[i]+b.a[i];
		for(int i=0;i<99;i++) c.a[i+1]+=c.a[i]/10,c.a[i]%=10;
		return c;
	}
	N mul2(){
		N b;
		for(int i=0;i<100;i++) b.a[i]=a[i]*2;
		for(int i=0;i<99;i++) b.a[i+1]+=b.a[i]/10,b.a[i]%=10;
		return b;
	}
	void out(){
		int top=99;
		while(top>0&&a[top]==0){
			top--;
		}
		while(top>=0) printf("%d",a[top--]);
		puts("");
	}
}po2[105],ans;
char s[205];
int p;
void solve(int k){
	if(s[p]=='0') ans=ans+po2[2*k],p++;
	else if(s[p]=='1'){
		p++;
	}else{
		p++;
		solve(k-1);
		solve(k-1);
		solve(k-1);
		solve(k-1);
	}
}
int main(){
	po2[0].a[0]=1;
	for(int i=1;i<=100;i++){
		po2[i]=po2[i-1].mul2();
	}
	scanf("%d",&k);
	scanf("%s",s+1);
	p=1;
	solve(k);
	ans.out();
	return 0;
}

题解 BZOJ 1201 [HNOI2005]数三角形

Thu, 12 Mar 2015 17:13:39 +0800

首先看正着的三角形。
每个点向右上方可以扩展的距离叫 mxl[],左上方叫mxr[]。
枚举三角形右下端点。
设右下为j.坐下为i。
需要满足底边i..j 联通，并且$j-i<=mxr[j]$且$j-i<=mxl[i]$。
把i挪到一边，把j挪到一边就是
$j-mxr[j]<=i$且$j<=mxl[i]+i$

题解 BZOJ 1845 [Cqoi2005] 三角形面积并

Wed, 11 Mar 2015 17:13:29 +0800

题目是求三角形面积并。$n<= 100$
就像hwd冬令营上讲的。
先$n^2$求交点。
然后在每两个交点之间求梯形中位线总长，然后求和。
ps.听说simpson也可以。。
…被卡精度了…

题解 BZOJ 1023 [SHOI2008]cactus仙人掌图

Wed, 11 Mar 2015 17:13:18 +0800

题意是仙人掌图直径。
大概做法就是树的直径dp求法加上基环树的技巧特技。
改了好长好长时间啊。
参考题解: http://z55250825.blog.163.com/blog/static/150230809201412793151890/
很细致。。

code:


#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m;
int H[50005],X[20000005],P[20000005],tot=1;
int dp[50005];
inline void add(int x,int y){
	P[++tot]=y;X[tot]=H[x];H[x]=tot;
}
int stk[50005],top;
int dfn[50005];
int low[50005],tim;
int ans;
int list[100005],t;
struct N{
	int pos,w;
	N(int a=0,int b=0){
		pos=a;w=b;
	}
};
struct dddl{
	N q[100005];
	int h,t;
	void clear(){
		h=t=0;
	}
	void pop(int x){
		if(q[h].pos==x) h++;
	}
	int ask(){
		return q[h].w;
	}
	void push(N a){
		while(h<t&&q[t-1].w<=a.w) t--;
		q[t++]=a;
	}
}q;
bool debug;
int fa[50005];
void dfs(int x){
	dfn[x]=low[x]=++tim;
	for(int i=H[x];i;i=X[i]){
		if(P[i]==fa[x])continue;
		if(fa[P[i]]&&dfn[P[i]]>dfn[x]){
			//find circle
			int k=P[i];
			t=0;
			while(k!=x) list[t++]=k,k=fa[k];
			list[t++]=x;
			int mx=0x80808080;
			for(int j=0;j<t;j++){
				list[j+t]=list[j];
				mx=max(mx,dp[list[j]]+min(j+1,t-j-1));
			}
			q.clear();
			q.push(N(0,dp[list[0]]));
			int adds=0;
			for(int j=1;j<2*t;j++){
				q.pop(j-t/2-1);
				adds++;
				ans=max(ans,dp[list[j]]+adds+q.ask());
				q.push(N(j,dp[list[j]]-adds));
			}
			dp[x]=max(dp[x],mx);
		}else if(!low[P[i]]){
			fa[P[i]]=x;
			dfs(P[i]);
			low[x]=min(low[x],low[P[i]]);
			if(low[P[i]]>dfn[x])ans=max(ans,dp[x]+dp[P[i]]+1),dp[x]=max(dp[x],dp[P[i]]+1);
		}else low[x]=min(low[x],dfn[P[i]]);
	}
	ans=max(ans,dp[x]);
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=0,k;i<m;i++){
		scanf("%d",&k);
		int last;
		scanf("%d",&last);
		for(int i=1,t;i<k;i++){
			scanf("%d",&t);
			add(last,t);add(t,last);
			last=t;
		}
	}
	dfs(1);
	printf("%d\n",ans);	
	return 0;
}

题解 BZOJ 1174 [Balkan2007]Toponyms

Tue, 10 Mar 2015 17:13:05 +0800

给你一个字符集合,你从其中找出一些字符串出来. 希望你找出来的这些字符串的最长公共前缀*字符串的总个数最大化.
这个问题很简单吧..用trie随便搞一搞。
但是Input是这么写的:
第一行给出数字N.N在[2,1000000] 下面N行描述这些字符串,长度不超过20000 。总字符不超过20000。

题解 BZOJ 2648 SJY摆棋子

Wed, 07 Jan 2015 17:12:59 +0800

双倍经验2716.
KD-Tree模板题。
注意一下"nth_element()"
参考： http://blog.csdn.net/jiangshibiao/article/details/34144829

code:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m,D,root,x,y;
struct N{
	int d[2];
	int l,r;
	int min[2];
	int max[2];
}a[1000005];
bool cmp(const N&a,const N&b){
	if(a.d[D]!=b.d[D]) return a.d[D]<b.d[D];
	return a.d[D^1]<b.d[D^1];
}
inline void up(int k,int s){
	a[k].min[0]=min(a[k].min[0],a[s].min[0]);
	a[k].min[1]=min(a[k].min[1],a[s].min[1]);
	a[k].max[0]=max(a[k].max[0],a[s].max[0]);
	a[k].max[1]=max(a[k].max[1],a[s].max[1]);
}
int bd(int l,int r,int d){
	D=d;
	int mid=(l+r)>>1;
	nth_element(a+l,a+mid+1,a+r+1,cmp);
	a[mid].min[0]=a[mid].max[0]=a[mid].d[0];
	a[mid].min[1]=a[mid].max[1]=a[mid].d[1];
	if(l!=mid) a[mid].l=bd(l,mid-1,d^1);
	if(r!=mid) a[mid].r=bd(mid+1,r,d^1);
	if(a[mid].l) up(mid,a[mid].l);
	if(a[mid].r) up(mid,a[mid].r);
	return mid;
}
const int inf=0x7fffffff;
int ans;
int getdis(int k){
	int ret=0;
	if(x<a[k].min[0]) ret+=a[k].min[0]-x;
	if(x>a[k].max[0]) ret+=x-a[k].max[0];
	if(y<a[k].min[1]) ret+=a[k].min[1]-y;
	if(y>a[k].max[1]) ret+=y-a[k].max[1];
	return ret;
}
void ask(int k){
	int d0=abs(a[k].d[0]-x)+abs(a[k].d[1]-y);
	if(d0<ans) ans=d0;
	int d1=(a[k].l)?getdis(a[k].l):inf;
	int d2=(a[k].r)?getdis(a[k].r):inf;
	if(d1<d2){
		if(d1<ans) ask(a[k].l);
		if(d2<ans) ask(a[k].r);
	}else{
		if(d2<ans) ask(a[k].r);
		if(d1<ans) ask(a[k].l);
	}
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d%d",&a[i].d[0],&a[i].d[1]);
	}
	root=bd(1,n,0);
	for(int i=0,t;i<m;i++){
		scanf("%d%d%d",&t,&x,&y);
		if(t==1){
			n++;
			a[n].min[0]=a[n].max[0]=a[n].d[0]=x;
			a[n].min[1]=a[n].max[1]=a[n].d[1]=y;
			int p=root,D=0;
			while(1){
				up(p,n);
				if(a[n].d[D]<a[p].d[D]){
					if(!a[p].l) {
						a[p].l=n;
						break;
					}else{
						p=a[p].l;
					}
				}else{
					if(!a[p].r) {
						a[p].r=n;
						break;
					}else{
						p=a[p].r;
					}
				}
				D^=1;
			}
		}else{
			ans=inf;
			ask(root);
			printf("%d\n",ans);
		}
	}
	return 0;
}

题解 BZOJ 1313 取整矩阵

Wed, 07 Jan 2015 17:12:48 +0800

题目

题里的矩阵是这样的：

A11 A12 A13… A1n-1 B1 A21 ………. A2n-1 B2 …… An-11 An-12….An-1n-1 Bn-1 C1 C2 C3… Cn-1 0

需要取整。每个元素大概有两种选择，一种是floor另一种是ceil。

这两种之间差了1。

网络流。

SPOJ QTREE系列

Tue, 06 Jan 2015 17:12:35 +0800

375.QTREE

给你一棵n个点$n-1$条边的树，支持修改一条边的边权、查询两点之间路径上的最大边权。

树链剖分+线段树基础题

code:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n;
int H[10005],X[20005],P[20005],tot, w[20005],E[20005],num[20005];
inline void add(int x,int y,int z,int no){
	P[++tot]=y;X[tot]=H[x];H[x]=tot;E[tot]=z;num[tot]=no;
}
int mx[10005<<2];
int to[10005],tt,siz[10005],son[10005],top[10005],dep[10005],f[10005],tim;
void dfs1(int x,int fa){
	siz[x]=1;son[x]=0;top[x]=0;
	int mx=0;
	for(int i=H[x];i;i=X[i]){
		if(P[i]==fa) continue;
		f[P[i]]=x;
		dep[P[i]]=dep[x]+1;
		w[P[i]]=E[i];
		to[num[i]]=P[i];
		dfs1(P[i],x);
		siz[x]+=siz[P[i]];
		if(siz[P[i]]>mx){
			mx=siz[P[i]];
			son[x]=P[i];
		}
	}
}
void dfs2(int x,int fa){
	if(!top[x]) top[x]=x;
	num[x]=++tim;
	E[tim]=w[x];
	if(son[x])top[son[x]]=top[x],dfs2(son[x],x);
	for(int i=H[x];i;i=X[i]){
		if(P[i]!=son[x]&&P[i]!=fa){
			dfs2(P[i],x);
		}
	}
}
int _c,_l,_r;
void ask(int o,int l,int r){
	if(_l<=l&&r<=_r){
		_c=max(mx[o],_c);
		return ;
	}else{
		int mid=(l+r)>>1;
		if(_l<=mid) ask(o<<1,l,mid);
		if(_r>mid) ask(o<<1|1,mid+1,r);
	}
}
int ask(int l,int r){
	_c=0;
	_l=l,_r=r;
	ask(1,1,n);
	return _c;
}
int Ask(int a,int b){
	int ret=0;
	while(top[a]!=top[b]){
		int x=top[a],y=top[b];
		if(dep[x]<dep[y]) swap(x,y),swap(a,b);
		ret=max(ask(num[x],num[a]),ret);
		a=f[x];
		
	}
	if(a!=b){
		if(dep[a]>dep[b]) swap(a,b);
		ret=max(ret,ask(num[a]+1,num[b]));
	}
	return ret;
}
char s[50];
void change(int o,int l,int r,int p,int k){
	if(l==r){
		mx[o]=k;
	}else{
		int mid=(l+r)>>1;
		if(p<=mid) change(o<<1,l,mid,p,k);
		else change(o<<1|1,mid+1,r,p,k);
		mx[o]=max(mx[o<<1],mx[o<<1|1]);
	}
}
void bd(int o,int l,int r){
	if(l==r){
		mx[o]=E[l];
	}else{
		int mid=(l+r)>>1;
		bd(o<<1,l,mid);
		bd(o<<1|1,mid+1,r);
		mx[o]=max(mx[o<<1],mx[o<<1|1]);
	}
}
int main(){
	scanf("%d",&tt);
	while(tt--){
		tim=tot=0;
		memset(H,0,sizeof H);
		scanf("%d",&n);
		for(int i=1,a,b,c;i<n;i++){
			scanf("%d%d%d",&a,&b,&c);
			add(a,b,c,i);
			add(b,a,c,i);
		}
		dep[1]=1;
		dfs1(1,-1);
		dfs2(1,-1);
		bd(1,1,n);
		int a,b;
		while(scanf("%s",s),s[0]!='D'){
			if(s[0]=='Q'){
				scanf("%d%d",&a,&b);
				printf("%d\n",Ask(a,b));
			}else{
				scanf("%d%d",&a,&b);
				change(1,1,n,num[to[a]],b);
			}
		}
	}
	return 0;
}

913.QTREE2

一棵树，两个操作。
1.询问a到b的最短路。
2.询问a到b的最短路上第k个点是啥。
倍增或者树链剖分。

SPOJ GSS 系列

Tue, 06 Jan 2015 17:12:20 +0800

SPOJ GSS系列是一系列序列维护的问题。
大部分用线段树，Splay等可以解决。

1043. GSS1

题意：给定长度为N的数串，M个询问查询[a,b]的的最大连续子段和。
线段树每个节点维护左端最大连续和，右端最大连续和，中间最大连续和。
就可以快速合并了。

题解 BZOJ 2527 [Poi2011]Meteors

Tue, 06 Jan 2015 17:12:09 +0800

一个叫做整体二分的东西。
单个国家可以二分。
为了节约时间，可以把情况相同的一同二分。
可以用树状数组做必要的统计。
时间复杂度不明。

code:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
int n,m;
vector<int> has[300005];
int p[300005];
typedef long long LL;
LL c[300005];
int k,T;
int l[300005],r[300005], a[300005];
int ans[300005];
int id[300005];
bool mark[300005];
int tmp[300005];
#define lowbit(x) ((x)&-(x))
void add(int x,int k){
	while(x<=m){
		c[x]+=k;
		x+=lowbit(x);
	}
}
LL ask(int x){
	LL ret=0;
	while(x>0){
		ret+=c[x];
		x-=lowbit(x);
	}
	return ret;
}
void Add(int x,int f){
	if(l[x]<=r[x]){
		add(l[x],f*a[x]);
		add(r[x]+1,-f*a[x]);
	}else{
		add(l[x],f*a[x]);
		add(1,f*a[x]);
		add(r[x]+1,-f*a[x]);
	}
}
void solve(int l,int r,int L,int R){
	if(L>R) return;
	if(l==r){
		for(int i=L;i<=R;i++){
			ans[id[i]]=l;
		}
	}else{
		int mid=(l+r)>>1;
		while(T<mid) T++,Add(T,1);
		while(T>mid) Add(T,-1),T--;
		for(int i=L;i<=R;i++){
			LL tmp=0;
			for(int j=0;j<has[id[i]].size();j++){
				tmp+=ask(has[id[i]][j]);
				if(tmp>=p[id[i]]) {
					mark[id[i]]=1;break;
				}
				
			}
		}
		int cc=0,tot=0;
		for(int i=L;i<=R;i++){
			if(mark[id[i]]){
				id[L+cc]=id[i];
				cc++;
			}else{
				tmp[tot++]=id[i];
			}
		}
		for(int i=L;i<L+cc;i++) mark[id[i]]=0;
		for(int i=0;i<tot;i++) id[L+cc+i]=tmp[i];
		solve(l,mid,L,L+cc-1),solve(mid+1,r,L+cc,R);
	}
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1,x;i<=m;i++) scanf("%d",&x),has[x].push_back(i);
	for(int i=1;i<=n;i++) scanf("%d",&p[i]);
	scanf("%d",&k);
	for(int i=1;i<=k;i++){
		scanf("%d%d%d",&l[i],&r[i],&a[i]);
	}
	k++;
	l[k]=1,r[k]=m,a[k]=0x3fffffff;
	for(int i=1;i<=n;i++) id[i]=i;
	solve(1,k,1,n);
	for(int i=1;i<=n;i++){
		if(ans[i]==k) puts("NIE");
		else printf("%d\n",ans[i]);
	}
	return 0;
}

NOIP2014 行记

Fri, 14 Nov 2014 15:53:38 +0800

不知道OI是啥或者信息学竞赛是啥的可以按Ctrl+W。

~~很早开始写的。。准备出分之后再发布。~~

谨以此文纪念我信息学竞赛的第一次正式考试。

考前一些天写的:

11.5

下午想转发一下去年的RP++的说说。发现去年没有发(或被删？)。。

题解 BZOJ 3707 圈地计算几何

Thu, 18 Sep 2014 15:53:33 +0800

题解：

简单的说就是给定平面上n个点，求这n个点组成三角形的最小面积。

如果分别枚举三个点的话是$O(n^3)$的，时间无法承受。

如果枚举了两个点a,b。设它们间的距离是L。如果以点a,b所在直线为y轴的话，与其他点所组成的三角形的面积$S=L*|x|/2$，x是其他点在这个坐标系中的横坐标。可以看出面积最小的就是离这个坐标系y轴最近的一个点。如果我们能够快速的得知最近的点的话，就可以将复杂度降低到$O(n^2)$。

题解 BZOJ 3714 [PA2014]Kuglarz

Sun, 14 Sep 2014 15:53:23 +0800

题解：

参考小胖的奇偶那道题。

那道题用的并查集维护。

如果知道 i..j 之间奇偶性的话，实际上知道的是 sum[j]-sum[i-1]的奇偶性（sum为前缀和）。

题解 BZOJ 3713 [PA2014]Iloczyn

Sun, 14 Sep 2014 15:53:16 +0800

题解:

鉴于F[44]> 1e9。于是可以把两两乘积算出来，枚举即可。

code:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<set>
//by zrt
//problem:
using namespace std;
int f[46];
set<int> s;
int main(){
    #ifdef LOCAL
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
    #endif
    f[0]=0;f[1]=1;
    for(int i=2;i<=45;i++){
        f[i]=f[i-1]+f[i-2];
    }
    int MAX=1e9;
    for(int i=0;i<=45;i++) s.insert(f[i]);
    for(int i=3;i<=45;i++){
        for(int j=i;j<=45;j++){
            if(f[i]*1LL*f[j]<=MAX){
                s.insert(f[i]*f[j]);
            }
        }
    }
    int t,x;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&x);
        if(s.count(x)){
            puts("TAK");
        }else{
            puts("NIE");
        }
    }
    return 0;
}

题解 BZOJ 3709 [PA2014]Bohater

Sun, 14 Sep 2014 15:53:08 +0800

题解：

若d[x]<a[x]，杀这个怪是有收益的，可以按d值从小到大杀。

若d[x]>a[x]，考虑反过来的过程，如果都杀完后血量是S，那么从后向前就是 S-a[i]+d[i]，相当于第一种情况，需要按a值从大到小杀。

题解 BZOJ 2442 [Usaco2011 Open]修剪草坪

Sun, 14 Sep 2014 15:53:01 +0800

题解：

定义F[i]为前 i-1 只奶牛工作效率的最大值。

sum[i]是Ei的前缀和。

有F[i]=max{F[j]+sum[i-1]-sum[j]} (i-j<k)

可以用单调队列维护 F[j]-sum[j] 。

题解 BZOJ 1735 [Usaco2005 jan]Muddy Fields 泥泞的牧场

Sun, 14 Sep 2014 15:52:55 +0800

题解：

不能盖住好地，那么宽为1的木板只能放在行、列连通块里。
所以行、列连通块对应左、右部中的点，泥地对应边。
求二分图最小覆盖就是答案。

二分图最小点覆盖==最大匹配

题解 BZOJ 1715 [Usaco2006 Dec]Wormholes 虫洞

Wed, 10 Sep 2014 15:52:47 +0800

题解：

题意是问你一个混合图是否存在负环。

spfa即可，开始时将所有点入队，求最短路，当最短路长度超过n时，说明有负环。

code:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
//by zrt
//problem:
using namespace std;
int n,m,w;
int H[505],X[6000],P[6000],E[6000];
int d[505],l[505];
int tot;
inline void add(int x,int y,int z){
    P[++tot]=y;X[tot]=H[x];H[x]=tot;E[tot]=z;
}
int tt;
queue<int> q;
int main(){
    #ifdef LOCAL
    freopen("in.txt","r",stdin);
//  freopen("out.txt","w",stdout);
    #endif
    scanf("%d",&tt);
    while(tt--){
        memset(H,0,sizeof H);
        tot=0;
        scanf("%d%d%d",&n,&m,&w);
        for(int i=1,x,y,z;i<=m;i++){
            scanf("%d%d%d",&x,&y,&z);
            add(x,y,z);
            add(y,x,z);
        }
        for(int i=1,x,y,z;i<=w;i++){
            scanf("%d%d%d",&x,&y,&z);
            add(x,y,-z);
        }
        memset(d,0,sizeof d);
        memset(l,0,sizeof l);
        while(!q.empty()) q.pop();
        for(int i=1;i<=n;i++) q.push(i);
        bool have=0;
        while(!q.empty()){
            int x=q.front();q.pop();
    //      printf("%d %d %d\n",x,d[x],l[x]);
            for(int i=H[x];i;i=X[i]){
                if(d[P[i]]>d[x]+E[i]){
                    d[P[i]]=d[x]+E[i];
                    q.push(P[i]);
                    l[P[i]]=l[x]+1;
                    if(l[P[i]]>n) {
                        have=1;goto ed;
                    }
                }
            }
        }
        ed:;
        if(have) puts("YES");
        else puts("NO");
         
    }
     
    return 0;
}

题解 COGS 577 蝗灾 [CDQ分治入门题]

Mon, 08 Sep 2014 15:52:30 +0800

题目链接

昨天mhr神犇，讲分治时的CDQ分治的入门题。

题意:

你有一个wxw正方形的田地。

初始时没有蝗虫。

给你两个操作：

1 x y z： (x,y)这个位置多了z只蝗虫。

题解 BZOJ 2878 [Noi2012]迷失游乐园

Mon, 08 Sep 2014 15:52:22 +0800

题目链接

题意：

给你一张n个点n-1或n条边的带权无向图。从每个点出发一直走下去，不能重复经过某个点。问走过的路径长度的数学期望是多少？ N<=100000。图中至多只有一个环，并且环长不超过20。

素数筛法

Fri, 15 Aug 2014 15:52:04 +0800

埃氏筛法：从2开始，找到第一个没有被筛的数，把它标记为素数，然后把它的2倍、3倍……筛掉。
复杂度O(nlogn)。

改进的埃氏筛法：从2开始，找到第一个没有被筛的数x，把它标记为素数，然后把它的x倍、x+1倍……筛掉。
复杂度O(nloglogn)。

题解 POJ 2425 A Chess Game

Fri, 15 Aug 2014 15:51:52 +0800

[题目链接][1]
题意：给定一个有向无环图（DAG），上面放有一些旗子，旗子可以重合，两个人轮流操作，每次可以把一个旗子从一个位置移动到相邻的位置，无法移动时输，询问先手是否必胜。

题解 POJ 1740 A New Stone Game

Fri, 15 Aug 2014 15:51:44 +0800

题目链接
题意：有n堆石子，两人轮流操作，每次每个人可以从一堆中拿走若干个扔掉（必须），并且可以从中拿走一些分到别的有石子的堆里（可选），当一个人不能拿时这个人输。给定状态，问是否先手必胜。

组合游戏学习

Thu, 14 Aug 2014 15:51:37 +0800

　　阅读了《由感性认识到理性认识——透析一类搏弈游戏的解答过程》、《解析一类组合游戏》、《组合游戏略述——浅谈SG游戏的若干拓展及变形》这三篇论文，对组合游戏以及SG函数有了更深的理解。这篇文章摘下了这三篇论文的部分重要内容，以及部分我对组合游戏的理解。

题解 POJ 2975 Nim 统计必胜走法个数

Thu, 14 Aug 2014 15:51:11 +0800

题目链接
题意介绍了一遍Nim取石子游戏，可以看上一篇文章详细介绍。
问当前状态的必胜走法个数，也就是走到必败状态的方法数。

我们设sg为所有个数的Xor值。
首先如果sg==0,它不可能有必胜走法，输出0.

题解 POJ 2505 A multiplication game 组合游戏

Thu, 14 Aug 2014 15:50:44 +0800

题目链接
题意：有一个数p=1，甲乙两人轮流操作，每次可以把p乘2-9中的一个数，给定一个n，当一个人操作后p>=n,那么这个人赢，问先手是否必胜。

题解 POJ 2484 A Funny Game

Thu, 14 Aug 2014 15:50:36 +0800

题目链接
题意：有n个硬币排成一圈，两个人轮流操作，每次可以取走一个或者相邻的连个硬币（只算最开始相邻的，取之后才相邻的不算），问先手必胜还是必败。

这个题可以证明若n>=3，则先手必败。
对称博弈
若n>=3，先手第一次必然把这个环拆成一个链，然后无论这条链长度的奇偶，后手总是可以把这条链分成两条相等的链，于是先手在一条链上做什么，后手就可以做什么。知道先手无法操作，后手胜。

题解 POJ 1067 取石子游戏

Thu, 14 Aug 2014 15:50:21 +0800

[题目链接][1]
有两堆石子，数量任意，可以不同。游戏开始由两个人轮流取石子。游戏规定，每次有两种不同的取法，**一是可以在任意的一堆中取走任意多的石子；二是可以在两堆中同时取走相同数量的石子。**最后把石子全部取完者为胜者。现在给出初始的两堆石子的数目，如果轮到你先取，假设双方都采取最好的策略，问最后你是胜者还是败者。（中文题面，感动ing）

题解 BZOJ 2337 [HNOI2011]XOR和路径

Wed, 13 Aug 2014 15:50:12 +0800

题目链接
题意：给一张n个点，m条边的带环无向连通图，求1-n路径上的权值Xor起来的期望。即无向有环图上求路径Xor期望。

因为是Xor操作，所以我们可以按二进制位分开考虑。
设考虑到了第k位，
那么现在这个图中的边权只有0或1。
期望问题反着设一般比较好列方程，所以我们设f[i]为从i这个点走到nXor值是1的期望，显然f[n]==0。
∑2^k*f[1]即为答案。
设P[i]为从i点出发的走每条边的概率，应该等于1/degree[i]。
我们可以列一组dp方程：
f[i]=∑P[i]*f[j](w[j,i]==0)+∑P[i]*(1-f[j])(w[j,i]==1)。
若w[j,i]==0则从这条边走过来Xor值不变，w[j,i]==1则从这条边走过来Xor值取反，(1-f[j])即为在jXor值是0的期望。

题解 BZOJ 1257 [CQOI2007] 余数之和

Wed, 13 Aug 2014 15:50:05 +0800

题目链接
题意：

　　给定n,k，求 ∑(k mod i) {1<=i<=n} 其中 n,k<=10^9。　　即 k mod 1 + k mod 2 + k mod 3 + … + k mod n的值。

我们先来看商之和。
　　给定n,k，求∑(k/i) {1<=i<=n} 其中/为整除。

题解 BZOJ 2242 [SDOI2011] 计算器

Tue, 12 Aug 2014 15:49:44 +0800

题目链接
noip级数论模版题了吧。
让求三个东西：

给定y,z,p,计算Y^Z Mod P 的值。
给定y,z,p，计算满足xy≡ Z ( mod P )的最小非负整数。
给定y,z,p，计算满足Y^x ≡ Z ( mod P)的最小非负整数。

其中P均为素数。
来分着处理。

题解 BZOJ 2725 故乡的梦

Fri, 01 Aug 2014 15:48:59 +0800

这个题在bzoj上好像是个权限题，想做的可以去Vani的博客下载测试数据。
这里有题面。

简单叙述一下题意：
给你一个n个点、m条边的带权无向图，S点和T点，询问Q次删一条给定的边的S-T最短路。

这几天做的POI的题

Sun, 22 Jun 2014 15:48:46 +0800

正规、严谨、精妙。 -POI

发现POI(波兰信息学奥赛)的题都很有意思。于是开刷bzoj上的poi题目(按ac人数降序。。)。顺手写一写题解，加深印象。

BZOJ 1103 : [POI2007]大都市meg

给一棵树，每次可以把树上的一些边标记了，问一个点与根之间需要走多少没有标记的边。

中考一周年纪念

Sat, 21 Jun 2014 15:48:32 +0800

“很难想象自己一年后的样子呢。”

如果没记错，一年前的今天应该是中考的日子，当时考完感觉很爽。(考完很爽的考试这辈子应该也没几次。) 过了几天后得知happyending，然后就玩了一个暑假，去了好多地方。

白书第五章(网络流)练习题

Sat, 21 Jun 2014 15:48:18 +0800

一切皆可网络流

图论题，转化到常见模型才是最重要的。

网络流

花了点时间把白书的网络流例题、习题都做了，覆盖的挺全的，主要是建图的技巧吧。一些简单的题解。

例题：

uva 11248 Frequency Hopping

网络扩容，bzoj上有一道类似的题，就是每次扩展最小割的边，看看够不够。两个优化：

HEOI2014 行记

Wed, 21 May 2014 15:47:22 +0800

HEOI2014 标程数据下载.

百度盘 http://pan.baidu.com/s/1qWx7YAo

有版权问题的话下面留言…

又到了一年一度的HEOI呢。我果然还是太弱了呢。

（省选是在5月17、18日举行）

Day0

报到日啊，省选今年在河北经贸大学举办，还算很便利的。

一个开始

Sat, 16 Feb 2013 21:28:57 +0800

恢复于： https://web.archive.org/web/20130420015442/http://zhangruotian.duapp.com/article/about-this-blog/

创建于2013年2月16日，mark一下，感谢BAE提供支持。

博客系统是用的是WordPress 3.5.1 原版，参考的“百度BAE安装原版WordPress教程_玻璃泉”这篇文章。

关于BAE的消息队列发送邮件这个东西折腾了半天，最终把评论提醒系统才弄好了。

Ruotian's I/O

2023年终总结

九寨沟之行

2022年终总结

Elden Ring

谷歌实习初体验

归档

💬 Chat

Peers

Message

Speak

Action

Source code

绕过清华校园网的二层隔离

PyScript 尝试

REPL

Python 协程 (async/await)

大端序与小端序的C++存取

大端序

关于

Links

Python 使用子进程调用函数并得到 exit code

Scrooge

🧶 Unravel Two

创建一个 Service/Systemctl 服务

使用有线网络连接树莓派

树莓派SD卡空间扩容

Apache 反向代理使用 Basic Auth

购入 Xbox Seires X

2021年终总结

奇思妙想

投掷计算器

Sublime Text 开启 GPU 渲染

Steve Jobs的经典演讲

Squid Game

Sony电视 Setup

BWT

一篇 Genome Biology 2009 的论文

滑板信用卡🛹

macOS使用Secoclient VPN时ssh, scp卡死的解决办法

Hello, macOS Big Sur!

使用体验

联想ThinkPad E480 CPU降频修复

重庆之行

修Macbook Pro 16' (2)

清华大学学费代交服务(2021)

☁️

清华大学美术学院毕业作品展

研究生

本科生

王兴已经两个月没发过饭否了

清洗HHKB键盘

小程序「THU洗衣」开通“流量主”

利用Service Worker实现CDN

无锡之行

小米对讲机1s

如何在清华吃汤圆？

小程序：THU洗衣

工商银行Clear卡

早安喵

宜家鲨鱼洗澡记

2020年终总结

《光遇》里那些好看的截图

脱单啦

国庆中秋，双节同庆

彩虹

洗衣机状态

修Macbook Pro 16'

开学典礼

新生报到

清华大学学费代交服务（2020）

观看了一次火箭发射

UOJ多评测机一键配置脚本

RSA来了，又走了

美国运通人民币信用卡

Apple Watch 使用初体验

大学四年回顾

Arduino连接两个MPU6050

关于 MPU 6050

Codeforces April Fools Day Contest 2020