对＠BeijingAir的一点怀疑

因为经常报出北京极端恶劣空气质量，美国驻华大使馆开的Twitter账号＠BeijingAir 最近备受关注。感谢这个账号，让我们了解了真实的情况。

我每天会多次查看＠BeijingAir，发现最近几天开始，准确地说是从12月14日18:00开始的那条推开始，每条空气质量指数(AQI)后都会加一句注解"at 24-hour exposure at this level"。例如某条推：

BeijingAir 12-20-2011; 11:00; PM2.5; 66.0; 151; Unhealthy (at 24-hour exposure at this level) // Ozone; 0.0; 0; No Reading

这句注解的意思是每个AQI都是在此级别下暴露24小时对人体健康的影响，从Good, Moderate, Unhealthy for Sensitive Groups到Unhealthy, Very Unhealthy和Hazardous皆如此。

可是实际情况呢？空气质量的影响怎么会24小时才有效果呢？若是很差的空气，恐怕几分钟就对健康有影响了，即使感觉不到。再说，恐怕没有哪个人会在户外连续暴露24小时之久；如果真如这条注解所言，这空气质量指数也就失去现实意义了。

联想到维基解密曾爆料过，说中国政府曾要求美国使馆让其北京空气质量检测数据的访问权限仅限使馆人员和美国公民。又考虑到最近媒体和民众对北京空气质量关注度很高，让政府狼狈应付 ，我谨慎怀疑这条"at 24-hour exposure at this level"是美国驻华使馆迫于中国政府压力做出妥协而加上去的 。正好看到Wikileaks的这条Cable最后说：

EmbOffs promised to relay MFA's concerns to relevant sections of the Embassy to see if a compromise could be identified.

难不成这就是妥协的结果？总之，我是不信"at 24-hour exposure at this level"这句话的，北京的空气该躲还得躲。

--
Sent from my mobile.

手机选号

最近要弄个联通的186手机号，想选个酷一点的号。首先想到的是平方数。186开头的平方数并不多，可惜可选的号码里面没有一个是平方数。后面又探索了一些特殊的数，但是这些数也很稀有，比如立方数、斐波那契数、三角数等，所以也没有找到。

后来就在圆周率Pi里面找。北京联通现在可以选的186号只有1860和1861开头的。在Pi的小数点后10亿位里面找1860和1861开头的手机号（11位数），根据联通可用号段进一步筛选，剩下了11万多个，一行一号地输出。我又把北京联通选号页面上的100多页近7000个手机号人工复制粘贴下来，也修改成每行一个的文件。然后开始找，大概花了两个半小时找了出来，只有64个。搜索的代码绝对是糙、快、猛：

#! /bin/bash
while read -r line1; do
    while read -r line2; do
        [[ $line1 == $line2 ]] && echo "$line1"
    done <pi
done <unicom

“快”是说写得快，不是算得快。里面pi就是Pi里面找出来的186[01]的号，unicom就是北京联通网站上的所有可选号。

现在放的号里面，不用预存话费的很少有不带4的。最后选了一个0最多的，因为我觉得数字0会让手机号显得比较“轻盈”。我的这个号在pi小数点后1亿多位就出现了。我不选要预存上千话费的“靓号”，比如186 0050 4000就在pi小数点后第881,551,945出现。

在找186号之前，我找了一下在用的电信133号（以前也是联通的），竟然在9亿多位的时候出现了，要知道总共只有10亿位，差点就出现不了了。总之，不是每个手机号都有这么幸运的。