鬼の领地 » 小程序

关于饭否的两个小脚本

upsuper — Wed, 02 Feb 2011 04:06:24 +0000

前几天因为某些原因，我把饭否上所有的好友和关注者全部清空了。当然，如果没有程序的帮忙，估计还不等我删完我也就后悔了。

我没有那么狠心的把饭否的消息给清空，因为消息是不可恢复的（而且也太多了），但好友和关注者是可以的。做事情都给自己留后路显然是我一贯的风格，不然的话，我大概早从我家阳台跳下去了……

两段脚本都不长，第一段是备份饭否的好友列表和关注者列表的，做的毫无泛用性，因为 bash 编程我并不很熟，只知道可以用 wget 来抓取。本来这种事情其实可以直接用 API 抓的，不过我想抓下我能直接看的东西，所以最终还是导出了 Cookie 抓网页。

第一段代码如下：

#!/bin/bash
 
mkdir friends
for i in {1..11}
do
    wget -k -e robots=off --load-cookies cookies.txt -P friends/ \
        http://fanfou.com/friends/upsuper/p.$i
    sed -i -e "s/http:\/\/fanfou\.com\/friends\/upsuper\///g" \
        "friends/p.$i"
    sed -i -e "s/http:\/\/fanfou\.com\/friends/p.1/g" \
        "friends/p.$i"
    sed -i -e "s/http:\/\/fanfou\.com\/followers/..\/followers\/p.1/g" \
        "friends/p.$i"
done
 
mkdir followers
for i in {1..8}
do
    wget -k -e robots=off --load-cookies cookies.txt -P followers/ \
        http://fanfou.com/followers/upsuper/p.$i
    sed -i -e "s/http:\/\/fanfou\.com\/followers\/upsuper\///g" \
        "followers/p.$i"
    sed -i -e "s/http:\/\/fanfou\.com\/friends/..\/friends\/p.1/g" \
        "followers/p.$i"
    sed -i -e "s/http:\/\/fanfou\.com\/followers/p.1/g" \
        "followers/p.$i"
done
 
for f in friends/* followers/*
do
    sed -i -e "s/http:\/\/fanfou\.com\/friends/..\/friends/g" \
        "followers/p.$i"
    sed -i -e "s/http:\/\/fanfou\.com\/followers/..\/followers/g" \
        "friends/p.$i"
done

想用的话，其中前两个 for 循环的数字必须改成自己的实际情况，里面的 upsuper 也要全数替换为自己的饭否ID，还有一个 cookies.txt 文件，可以用 Firefox 的 Cookie Exporter 插件从浏览器导出。

这个脚本非常值得赞美的一点是自动修复了所有的链接，使得在本地的页面全部可以通过链接进入，不在本地的页面也可以链接到饭否。这也是我调了很长世间的东西……

第二段脚本就是颇具破坏性的了：

#!/usr/bin/env python
# - * - coding: utf8 - * -
 
import sys
import json
 
from socket import timeout
from urllib import urlencode
from base64 import b64encode
from getpass import getpass
from httplib2 import Http
 
username = raw_input('Username: ')
password = getpass('Password: ')
auth = 'Basic ' + b64encode('%s:%s' % (username, password))
headers = { 'Authorization': auth }
 
h = Http(timeout=1)
sys.stdout = sys.stderr
 
# Delete friends
print 'Delete friends:'
resp, content = h.request('http://api.fanfou.com/friends/ids.json', headers=headers)
friends = json.loads(content)
for friend in friends:
    query_str = urlencode({ 'id': friend.encode('utf8') })
    print friend.encode('utf8'),
    try:
        resp, content = h.request('http://api.fanfou.com/friendships/destroy.json?' + query_str, 'POST', headers=headers)
    except timeout:
        print 'timeout'
    else:
        print resp.status
 
# Delete followers
print 'Delete followers:'
resp, content = h.request('http://api.fanfou.com/followers/ids.json', headers=headers)
followers = json.loads(content)
for follower in followers:
    query_str = urlencode({ 'id': follower.encode('utf8') })
    print follower.encode('utf8'),
    try:
        resp, content = h.request('http://api.fanfou.com/blocks/create.xml?' + query_str, 'POST', headers=headers)
    except timeout:
        print 'timeout',
    else:
        print resp.status,
    try:
        resp, content = h.request('http://api.fanfou.com/blocks/destroy.xml?' + query_str, 'POST', headers=headers)
    except timeout:
        print 'timeout'
    else:
        print resp.status

这段 Python 脚本其实颇具我写饭否应用的风格。这段脚本还是比较有泛用性的，进入的时候会提示输入饭否id和密码，然后他也会自动备份一个 JSON 格式的饭否好友列表和关注者列表。只是处于周全我才做的这一步，其实是完全没有必要的，因为我已经备份过了，而且我也不准备让程序帮我全部还原……

正好清理一下两个列表……

原理就不多做介绍了。第二段脚本在运行开始清理好友的时候，可能会连续出现 timeout 的情况不必在意，因为饭否似乎会在这个 API 执行的时候处理 timeline 的缓存。不过据我观察，脚本运行结束后，似乎会有部分好友留下，数量对我来说少于十个，可以快速手动清除……

出于观察输出的考虑，我没有用多线程。脚本删除关注者的原理是加入黑名单再取消黑名单，如果加入了黑名单在取消的时候没有成功，那就意味着误加了，这是找不到的，因为黑名单的列表似乎无法察看……

就这样了……我现在已经逐渐的开始恢复饭否活动了……不过还不准备完全回到过去……以上

突然觉得，第一段脚本可以用 make 实现？

Rhythmbox 封面艺人插件的中国版

upsuper — Tue, 19 Oct 2010 05:17:05 +0000

上一篇日志里面我做了一个从 VeryCD 为 Rhythmbox 提取专辑封面的程序，后来想到豆瓣其实也有这个东西，于是又做了一个豆瓣的。然后想，既然这样我干脆把原来的封面艺人插件 fork 出来，自己开一个新项目，专门做封面艺人插件的 China Mod~

现在已经开好了，就在 Google Code 上：http://code.google.com/p/rhythmbox-artdisplay-cn/

其实说 Mod 也不算，就是加了几个扩展进去……主要是考虑到如果不 fork 似乎很单独安装，所以就 fork 出来了……那个插件看过去也很长时间没有更新了呢？

第一个版本0.1可以点击这里下载：rhythmbox-artdisplay-cn-0.1.tar.bz2，暂时还没有做任何安装包，deb 包和 arch 的包应该会随后奉上~

安装很简单，常规的 make 和 make install 就可以了~

如果各位有什么改进意见，欢迎提出~

利用 VeryCD 给 Rhythmbox 添加专辑封面

upsuper — Mon, 18 Oct 2010 11:54:14 +0000

我们知道 Rhythmbox 里面有一个插件叫做“封面艺人”，可以通过互联网下载对应专辑的封面图片。可是，那个插件使用的似乎是美国的 Amazon，中文歌曲基本找不到对应封面，怎么办呢……

想到我的很多专辑都是从 VeryCD 上下载的，那些下载页面通常都会有专辑封面图片，于是我就写了个封面艺人的扩展，让其可以从 VeryCD 上搜索专辑，并下载相应专辑封面。这个扩展可以在这里下载：verycd-cover.tar.bz2 (5KB)

下载解压后，运行其中的 compile.py，然后把除了 compile.py 以外的四个文件拷贝到 Rhythmbox 的封面艺人插件的目录就可以了，在 Ubuntu 下，那个目录是 /usr/lib/rhythmbox/plugins/artdisplay。

已知问题

写这个总共没花多少时间，所以也不能指望十分精品。实现的基本方式就是通过 VeryCD 的搜索 API 获得匹配的 ID，然后进去寻找是否有封面图片。

这种做法有很多局限性，比如可能有的专辑不是独立为一个页面的，这样可能获取的封面是错误的……

此外，由于 Rhythmbox 对于不写艺术家一项的歌曲，艺术家会返回“未知”，如果直接用“未知”去找，十有八九是找不到的，所以我在代码里着了点特判，不知道有没有更优美的方法……

话说，我发现豆瓣音乐似乎也有不少资源，这样还可以写一个从豆瓣获取封面的扩展~

寻找更高仿的 ID 第二季

upsuper — Sat, 11 Sep 2010 07:31:17 +0000

继上一篇文章之后，我又下大力气对这个程序做了许多修改，在精确度和速度方面似乎都有些许提高。在此推出第二季~

使用真正的12px宋体

在上一次的程序中使用的 PIL 似乎是因为不支持宋体 ttc 文件中对于小字体下优化的点阵形式，才在选择小于 19px 的字号时不能正确渲染汉字。考虑到这一点，我就想到把 ttc 文件里面 12px 的点阵字体单独提取出来使用，毕竟贴吧上面显示 ID 都是用这个字号显示的。

使用 FontForge 提取出来了 simsun-12.bdf 文件，就是宋体 12px 下的点阵。参考 PIL 的手册，发现 PIL 不能直接使用 .bdf 文件，需要使用一个叫做 pilfont 的脚本转换成专有的 .pil 文件才行。我想转换就转换呗。simsun-12.bdf 一个 2.4MB 的文件，转换完就剩不到 100KB，我就觉得肯定有问题，用 PIL 导入，发现还是不能渲染中文。后来知道，这个 .pil 文件根本不支持非拉丁字母的字符，它的储存空间限定了 256 个字符……

无奈了，这意味着 PIL 完全无法支持中文点阵了……

当然，办法总归是有的，那就是——抛弃 PIL！为什么我能有这样的想法呢，因为看到 .bdf 文件是 UNIX 标准的。UNIX 标准意味着什么呢？记不记得 UNIX 有一个非常好的传统叫做，尽量使用纯文本。是的，这虽然让有些文件会变得太大，不过同时也让这些东西更容易被其他程序读取，而 .bdf 恰好即使这么一种文件。

这样读取 .bdf 点阵字体文件的程序自己写不就好了，什么额外的库都不需要……当然，纠错性极弱就是了~

line_count = 0
 
def read_split(f):
    global line_count
    line_count += 1
    line = f.readline()
    if not line:
        return False
    return line[:-1].split()
 
with open('simsun-12.bdf', 'r') as f:
    chars = {}
    line = ['']
    try:
        # 获得字符总数
        while line[0] != 'CHARS':
            line = read_split(f)
        # 读取所有字符
        for i in range(int(line[1])):
            assert read_split(f)[0] == 'STARTCHAR'
            line = read_split(f)
            # 编码
            assert line[0] == 'ENCODING'
            char = unichr(int(line[1]))
            line = read_split(f)
            # 绘制参数
            assert line[0] == 'BBX'
            width, height, x1, y1 = [int(x) for x in line[1:]]
            x0, y0 = (1 + x1, 12 - y1 - height)
            # 准备绘制文字
            base_image = [[0 for x in range(14)] for y in range(14)]
            # 读取并绘制文字
            line = read_split(f)
            assert line[0] == 'BITMAP'
            for y in range(y0, y0 + height):
                line = read_split(f)[0]
                bits = int(line, 16)
                bits >>= len(line) * 4 - width
                for x in range(x0 + width - 1, x0 - 1, -1):
                    base_image[y][x] = bits & 1
                    bits >>= 1
            chars[char] = base_image
            # 结束这个文字
            line = read_split(f)
            assert line[0] == 'ENDCHAR'
        line = read_split(f)
        assert line[0] == 'ENDFONT'
    except AssertionError:
        print 'line', line_count
        raise

这样所有的字符就被读入，并变成一个单色像素二位数组了。

当然，这个性能很低，在我的机器上转换读取文件的2W+字符大概需要 18s，这可能也是为什么 PIL 要选择进行转换。事实上，使用纯文本储存一直以来都给 UNIX 风格的这一类软件带来一定性能缺陷。不过其实，这很值得，因为方便。

不过这个时间确实是太长了，更何况我们到目前为止还什么都没处理。怎么办呢？两个想法：一、优化代码；二、保存处理的数据。

第一种，基本上是没什么希望了，而且即使能优化，估计效果也不会太好，可能省个几秒封顶了。第二个显然不是个坏想法。

Python 在数据的持续化方面还是有很多现成的东西的，比如 pickle 什么的。不过那个速度太慢，而且是纯文本！好吧，偶尔我也会不喜欢纯文本，因为在这里意义不大……因此选择了 marshal。marshal 也是一个用于数据持续化的库，不过仅能对 Python 的内部类型进行。我会看中它最重要的原因就是它的应用范围极其有限，只能持续化内部类型。如果一个 Python 标准库，它有很明显的限制，却没有标明不推荐或在新版中被剔除，说明它必然有一个其他库不可及的优势。对于 marshal，我猜它的优势就是效率。

使用 marshal 就很简单了……

62
63

import marshal
marshal.dump(chars, open('base_image', 'wb'))

这样后面的处理不需要不断重复这个低效的步骤了~

处理文字图像

原来是用 PIL 处理文字图像，现在抛弃 PIL 了，就得自己写了……不过这样也很好，自由发挥的空间很大了~

我猜用的是 Matrix67 大牛的那种在附近留阴影的方法，不过似乎我写的不够好就是了，怎么测试效果都不大理想。除此之外，纯 Python 实现的算法效率和 PIL 这种包装还是没得比，很简单的算法却慢的不得了……

下面是目前的处理代码：

#!/usr/bin/python
# - * - coding: utf8 - * -
 
import marshal
import copy
 
base_chars = marshal.load(open('base_image', 'rb'))
chars = {}
for char, image in base_chars.iteritems():
    new_image = []
    for y in range(14):
        new_row = []
        for x in range(14):
            if image[y][x]:
                value = 81
            else:
                value = 0
            if y > 0 and image[y-1][x]:
                value += 4
            if y < 13 and image[y+1][x]:
                value += 4
            if x > 0 and image[y][x-1]:
                value += 4
            if x < 13 and image[y][x+1]:
                value += 4
            new_image.append(value)
    chars[char] = new_image
 
marshal.dump(chars, open('advanced_data', 'wb'))

处理效果不是很理想就是了，耗时大概也是 30s+。

寻找相似字符

其实这个部分就是一样的了……直接贴代码好了……

#!/usr/bin/python
# - * - coding: utf8 - * -
 
import marshal
 
chars = marshal.load(open('advanced_data', 'rb'))
 
def image_diff(image1, image2):
    ret = 0
    for v1, v2 in zip(image1, image2):
        ret += (v1 - v2) ** 2
    return ret
 
remember_number = 5
try:
    searched = marshal.load(open('searched_chars', 'rb'))
except IOError:
    searched = {}
while True:
    input_string = unicode(raw_input('ID: '), 'utf8')
    if not input_string:
        break
    for char in input_string:
        if char not in chars:
            continue
        if char not in searched:
            diff_data = []
            image = chars[char]
            for c, v in chars.iteritems():
                if c == char:
                    continue
                diff_data.append((image_diff(v, image), c))
            diff_data.sort()
            searched[char] = diff_data[:remember_number]
        print char,
        for item in searched[char]:
            print u'({0}, {1})'.format(*item),
        print
    print
 
marshal.dump(searched, open('searched_chars', 'wb'))

和上次不同的是查找过的字会被保存下来，效率可以高一些……

继续改进

现在的主要问题就是如何提高相似度的识别精度了……目前的想法是通过逐像素比对测试两个字的相似度，最多加一些模糊化什么的处理。doggy 提出一个想法是计算连通区域面积的比例，我个人认为不大可行……我的想法是识别文字的笔画，把文字的骨架弄出来，然后对比什么的，可能效果更好吧？

不知道各位还有没有其他什么想法？

下载

最后是完整版的目前的程序以及已经生成的数据，各位有兴趣可以直接拿来试用了……find-similar-character.tar.bz2 (1.7MB)

参考资料

寻找更高仿的 ID

upsuper — Thu, 09 Sep 2010 10:21:43 +0000

今天大学军训完了，不想做什么正经事，就想到前一段时间想做的寻找相似汉字的程序，用以寻找更高仿的贴吧 ID。用程序来寻找相似汉字，从另一个角度，也是从 Matrix67 大牛的一篇日志里得到的启发。不过 Matrix67 大牛使用的是 Mathematica 来寻找，我不大会 Mathematica，就想用我熟悉的 Python 来解决，毕竟 Python 是一个很强大的东西~

其实寻找的思路很简单，就是把某个汉字当作图片弄出来，让后对比两个图片的相似程度。因此做这个程序的第一步就是研究如何用 Python 处理图片和文字。Python 有一个非常著名的第三方库，名叫 Python Imaging Library，简称 PIL，就是专门用来处理图片的。

文字 to 图像

PIL 可以很轻松的将文字转换为图像，并且提供了虽然不能说是强大，但暂时够用的图像处理函数。

处理文字生成的图像，显然和彩色没有太大关系，因此可以使用灰度图像节省计算需要的空间和时间。此外我们知道，文字到图像有一个中间媒介，就是字体。因此我们必须先确定我们需要的字体和大小。用 Firebug 考查贴吧用于显示 ID 的字体和大小，发现是宋体 12px。这样参考 PIL 的手册，就有了最基本的文字到图片的代码了：

import Image, ImageDraw, ImageFont
image = Image.new('L', [14, 14], 255)
draw = ImageDraw.Draw(image)
font = ImageFont.truetype('simsun.ttc', 12)
draw.text((1, 1), u'好', font=font)
image.save('hw.png')

看过去这个代码是正确的，不过在查看 hw.png 后，我发现并没有出现预期的“好”字，而是一堆混乱的东西，不知何故。后经过不断实验，发现只有当初始的文字大小设定为19或以上时，文字才可以被正确地绘制出来，_-b

于是不得不修改这段代码，让它先绘制一个大的，再缩放成需要的大小：

import Image, ImageDraw, ImageFont
image = Image.new('L', [28, 28], 255)
draw = ImageDraw.Draw(image)
font = ImageFont.truetype('simsun.ttc', 24)
draw.text((2, 2), u'好', font=font)
# 抗锯齿方式缩放
image = image.resize((14, 14), Image.ANTIALIAS)
image.save('hw.png')

效果还不错。

枚举汉字

然后来找一种方式来枚举汉字。参考维基百科的 GB2312 条目，确定了 GB2312 中汉字区的编码，为第一个字节 0xB0-0xF7，第二个字节 0xA1-0xFE，其中 D7FA-D7FE 是空的，这样共有6763个汉字拿来搞。

其实还是挺少的就是了……

于是枚举汉字的代码：

for i in range(0xB0, 0xF8):
    chr_i = chr(i)
    for j in range(0xA1, 0xFF):
        if i == 0xD7 and j >= 0xFA and j <= 0xFE:
            continue
        char = unicode(chr_i + chr(j), 'gb18030')
        print char,

初步成果

现在结合上面两段，得到了下面的初步代码：

#!/usr/bin/python
# - * - coding: UTF-8 - * -
 
import Image, ImageDraw, ImageFont, ImageFilter
 
font = ImageFont.truetype('simsun.ttc', 24)
 
def get_char_data(char):
    image = Image.new('L', [28, 28], 255)
    draw = ImageDraw.Draw(image)
    draw.text((2, 2), char, font=font)
    image = image.resize((14, 14), Image.ANTIALIAS)
    return list(image.getdata())
 
def diff_of_data(data1, data2):
    ret = 0
    for i in zip(data1, data2):
        ret += abs(i[0] - i[1])
    return ret
 
chars = {}
for i in range(0xB0, 0xF8):
    chr_i = chr(i)
    for j in range(0xA1, 0xFF):
        if i == 0xD7 and j >= 0xFA and j <= 0xFE:
            continue
        char = unicode(chr_i + chr(j), 'gb18030')
        chars[char] = get_char_data(char)
 
remember_number = 5
input_string = unicode(raw_input('ID: '), 'utf8')
for input_char in input_string:
    if input_char not in chars:
        print input_char
        continue
    input_data = chars[input_char]
    diff_list = []
    for char, data in chars.iteritems():
        if char == input_char:
            continue
        diff_list.append((diff_of_data(data, input_data), char))
    diff_list.sort()
    print input_char,
    for diff, char in diff_list[:5]:
        print u'({0}, {1})'.format(diff, char).encode('utf8').ljust(11),
    print

发现效果还有待提高，并且每次测试一个新 ID 都要做一次预处理，每次预处理大概需要6s左右。

目前代码

对上面的代码做了一些修改，包括对文字图像做一些处理，以使其查找出的相似文字能更接近人的感觉，并且让一次预处理可以做多次比对，得到了下面代码：

#!/usr/bin/python
# - * - coding: UTF-8 - * -
 
import Image, ImageDraw, ImageFont, ImageFilter
 
font = ImageFont.truetype('simsun.ttc', 24)
 
def get_char_data(char):
    image = Image.new('L', [28, 28], 255)
    draw = ImageDraw.Draw(image)
    draw.text((2, 2), char, font=font)
    image = image.resize((14, 14), Image.ANTIALIAS)
    # XXX 期待更好的处理方式
    image = image.filter(ImageFilter.SHARPEN).filter(ImageFilter.SMOOTH)
    return list(image.getdata())
 
def diff_of_data(data1, data2):
    ret = 0
    for i in zip(data1, data2):
        ret += abs(i[0] - i[1])
    return ret
 
chars = {}
for i in range(0xB0, 0xF8):
    chr_i = chr(i)
    for j in range(0xA1, 0xFF):
        if i == 0xD7 and j >= 0xFA and j <= 0xFE:
            continue
        char = unicode(chr_i + chr(j), 'gb18030')
        chars[char] = get_char_data(char)
 
remember_number = 5
while True:
    input_string = unicode(raw_input('ID: '), 'utf8')
    if not input_string:
        break
    for input_char in input_string:
        if input_char not in chars:
            print input_char
            continue
        input_data = chars[input_char]
        diff_list = []
        for char, data in chars.iteritems():
            if char == input_char:
                continue
            diff_list.append((diff_of_data(data, input_data), char))
        diff_list.sort()
        print input_char,
        for diff, char in diff_list[:5]:
            print u'({0}, {1})'.format(diff, char).encode('utf8').ljust(11),
        print
    print

继续改进

目前这个程序测试的效果还行，不过还有待进一步改进。可改进的地方包括上面标出 XXX 的部分，即对文字图像进行处理的部分，或许可以有更好的处理方式。除此之外，如果能直接绘制出 12px 的字，必然也会有更好的效果。如果增大字库，从 GB2312 增加到 GB18030 的字库，可以用来对比的文字也会更多，也就更可能找到相似的字了~

这个，还有继续改进的余地嗯~有时间就继续努力咯……

参考资料

百度贴吧的世界杯竞猜

upsuper — Thu, 17 Jun 2010 08:12:24 +0000

貌似很久没发文章了，无聊冒出来发一个……

这次发的东西是很无聊的，大家知道最近世界杯在踢。不过呢，作为一个对体育几乎毫不关心的人，这事基本上也跟我没什么关系……不过我要说的呢，是关于世界杯竞猜的“赌博”游戏，貌似不止百度再做，还有网易、校内什么的都在搞。不过鉴于百度这个，比较经常上，我就也来参加了~

当然，发在这里的东西总归要有一些技术性……这次也不例外的……

我看到网易的竞猜不仅有各种形式，而且最重要的是有赔率！而百度只显示投注金额，显然不够专业，于是我就拿 Greasemonkey 插件，可以在百度投注金额的下面显示当前赔率~就像下面这样：

安装了 Firefox 的 Greasemonkey 插件的，可以点击这里安装这个脚本：tieba_guess.user.js

哦，有人大概会注意到我在前面那张图里面下注给了韩国。倒不是我多喜欢棒子，只是因为我几乎完全不懂足球，所以就根据赔率压冷门……话说压冷门有的时候其实很划算的，比如昨晚我压瑞士，赔率40！赚翻了~强大吧~

QR code 生成器

upsuper — Fri, 30 Apr 2010 15:04:52 +0000

前几天在找 Android 应用的时候，又见到了 QR code。说到 QR code，记得2006年，我还在维基活跃的时候，有一个人在维基的 QQ 群上发了这个的东西，然后大家就觉得很好玩。很早的事情了。寒假回福州的时候，在乌山公园，不少地方也有这个二维码，是移动设置的“物联网试验区”。

二维码着实是一个很有趣的东西。据我去过日本的同学讲，这个东西在日本已经十分普及了。在维基上查了一下，QR code 是开放的，有专利但没有被行使。而且识别算法似乎已经十分稳定成熟了。Android 手机里面安装一个 Barcode Scanner 就可以利用摄像头扫描二维码了，非常方便~

二维码可以干嘛呢？它可以储存网址、储存电话、储存名片信息，扫描一下就可以读取出来，将人从人工的复制方式中解脱出来~

于是我便也做了一个 CGI 的 QR code 生成器。

这个生成器用了开源的 libqrencode 和 libmagick 这两个库，都可以直接在 Ubuntu 里面 apt-get 下来，第一个是生成二维码数据，第二个是输出为图片。使用二进制 CGI 还涉及了我可爱的服务器商的支持问题等等……

生成器源代码

核心程序没什么好说的，其实就是查文档的问题而已。参考了一下 libqrencode 和 libmagick 的文档就可以了。

直接贴代码：

#include 
#include 
#include 
#include 
#include 
#include 
 
#define DEFAULT_TEXT "http:/""/upsuper.org/"
 
int main(int argc, char **argv) {
    int i, j;
    int scale;
    QRecLevel level;
    char *text;
 
    // 设置默认参数
    scale = 8;
    level = QR_ECLEVEL_Q;
    text = (char *)malloc(strlen(DEFAULT_TEXT) + 1);
    strcpy(text, DEFAULT_TEXT);
 
    // 读取参数
    const char *orig;
    char *query_string;
    char *p;
    int len;
    orig = getenv("QUERY_STRING");
    len = strlen(orig);
    query_string = (char *)malloc(len + 2);
    *query_string = '&';
    strcpy(query_string + 1, orig);
    for (p = query_string; *p; ++p) {
        if (*p != '&')
            continue;
        if (strncasecmp(p + 1, "scale=", 6) == 0) {
            // 读取缩放级
            scale = atol(p + 7);
        } else if (strncasecmp(p + 1, "text=", 5) == 0) {
            // 读取文本
            char *p_last, *p_src, *p_dest;
            for (p_last = p + 6; *p_last && *p_last != '&'; ++p_last);
            text = (char *)realloc((void *)text, p_last - p - 6 + 1);
            // urldecode
            for (p_src = p + 6, p_dest = text;
                    p_src < p_last;
                    ++p_src, ++p_dest) {
                if (*p_src == '+') {
                    *p_dest = ' ';
                } else if (*p_src == '%') {
                    int code;
                    if (sscanf(p_src + 1, "%2x", &code) != 1)
                        code = '?';
                    *p_dest = code;
                    p_src += 2;
                } else {
                    *p_dest = *p_src;
                }
            }
            *p_dest = '\0';
        } else if (strncasecmp(p + 1, "level=", 6) == 0) {
            // 读取纠错级别
            switch (*(p + 7)) {
                case 'l': case 'L': level = QR_ECLEVEL_L; break;
                case 'm': case 'M': level = QR_ECLEVEL_M; break;
                case 'q': case 'Q': level = QR_ECLEVEL_Q; break;
                case 'h': case 'H': level = QR_ECLEVEL_H; break;
            }
        }
    }
 
    // 生成 QR code
    QRcode *qrcode;
    int width;
    qrcode = QRcode_encodeString8bit(text, 0, level);
    width = qrcode->width;
    // 初始化 Image Magick
    ImageInfo *image_info;
    ExceptionInfo exception;
    InitializeMagick(*argv);
    GetExceptionInfo(&exception);
    // 初始化图片
    Image *image, *resized_image;
    image_info = CloneImageInfo((ImageInfo *)NULL);
    image = AllocateImage(image_info);
    image->columns = image->rows = width;
    // 生成图片
    PixelPacket *pixel_packet;
    int t_bit;
    for (i = 0; i < width; ++i) {
        pixel_packet = SetImagePixels(image, 0, i, width, 1);
        if (! pixel_packet)
            break;
        for (j = 0; j < width; ++j, ++pixel_packet) {
            t_bit = ! (qrcode->data[i * width + j] & 1);
            pixel_packet->red   = MaxRGB * t_bit;
            pixel_packet->green = MaxRGB * t_bit;
            pixel_packet->blue  = MaxRGB * t_bit;
        }
        if (! SyncImagePixels(image))
            break;
    }
    QRcode_free(qrcode);
    // 处理图片
    resized_image = ScaleImage(image,
            width * scale, width * scale, &exception);
    DestroyImage(image);
    CompressImageColormap(resized_image);
    // 输出图片
    puts("Content-Type: image/png");
    puts("");
    strcpy(resized_image->filename, "png:-");
    WriteImage(image_info, resized_image);
    // 清理
    DestroyImage(resized_image);
    DestroyImageInfo(image_info);
    DestroyExceptionInfo(&exception);
    DestroyMagick();
 
    return 0;
}

话说最近喜欢用 C，C++ 什么的最讨厌了……

另外写了个 Makefile：

LDFLAGS:=-lqrencode `Magick-config --ldflags --libs`
CFLAGS:=`Magick-config --cflags`
 
all: qrcodetest
qrcodetest: qrcodetest.c

下面就是把这个 cgi 部署到服务器的问题了……

部署到服务器

我的这个空间商吧，我猜他装了 libmagick 的库，不过肯定没装开发库，也就是 libmagick-dev 这种东西，这样就编译不了了，比较麻烦。探究了一下发现安装了 gcc~很是高兴啊~而且 cgi 支持 bash 的。

这样就好玩了~我在我的访问不到的根目录下建立了一个 usr 文件夹，就已这个作为存在的核心了~

我写了 bash cgi 程序直接从服务器上获取源代码包，在服务器上解压编译 ImageMagick，也就是下面这段代码：

#!/bin/bash
 
export PATH=$PATH:/xxxx/upsuper.org/usr/bin
export CPATH=/xxxx/upsuper.org/usr/include
export LIBRARY_PATH=/xxxx/upsuper.org/usr/lib
 
echo Content-Type: text/plain
echo
wget http://downloads.sourceforge.net/project/imagemagick/ImageMagick/00-6.6.1/ImageMagick-6.6.1-5.tar.bz2
mkdir imagick
mv ImageMagick-6.6.1-5.tar.bz2 imagick
cd imagick
tar -jxvf ImageMagick-6.6.1-5.tar.bz2
cd imagick/ImageMagick-6.6.1-5
./configure --prefix=/xxxx/upsuper.org/usr
make
make install
make clean
make distclean

至于 qrencode，也不过就是把里面的一些东西改动一下就可以了~由于 make 的优美特性，即使中间被中断也可以断点继续，很好很好~

最后编译好以后，再把我的 cgi 代码传上去，编译，通过了。把生成的文件拷贝出来加上 .cgi，却不能运行……

弄了个 bash 脚本上去加 ldd，发现 libqrencode.so 找不到，libmagick 看起来确实是有预装，在 /usr/lib 里面就轻易找到了。不甘心啊，这个东西应该是有办法解决的……

研究了一下，想起那个环境变量 LD_LIBRARY_PATH。我在 bash 的 cgi 里面加上

1	export LD_LIBRARY_PATH=/xxxx/upsuper.org/usr/lib

之后，ldd 里面所有库都可以找到了~真好

然后我就想啊，我们可不可以让他自动设置 LD_LIBRARY_PATH 这个环境变量呢？我就查了一下，发现在 .htaccess 里面加入 SetEnv 是可以的。

于是我就

1	SetEnv LD_LIBRARY_PATH /xxxx/upsuper.org/usr/lib

结果不行。弄了个 bash cgi：

#!/bin/bash
 
echo Content-Type: text/plain
echo
env

发现这个东西没有被加入环境变量当中。郁闷……在网上查了半天，没有发现任何人遇到过这种情况。我就想，不会吧，我这么倒霉？

然后我试了一下网上普遍喜欢用 SetEnv 的样例：

1	SetEnv TZ Asia/Chongqing

再打开刚才的查环境变量专用 cgi，发现出现了！

这就奇怪了。这说明 Apache 的 env 模块非常正常，但设置的东西却是选择性的进入 cgi。我猜一定是有什么白名单之类的东西。

于是我就跑去问服务商的技术客服……（又要用英语……）

话说这个主机商的技术客服看起来十分繁忙，我每次去都要等很久很久才有人应答。不知道是我的英语描述水平的问题，还是对方无法理解这种高深的问题，他直接告诉我这是我自己代码的问题，跟他没关系……然后居然就直接关断了聊天，我靠，这什么态度！我就囧了……又连了一次，我非常耐心的跟那个客服解释了一遍我想要的效果，他终于明白了！而且非常耐心的帮我解决问题。最后他告诉我，因为服务器使用了 suEXEC 这个模块，对进入 cgi 的环境变量做了一个白名单，所以不行……

唉，不行就不行呗，早知这样，我也不是没有办法的……

在 suexec 的代码里查到了环境变量白名单，发现所有 HTTP_ 和 SSL_ 开头的都是白名单范围的。于是我就写下了这样一个 .htaccess：

RewriteEngine On
 
RewriteCond %{REQUEST_URI} \.cgi$
RewriteRule ^(.*)$ /setenv.sh?cgi=$1 [QSA]
 
RewriteCond %{QUERY_STRING} ^cgi=([^&]+)
RewriteRule . - [E=HTTP_CGI:%1]

与之配套的 setenv.sh：

#!/bin/bash
 
export PATH=$PATH:/xxxx/upsuper.org/usr/bin
export CPATH=/xxxx/upsuper.org/usr/include
export LIBRARY_PATH=/xxxx/upsuper.org/usr/lib
export LD_LIBRARY_PATH=/xxxx/upsuper.org/usr/lib
 
./$HTTP_CGI

终于可以了……

这个论原理也很简单，就是把对 cgi 的请求全部转化到一个专门的设置环境变量的 bash 脚本里，然后让它调用相应的 cgi。对了，这个 setenv.sh 必须有执行权限才行……

这一切终于结束了，然后这个生成器也可以用了~

最后是一个小小的东西：

QR code 生成器

参考资料

Erlang 版 0-1 背包算法

upsuper — Sun, 04 Apr 2010 12:30:19 +0000

因为以后可能要用到，于是今天用 Erlang 写了一个 0-1 背包算法，这算是我第一次用 Erlang 写算法，也是我第二次写 Erlang 程序……

以前都使用正常语言写，各种算法自然不会难写。所谓正常的语言，呃……也就是我们平时用的最多的语言了，像 C++、Python 什么的在我看来都还算正常的语言。那么 Erlang 到底哪里不正常了呢？

其实 Erlang 只要两个不正常的地方就足够囧死人了……那就是：

没有循环，必须用递归代替
数组不能随机访问，只能从头部读写

其实说数组不能随机访问吧，应该也是能的，只不过时间恐怕就不是 O(1) 了罢了。

首先我们回忆一下一维 0-1 背包算法是怎么写的？

在这里就不分析这个算法了，不知道 0-1 背包算法的，呃……算导里竟然没有……那就自己 Google 好了……

先贴一个随手写的伪代码：

for i = 1 to n
    for j = size - s[i] to 0
        if b[j] + v[i] > b[j+s[i]]
            b[j+s[i]] = b[j] + v[i]

唉……要是 Erlang 能像这个伪码这么简单就好了……

不过考虑到 Erlang 的特性，最后没有使用伪码这种形式，主要的改动有两点，一是把连续的大小数组改成了可能达大的大小的列表，二是使用双列表最后归并。至于记录路径什么的，都是小事……

最后贴代码

-module(knapsack).
-export([calc/2]).
 
calc(Objs, Size) ->
    K = [{0, 0, []}],
    calc(Objs, K, Size).
 
calc([], K, _) -> final(K);
calc([Obj|Objs], K, Size) ->
    NewK = calc(Obj, K, Size, []),
    FinalK = merge(lists:reverse(K), NewK),
    calc(Objs, FinalK, Size).
 
final(K) -> lists:reverse(final(K, {0, []})).
final([], {_, NameList}) -> NameList;
final([{_, Val, List}|K], {MaxVal, _} = Cur) ->
    if
        Val > MaxVal -> final(K, {Val, List});
        true -> final(K, Cur)
    end.
 
calc(_, [], _, NewK) -> NewK;
calc({_, ObjSize, _}, [{CurSize, _, _}|_], Size, NewK)
        when CurSize + ObjSize > Size ->
    calc(nothing, [], 0, NewK);
calc({ObjName, ObjSize, ObjVal} = Obj,
        [{CurSize, CurVal, CurList}|K], Size, NewK) ->
    calc(Obj, K, Size,
        [{CurSize + ObjSize, CurVal + ObjVal, [ObjName|CurList]}|NewK]).
 
merge(A, B) -> merge(A, B, []).
merge([], [], K) -> K;
merge(A, [], K) -> lists:reverse(A, K);
merge([], B, K) -> lists:reverse(B, K);
merge([{SizeA, ValA, ListA} = TA|LeftA] = A,
      [{SizeB, ValB, ListB} = TB|LeftB] = B, K) ->
    if
        SizeA > SizeB -> merge(LeftA, B, [TA|K]);
        SizeA < SizeB -> merge(A, LeftB, [TB|K]);
        ValA > ValB -> merge(LeftA, LeftB, [TA|K]);
        ValA < ValB -> merge(LeftA, LeftB, [TB|K]);
        length(ListA) < length(ListB) ->
            merge(LeftA, LeftB, [TA|K]);
        true ->
            merge(LeftA, LeftB, [TB|K])
    end.

用法解释：执行 knapsack:calc(物件列表, 背包大小)，其中物件列表中每个物件以元组 {物件名称, 物件大小, 物件价值} 表示，最终返回物件名称的列表，价值相同情况下这个算法会返回包内物件平均价值最高的一种选法 (数量最少)。此外，这个算法保证输出列表中物件的顺序与输入顺序相同。

我知道代码风格很丑……而且还没有精心调试过……因为终于写出来太高兴了，所以就马上贴了，如果各位觉得不太对，请及时指出！谢谢……

人人网的好友交集小工具

upsuper — Sat, 09 Jan 2010 05:06:30 +0000

昨晚花了几个小时写了个计算人人网（其实我还是更喜欢叫他校内）当中，任意两个用户的好友之间交集的在线小工具，可以到我的实验室里看看这个小工具：人人网好友交集。

话说能写出这样的工具，主要有赖于我的空间提供商将系统换为 FreeBSD 后可以解析 Python 了，而且也没有限制 CGI，所以就成功了~

由于必须从人人网，而非本地，的用户页面，而非接口，获取数据，速度自然不可能快。所以我才有了那个提示：这一过程可能非常缓慢。不过事实上，经过测试，在我的网站的服务器上，速度还是相当了得的，基本上 2s-5s 就能出结果！想我自己的机子上好的时候都要 7s，不好的时候就根本出不来了……

如果单线程进行获取，那就真可以等死人了……所以用了5线程并发，每个线程维护一个链接……此外还连接的是手机人人网，以换来较小的传输流量，较快的解析速度，和较容易的处理方式……

最后，看代码之前声明一下，这个是以 AGPLv3 协议发布的，根据这一协议，如果你修改了这个代码并且使用你修改过的代码为他人服务，你所修改的代码也必须公开，并且以 AGPLv3 协议发布。

由于前台界面的东西大家都有办法获取，我就不给了，这里直接给出那个后台处理程序：

') class RequestThread(Thread): def__init__(self, id, sid, friends): Thread.__init__(self)self.__page = '/getfriends.do?curpage=%%d&id=%d&sid=%s'%(id, sid)self.__friends = friends self.__conn = HTTPConnection('3g.renren.com') def run(self): global curpage, stop_sign whilenot stop_sign: # 获取当前页面 curpage_lock.acquire() page = curpage curpage += 1 curpage_lock.release() # 连接获取数据 conn = self.__conn conn.request('GET', self.__page %(page, ), None, {'User-Agent': USER_AGENT }) data = conn.getresponse().read() # 处理获取的数据 friend_iter = parse_re.finditer(data) t_stop_sign = Truefor f in friend_iter: t_stop_sign = Falseid = int(f.group(1))self.__friends[id] = (f.group(3), f.group(2)) stop_sign = t_stop_sign def readFriends(id): # 初始化多线程global curpage, curpage_lock, stop_sign curpage = 0 curpage_lock = Lock() stop_sign = False threads = [] friends = {} # 创建线程for i inxrange(THREADS_NUM): threads.append(RequestThread(id, sid, friends))# 开始执行线程for t in threads: t.start()# 等待线程结束for t in threads: t.join() return friends def main(): # 初始化CGI输出print'Content-Type: text/plain'print # 获取ID信息 form = cgi.FieldStorage() id1, id2 = 0, 0if form.has_key('id1'): id1 = int(form['id1'].value)if form.has_key('id2'): id2 = int(form['id2'].value)ifnot(id1 and id2): json_dump({'error': True})return # 登入人人网 conn = HTTPConnection('3g.renren.com') conn.request('POST', '/login.do', urllib.urlencode({'email': RENREN_USER, 'password': RENREN_PWD}), {'Content-Type': 'application/x-www-form-urlencoded', 'User-Agent': USER_AGENT, }) response = conn.getresponse() data = response.read() conn.close() match = re.search(r'sid=([0-9a-f]+)', data, re.I)global sid sid = match.group(1) # 读取好友 friends1 = readFriends(id1) friends2 = readFriends(id2)iflen(friends1)>len(friends2): friends1, friends2 = friends2, friends1 # 判断交集 intersection = []for i in friends1.iterkeys(): if i in friends2: intersection.append((i, friends2[i][0], friends2[i][1])) json_dump({'error': False, 'count': len(intersection), 'data': intersection}) if __name__ == '__main__': try: main()except: json_dump({'error': True})

#!/usr/bin/env python
# - * - coding: UTF-8 - * -
 
# Copyright (C) 2010 Upsuper 
# License: AGPLv3
 
from threading import Thread, Lock
from httplib import HTTPConnection
import urllib
import cgi, json
import sys, re
 
# 常数设置
THREADS_NUM = 5             # 最大线程数
RENREN_USER = ''            # 人人网用户名
RENREN_PWD  = ''            # 人人网密码
USER_AGENT  = 'urenren 0.1' # 提交给人人网的 User-Agent
 
# JSON 输出函数
json_dump = lambda v: json.dump(v, sys.stdout)
# 编译匹配用正则表达式
parse_re = re.compile(r'

\?]+\?id=(\d+)[^"]*">' r'

' r'([^<]+)

代码中的配置常数 RENREN_USER 必须是一个已有的人人网账户的帐户名或邮箱地址，反正就是可以用来登入人人网的东西，RENREN_PWD 也必须是前面的账户对应的密码。因为人人网不允许未登入用户查看各种页面（大约是为了防止搜索引擎的人肉吧……不过其实我们有办法让搜索引擎索引到那些内容的~嘿嘿，突然有了一些邪恶的想法~），所以只能弄一个帐号来登入了。

人人网的手机版为了适应部分手机浏览器不支持 Cookies，用了 SID 的策略，这样我也就可以放心大胆地不管 Cookies 了~这个程序还是蛮好写的说~也只有 Python 能提供如此快速的程序构建~

方便使用 VC6 编译器的小脚本

upsuper — Wed, 30 Dec 2009 10:45:56 +0000

昨天的那篇日志跨编译器的 C 语言 NaN 支持当中讲到了如何在 Linux 下用 wine 执行 VC6 编译器编译程序，不过总觉得还是有那么些麻烦，要把待编译的文件复制到 VC6 的安装目录，还要写那么长一串东西。要是能像调用 GCC 那么方便就好了~

于是就有了下面这个小脚本：

#!/bin/bash
# - * - coding: UTF-8 - * -
 
VC6_DIR="这里写上VC6的安装地址"
 
BIN="$VC6_DIR/VC98/Bin"
export INCLUDE="$VC6_DIR/VC98/Include"
export LIB="$VC6_DIR/VC98/Lib" 
 
ARGS=
 
while getopts "o:c" optname
do
  case "$optname" in
  "o")
    ARGS="$ARGS /o$OPTARG"
    ;;
  "c")
    ARGS="$ARGS /c"
    ;;
  esac
done
 
wine "$BIN/CL.EXE" $ARGS ${@:$OPTIND}

然后把他放在 PATH 里面的某个目录下 (我放在了用户级的 /home/upsuper/bin 里，这个似乎要自己添加就是了)，然后给这个文件加上可执行属性，最后只要在需要的地方执行：

1	vc6 xxx.cpp

就解决了~

不过从这个脚本中也可以看出，它的功能还不太完善，不对，是很不完善。目前支持设置输出文件名和阻止执行连接。我很想加入很多其他的参数，不过不知道该怎么弄……

参考资料：

鼠标控制音乐播放的小程序

upsuper — Tue, 03 Nov 2009 12:06:03 +0000

看这个标题一定很奇怪：难道我以前控制音乐播放不用鼠标么？这个文章的标题看起来像鼠标刚刚发明的推广期的文章……不过，当然不是这样的！

其实只是想：如何把我的小本合上放书架上，当作一个音乐播放器+功放，并用我的无线6键鼠当遥控器遥控控制之。

想想其实还是蛮有意义的功能，这样我做作业的时候可以不用戴耳机，不用用MP3，直接把本当播放器；同时，我不会看到屏幕上的东西，可以安心做作业~再看看我的6键无限鼠，那额外的功能键平时根本不用，也想不出能有什么用……这么好的东西就这样被我浪费了……（话说，拿本当音乐播放器是不是更浪费？）

说干就干！

首先提出构想：左键用于暂停和播放，滚轮调节音量，侧边的两个功能键用来切换上一首和下一首。至于右键和中键……再说吧，说不定以后可以扩展更多功能？说不定以后高兴了弄个鼠标手势什么的~嘿嘿

接下来查找资料。印象中我的 Audacious 是可以用 D-Bus 控制的。简单地查阅了一下相关资料，发现了一个叫做 MPRIS 的播放器控制接口。为此，我还专门学习了一下 python-dbus 的使用。

插一句话：python-dbus 怎么没有中文教程啊！英文教程看得还是蛮吃力的……

连接播放器的代码很简单：

1
2
3

import dbus
bus = dbus.SessionBus()
player = bus.get_object('org.mpris.audacious', '/Player')

其中 org.mpris.audacious 指的是 audacious，支持 MPRIS 的其他播放器还可以有如下：

org.mpris.bmp
org.mpris.vlc
org.mpris.xmms2

至于那个 /Player 指的是我们要获取的是一个叫 Player 的对象，用于控制播放器。在 MPRIS 中除了 Player 以外，还有根对象 (‘/’) 和 TrackList 对象 (‘/TrackList’) 两个，具体的用法其文档里详细说明。这里我只用到了 Player 对象……

看看 Player 对象的方法还挺多，这里我用了这几个：

Next() – 切换到下一首歌曲
Prev() – 切换到上一首歌曲
Pause() – 暂停/继续
Play() – 播放歌曲
VolumeGet() – 获取当前音量
VolumeSet() – 设置音量

根据前面写的需求，这些已经足够了~

怎么用这些？很简单嘛！就像平时在 Python 里面调用类的方法一样~类型转换？不，那完全不需要你关心！

Python 程序拥有控制播放器的能力了，可是鼠标呢？如何捕获鼠标的动作呢？

想想……屏幕都关了……你还能指望他为你显示什么呢？干脆建立一个窗口把整个屏幕盖住算了！然后让他截获鼠标事件。

快速学习了一下 pygtk 的用法（其实就是照着 Hello World 打了一遍），差不多就知道怎么用了。窗口的建立代码大约如下：

window = gtk.Window(gtk.WINDOW_TOPLEVEL)
window.fullscreen()
window.show()
gtk.main()

这段代码很好理解了，第一行就是创建窗口对象，第二行令其全屏，第三行设置其显示，第四行开始 GTK 的主循环。

下面就是鼠标事件的问题了~添加事件用的是类似下面的代码：

1 2	window.connect('key-press-event', key_press) window.connect('destroy', lambda widget, data = None: gtk.main_quit())

第一个是注册了键盘按键的事件捕获，将 key-press-event 事件绑定上 key_press 这个函数，而把窗口的销毁事件 destroy 连接到退出 GTK 主循环。说一句废话，连接事件的语句要放在 gtk.main() 之前……key_press 主要用来实现当点击 escape 时退出：

1
2
3

def key_press(widget, event, data = None):
  if gdk.keyval_name(event.keyval) == 'Escape':
    window.destroy()

好了，下面是鼠标点击事件 button-press-event：

1
2
3

def button_press(widget, event, data = None):
  print event
window.connect('button-press-event', button_press)

试试？不行的！差了一番后，发现 Window 这个东西默认是不打开鼠标点击的事件捕获的，逼近大家很少会点窗口本身，都是点里面的按钮什么的……下面的代码要求窗口打开鼠标点击事件的捕获：

1	window.add_events(gtk.gdk.BUTTON_PRESS_MASK)

接下来就可以了~event.button 表示的是点击的键是哪个。检查了一下，左键是1，右键是3，中键是2，功能键1是8，功能键2是9。滚轮呢？

原来滚轮有自己的事件 scroll。

1
2
3

def scroll(widget, event, data = None):
  print event
window.connect('scroll-event', scroll)

直接可以用了。识别滚轮方向的是 event.direction，应该是 gtk.gdk.SCROLL_UP、SCROLL_DOWN、SCROLL_LEFT、SCROLL_RIGHT 其中的一个。

下面根据最初的设想，把他们都拼起来就形成了最初的版本：

#!/usr/bin/env python
 
import pygtk
pygtk.require('2.0')
import gtk
from gtk import gdk
 
import dbus
 
def button_press(widget, event, data = None):
  button = event.button
  if button == 1:
    if player.GetStatus() == 2:
      player.Play()
    else:
      player.Pause()
  elif button == 9:
    player.Prev()
  elif button == 8:
    player.Next()
 
def scroll(widget, event, data = None):
  if event.direction == gdk.SCROLL_UP:
    player.VolumeSet(player.VolumeGet() + 3)
  elif event.direction == gdk.SCROLL_DOWN:
    player.VolumeSet(player.VolumeGet() - 3)
 
def key_press(widget, event, data = None):
  if gdk.keyval_name(event.keyval) == 'Escape':
    window.destroy()
 
bus = dbus.SessionBus()
player = bus.get_object('org.mpris.audacious', '/Player')
 
window = gtk.Window(gtk.WINDOW_TOPLEVEL)
window.connect('destroy', lambda widget, data = None: gtk.main_quit())
window.add_events(gdk.BUTTON_PRESS_MASK)
window.connect('button-press-event', button_press)
window.connect('key-press-event', key_press)
window.connect('scroll-event', scroll)
window.fullscreen()
window.show()
gtk.main()

OK 看起来不错了，也能控制了。可是，作为一个程序员，怎么能不考虑各种意外呢？试想，如果我们控制着控制着，突然，一个不明真相的窗口弹了出来怎么办呢？！这是一个问题……

最容易想到的方法：把这个窗口永久置顶！

事实上也确实找到了这么个方法，就是把窗体注册成 dock。我们知道 dock 是需要永久置顶的~

用如下方法可以让窗体变成 Dock：

1	window.set_type_hint(gtk.gdk.WINDOW_TYPE_HINT_DOCK)

但是问题又来了：dock 似乎没法接受键盘输入。因为 dock 只是形式上覆盖，并且可以拦截一切在他上面发生的鼠标事件，可键盘焦点就未必在他那里了~

那么我们就不要键盘事件了吧~于是我就把键盘退出给删了，改成了双击退出~查了下资料，双击事件也通过 button-press-event 发送的，必须根据 event.type 来判断是单击、双击还是三击。此外，双击和三击中的每一次点击都会引发一次单击事件。

OK，于是我们把 button_press 函数改成下面这样：

def button_press(widget, event, data = None):
  if event.type == gdk._2BUTTON_PRESS:
    window.destroy()
    return
  elif event.type == gdk._3BUTTON_PRESS:
    return
  button = event.button
  if button == 1:
    if player.GetStatus() == 2:
      player.Play()
    else:
      player.Pause()
  elif button == 9:
    player.Prev()
  elif button == 8:
    player.Next()

然后删掉键盘处理，最后就成了这样：

#!/usr/bin/env python
 
import pygtk
pygtk.require('2.0')
import gtk
from gtk import gdk
 
import dbus
 
def button_press(widget, event, data = None):
  if event.type == gdk._2BUTTON_PRESS:
    window.destroy()
    return
  elif event.type == gdk._3BUTTON_PRESS:
    return
  button = event.button
  if button == 1:
    if player.GetStatus() == 2:
      player.Play()
    else:
      player.Pause()
  elif button == 9:
    player.Prev()
  elif button == 8:
    player.Next()
 
def scroll(widget, event, data = None):
  if event.direction == gdk.SCROLL_UP:
    player.VolumeSet(player.VolumeGet() + 3)
  elif event.direction == gdk.SCROLL_DOWN:
    player.VolumeSet(player.VolumeGet() - 3)
 
bus = dbus.SessionBus()
player = bus.get_object('org.mpris.audacious', '/Player')
 
window = gtk.Window(gtk.WINDOW_TOPLEVEL)
window.connect('destroy', lambda widget, data = None: gtk.main_quit())
window.add_events(gdk.BUTTON_PRESS_MASK)
window.connect('button-press-event', button_press)
window.connect('scroll-event', scroll)
window.set_type_hint(gdk.WINDOW_TYPE_HINT_DOCK)
window.fullscreen()
window.show()
gtk.main()

完成了！

以后就可以拿我的小本来当音乐播放器+功放啦~

参考资料：

Python 物理实验数据处理类

upsuper — Tue, 20 Oct 2009 09:06:09 +0000

非常无语的，同济大一什么物理理论都没学，需要的高数基础也都没有学过，竟然要我们先学物理实验，而且要处理数据，计算平均值和标准误差什么的，十分麻烦……

上次处理数据弄了个零碎的小程序来，很麻烦，这次就把他统一编为一个模块，方便以后使用……这个模块里面的说有公式依据是同济大学出版社《物理实验教程》第1版第二章的内容。

如何使用这个模块呢？首先，作为一个 Python 模块，肯定要用 Python……

然后我用我第一次做实验的部分数据演示一下：

直接输入数据，标准格式输出~

最后是代码：

#!/usr/bin/env python
# -*- coding: utf-8 -*-
 
import math
sqrt = math.sqrt
 
def gen_imet(func):
  def retfunc(self, other):
    ret = func(self, other)
    self.ave = ret.ave
    self.u = ret.u
    self.digs = ret.digs
    return self
  return retfunc
 
def get_digs(a, b):
  if a == 0: return b
  elif b == 0: return a
  else: return min(a, b)
 
class PhyStat:
  """
  物理实验数据
 
  PhyStat(ave, u = 0.0, digs = 0)
  ave   表示平均值，如果传入一列数则自动计算平均值及其标准误差
  u     表示标准误差，如果 ave 传入数据列表则此处表示仪器精度
  digs  表示有效数字个数，0 表示不限制
  """
 
  def __init__(self, ave, u = 0.0, digs = 0):
    """
    对于 ave 为数据列表的情况，使用下面公式计算平均值：
    ave = Σx[i] / n
    根据下面公式计算标准误差：
    S = Sqrt(Σ(x[i]-ave)^2 / (n - 1))
    u = Δinst / Sqrt(3)
    U = Sqrt(S^2 + u^2)
    """
    self.digs = digs
 
    if isinstance(ave, tuple) or isinstance(ave, list):
      data = ave
      n = len(data)
      ave = math.fsum(ave) / n
      self.ave = ave
 
      s = 0.0
      for i in data:
        s += (i - ave) ** 2
      s = sqrt(s / (n - 1))
      u /= sqrt(3)
      self.u = sqrt(s ** 2 + u ** 2)
    else:
      self.ave = float(ave)
      self.u = abs(float(u))
 
  def __repr__(self):
    ave = self.ave
    exp = int(math.floor(math.log10(ave)))
    exp10 = 10 ** exp
    ave = ave / exp10
    u = self.u / exp10
    digs = self.digs
    if digs == 0:
      ret = '%f ±%f'
    else:
      ret = '%%.%df ±%%.%df' % (digs, digs)
    if exp:
      ret = '(%s) E %d' % (ret, exp)
    return ret % (ave, u)
 
  __neg__ = lambda self: PhyStat(-self.ave, self.u, self.digs)
  __pos__ = lambda self: self
  __abs__ = lambda self: PhyStat(abs(self.ave), self.u, self.digs)
 
  def __add__(self, other):
    if isinstance(other, PhyStat):
      ave = self.ave + other.ave
      u = sqrt(self.u ** 2 + other.u ** 2)
      digs = get_digs(self.digs, other.digs)
      return PhyStat(ave, u, digs)
    else:
      ave = self.ave + other
      return PhyStat(ave, self.u, self.digs)
 
  __radd__ = __add__
  __sub__ = lambda self, other: self + (-other)
  __rsub__ = lambda self, other: other + (-self)
 
  def __mul__(self, other):
    if isinstance(other, PhyStat):
      a = self.ave
      b = other.ave
      ave = a * b
      u = ave * sqrt((self.u / a) ** 2 + (other.u / b) ** 2)
      digs = get_digs(self.digs, other.digs)
      return PhyStat(ave, u, digs)
    else:
      ave = self.ave * other
      u = self.u * other
      return PhyStat(ave, u, self.digs)
 
  __rmul__ = __mul__
 
  def __div__(self, other):
    if isinstance(other, PhyStat):
      return self * other ** -1
    else:
      return (1.0 / other) * self
 
  __rdiv__ = lambda self, other: other * self ** -1
 
  def __pow__(self, other):
    a = self.ave
    ave = a ** other
    u = abs(other * self.u / a * ave)
    return PhyStat(ave, u, self.digs)
 
  __iadd__ = gen_imet(lambda self, other: self + other)
  __isub__ = gen_imet(lambda self, other: self - other)
  __imul__ = gen_imet(lambda self, other: self * other)
  __idiv__ = gen_imet(lambda self, other: self / other)
  __ipow__ = gen_imet(lambda self, other: self ** other)
 
  __int__ = lambda self: int(self.ave)
  __long__ = lambda self: long(self.ave)
  __float__ = lambda self: self.ave

如果大家发现什么 bug，或者有什么改进建议，欢迎提哈~

ujudger 0.1 发布

upsuper — Fri, 10 Jul 2009 02:54:48 +0000

曾经的曾经，很早很早以前，我发布了一个叫 ujudger-base 0.1 的东西，是用 C++ 写成的简单测评程序，仅能测试程序运行时间和内存耗用（还不准确）。

这次发布的是一个完整的测评机，支持传统题型、交互题和提交答案题，可导出测评结果。整个程序使用 Python 编写，直接把 ujudger.py 复制到待测目录下并编写 config 文件即可测评。

不过现在这个程序只有一个文字界面，非常简陋，而且安全性不佳。但现在正在撰写的 ujudger 1.0 (版本跨越大，意味着变化之巨大……) 将就这些问题有重大改进，敬请期待。

另外，有人询问是否准备让这个测评机与另一个 Windows 下的开源测评机 Cena 兼容。我的回答是，由于设计差异，可能比较困难。对比 Cena 的自定义校验器，ujudger 的自定义校验器更强大，同时 Cena 的校验器可以通过一点简单的修改或加一个脚本外壳就可以在 ujudger 中使用。在 1.0 中计划兼容 Cena 的校验器。至于配置文件，1.0 中配置文件的格式会有较大变动，不过可能可以通过转换程序将 Cena 的配置文件转为 ujudger 的。

我在 Google Code 上开了一个项目，感兴趣的同学可以关注关注~地址：http://ujudger.googlecode.com/

最后，ujudger 0.1 的下载地址也在 Google Code 上：ujudger.py

Ubuntu 下直接挂载 iso

upsuper — Tue, 17 Feb 2009 03:11:50 +0000

感觉挂载 iso 每次要建文件夹、挂载，用完还要删除、删文件夹，很麻烦。在网上找了一个帖子：用于使用鼠标右键挂载iso的脚本!，觉得看起来不错，不过好像不能自动卸载。经过一番学习，给出了如下自动挂载/卸载脚本：

#!/bin/bash
[ $UID == 0 ] || cmd="gksu"
if [ ! -f $1 ]; then
	exit 2
fi
dir="/media/`basename $1 .iso`"
if [ ! -d $dir ]; then
	$cmd "mkdir $dir"
	$cmd "mount -t iso9660 -o loop $1 $dir"
else
	$cmd "umount $dir"
	$cmd "rm -r $dir"
fi

写的很短很简陋，如果把 cmd= 后面改成 kdesudo 就可以在 KDE 下用了。

我把这个设置为 iso 文件的默认打开方式，这样以后双击 iso 文件，马上挂载，再双击，自动卸载，太方便了！用 Linux 的都可以试试~

WAP版百度投票

upsuper — Mon, 16 Feb 2009 02:00:31 +0000

百度贴吧的投票不支持手机，连 UCWEB 和 Opera Mini 都看不了，于是就自己写了个，不过只能查看，不能投票就是了……

其实我一直很奇怪，连我当时的音乐投票站点用的全 AJAX 那俩浏览器都能完全支持，为什么百度这个却不行呢？经过研究，我在投票的地方发现如下代码：

1	document.write('+r_url+'&t='+Math.random()+'">')

这段代码出现在投票信息出现的地方，很明显，这就是投票信息。所以就打开 http://toupiao.baidu.com/sys/dtview/1/299ffee9f95d723176f722e6 看到了那个投票信息，再研究代码，便写出了下面程序：

header('Content-Type: text/vnd.wap.wml;charset=UTF-8');
 
$id = $_GET['id'];
$id = preg_replace('/[^0-9a-f]/i', '', $id);
 
if (! $id) exit('错误的投票 id ！');
 
$file = './cache/'.$id;
 
if (file_exists($file)) {
	include $file;
 
	$curl = curl_init();
	curl_setopt($curl, CURLOPT_URL, 'http://toupiao.baidu.com/sys/rsview/1/'.$id.'?t='.rand());
	curl_setopt($curl, CURLOPT_HEADER, 0);
	curl_setopt($curl, CURLOPT_RETURNTRANSFER, 1);
	$content = curl_exec($curl);
	curl_close($curl);
 
	preg_match('/vtVoterSum:(\d+),/i', $content, $match);
	if (! $match) exit('错误！');
	$total = $match[1];
 
	preg_match_all('/{vtVoteNum:(\d+),vtVotePt:(\d+\.\d+)}/', $content, $match, PREG_SET_ORDER);
	if (! $match) exit('错误！');
	for ($i = 0; isset($match[$i]); ++$i) {
		$item[$i][1] = $match[$i][1];
		$item[$i][2] = $match[$i][2].'%';
	}
 
} else {
	$curl = curl_init();
	curl_setopt($curl, CURLOPT_URL, 'http://toupiao.baidu.com/sys/dtview/1/'.$id.'?t='.rand());
	curl_setopt($curl, CURLOPT_HEADER, 0);
	curl_setopt($curl, CURLOPT_RETURNTRANSFER, 1);
	$content = iconv('gb2312', 'utf-8', curl_exec($curl));
	curl_close($curl);
 
	preg_match_all('/vote_desc_array\[\d+\]=new Array\("(.+?)",\d+\);/i', $content, $match, PREG_SET_ORDER);
	if (! $match) exit('错误！');
	$fp = fopen($file, 'wb');
	fwrite($fp, "\n");
	for ($i = 0; isset($match[$i]); ++$i) {
		$item[$i][0] = $match[$i][1];
		fprintf($fp, '$item[%d][0]=%s;', $i, var_export($item[$i][0], true));
	}
	fwrite($fp, "\n?>");
	fclose($fp);
 
	preg_match('/(\d+)<\/span>/i', $content, $match);
	//if (! $match) exit('错误！');
	$total = $match[1];
 
	preg_match_all('/{optVoNum:(\d+),vtVotePt:(\d+\.\d+)}/i', $content, $match, PREG_SET_ORDER);
	if (! $match) header('Location: index.php?id='.$id);
	for ($i = 0; isset($match[$i]); ++$i) {
		$item[$i][1] = $match[$i][1];
		$item[$i][2] = $match[$i][2].'%';
	}
}
 
include './userface.php';

至于 userface.php 就自己写 wml 界面就好了~

我的 WAP 上传通道

upsuper — Fri, 31 Oct 2008 10:07:01 +0000

这年头……手机连接线也要大几十块……我穷人买不起……于是就想利用我的空间建一个 WAP 上传通道，记得上一次放了一个下载通道。这样我要让手机和电脑互传就没有障碍了！

去网上找了一下，很短的时间写了个小程序：


// 前面就是 xhtml 的需要了
header("Content-Type: application/xhtml+xml");
echo "";
?>
DOCTYPE html
	PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
	"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="en" lang="en">
<head>
<meta http-equiv="content-type" content="application/xhtml+xml;charset=UTF-8" />
<title>WAP上传通道title>
head>
<body>

if ($_FILES['file'])
{
	// 验证密码，不想让谁都能随便上传……
	if (md5($_POST['password']) != '********************************') {
		echo "密码错误！";
	}
	else {
		define('UPLOAD_DIR', dirname(__FILE__) . '/upload');
		@mkdir(UPLOAD_DIR);
		if (move_uploaded_file($_FILES['file']['tmp_name'], 
			UPLOAD_DIR . '/' . basename($_FILES['file']['name']))) {
			echo basename($_FILES['file']['name']) . 
				" ({$_FILES['file']['size']}bytes) 上传成功！";
		} else {
			echo "上传失败！({$_FILES['file']['error']})
";
			switch ($_FILES['file']['error']) {
				case UPLOAD_ERR_INI_SIZE:
				case UPLOAD_ERR_FORM_SIZE:
					echo "大小超过限制的值"; break;
				case UPLOAD_ERR_PARTIAL:
				case UPLOAD_ERR_NO_FILE:
					echo "上传文件不完整"; break;
				case UPLOAD_ERR_NO_TMP_DIR:
				case UPLOAD_ERR_CANT_WRITE:
					echo "其他错误"; break;
			}
		}
	}
}
echo "
";
?>
<form name="upload" method="post" action="index.php" enctype="multipart/form-data">
上传密码：<input type="password" name="password" /><br />
<input type="file" name="file" /><br />
<input type="submit" value="提交" />
form>
body> 
html>

马上试验了一下，就成功了！

很好很强大！拿来分享一下~

我的 PHP 测试小程序

upsuper — Thu, 02 Oct 2008 01:52:11 +0000

最近研究 PHP，为了方便，写了个小小的 PHP 页面用来测试简单的 PHP 代码段，为了看了更方便。虽然可能有更好的工具，不过也胡乱贴出来，呵呵~

首先要先在这个文件的文件夹下建立一个 script_test.php 文件，还要把这个文件设为 0777（这是 Linux 的属性，用 Windows 的不用鸟……Windows 的权限设置过于宽松了……）。

然后是代码：


$script = stripslashes($_REQUEST['script']);
if ($script !== false) {
	$ftest = fopen("script_test.php", "wb");
	fwrite($ftest, $script);
	fclose($ftest);
}
?>
 if ($_REQUEST['ret'] == 'js') : ?>

 else : ?>
<style>textarea { width: 100%; height: 400px; }style>
<table width="100%">
<tr valign="top"><td width="50%">
<form action="?ret=js" method="POST" target="submit_iframe">
<textarea id="script" name="script">
 if ($script !== false) : ?> $script ?>php endif; ?>
textarea><br />
<input type="submit" value="Test" />
form>td>
<td width="50%">php if ($script !== false) : ?>
<iframe id="show_test" src="script_test.php"
 style="border: 1px black solid; width: 100%; height: 400px;">iframe>
 endif; ?>
td>tr>table>
<iframe id="submit_iframe" name="submit_iframe" style="display: none;">iframe>
 endif; ?>

其中那些高度对于我的显示器效果是不错的，可以根据自己需要调整……

真彩色展示

upsuper — Fri, 15 Feb 2008 01:21:34 +0000

比较无聊的我用VB写了一个能生成所有24位真彩色的程序，可以循环播放真彩色。在这个程序中，你开始会觉得看不出来有变化，不知不觉间就完全变掉了。每轮一共512桢，每桢256*256，向右绿色增加，向下蓝色增加，播放为红色的变化。

如果感兴趣可以在这里下载：allTrueColor.exe (28KB)

Windows 9x需要有VB6运行库支持。另外运行本程序可能导致CPU占用率极高（我用双核，其中一个核可以达到100%）。这个程序的变化速度取决于电脑的运行速度，但通常是比较快速的，这完全得益于使用Windows API代替VB自带的函数带来的效率。

下面是源代码

Option Explicit
Dim isBegin As Boolean, isAdd As Boolean
Dim i As Integer, j As Integer, k As Integer
Private Declare Function SetPixel Lib "gdi32" _
    (ByVal hdc As Long, _
    ByVal x As Long, ByVal y As Long, _
    ByVal crColor As Long) As Long
 
Private Sub cmdStart_Click()
    Dim hdc As Long
    hdc = pic.hdc
    isBegin = True
    Do
        For j = 0 To 255
            For k = 0 To 255
                SetPixel hdc, j, k, RGB(i, j, k)
            Next k
            DoEvents
        Next j
        If Not isBegin Then Exit Sub
        'SavePicture pic.Image, _
            'App.Path & "\" & CStr(i) & ".bmp"
        pic.Refresh
        DoEvents
        If isAdd Then
            i = i + 1
            If i = 255 Then isAdd = False
        Else
            i = i - 1
            If i = 0 Then isAdd = True
        End If
    Loop
End Sub
 
Private Sub cmdStop_Click()
    isBegin = False
End Sub
 
Private Sub Form_Load()
    isAdd = True
End Sub

整个窗体只有3个控件：一个PictureBox控件pic，长宽256px，记得要把ScaleMode设为Pixel；剩下的就是两个按钮cmdStart和cmdStop了

自己写的QuickSearch

upsuper — Thu, 14 Feb 2008 05:53:13 +0000

今天研究Discuz!代码的时候要搜索一个函数在哪里被调用，结果不懂怎么办。

Windows的搜索貌似某个服务被我禁了一直以来都不能用，而Google的Desktop又不收编这些东西……

花了一点时间写了一个QuickSearch，代码也很短，效率也很高。用法就是“QuickSearch 关键字”。就会遍历当前目录及其子目录下的所有的文件，查找包含关键字的文件。

如果你也喜欢可以在这里下载：QuickSearch.exe (16KB)

由于使用VS2005编译，需要.NET Framework 2.0支持。源代码在下方，可以用VS2003编译。

蛮好用的说，代码也很短：

using System;
using System.IO;
 
namespace QuickSearch
{
    class Program
    {
        static string curdir = Directory.GetCurrentDirectory();
        static int lcurdir = curdir.Length + 1;
        static string searchText;
 
        static void Main(string[] args)
        {
            searchText = args[0];
            SearchDir(curdir);
        }
 
        static string filecontent;
 
        static void SearchDir(string path)
        {
            string[] files = Directory.GetFiles(path);
            foreach (string file in files)
            {
                filecontent = File.ReadAllText(file);
                if (filecontent.IndexOf(searchText) > 0)
                    Console.WriteLine(file.Substring(lcurdir));
            }
            string[] dirs = Directory.GetDirectories(path);
            foreach (string dir in dirs)
                SearchDir(dir);
        }
    }
}

blackflip，Black Flip

upsuper — Wed, 14 Nov 2007 08:19:50 +0000

有一个很好玩的游戏，是Flash+AJAX的：http://blackflip.org/

玩法就是一笔划，然后划过的格子都会被反色，我们的目标是，使得每一行的颜色相同~

于是，经过艰辛的摸索，我编写了一个很好的工具blackflip_helper。

这个工具的原理就是每一行反色其实走法是一样的！只不过这样更容易看出来。

经过思考，发现这个游戏似乎使用搜索是非常麻烦的，而基于状态压缩的动态规划似乎也快不了多少，都是次方级的。所以这个游戏的解法有待继续研究哈。