mathematics and me

We must know, We will know。 -----Hilbert

2014年5月14日星期三

Project Euler答题(不定时更新)

继续答题。

2014.5.14

第11题(python):

grid = """08 02 22 97 38 15 00 40 00 75 04 05 07 78 52 12 50 77 91 08
          49 49 99 40 17 81 18 57 60 87 17 40 98 43 69 48 04 56 62 00
          81 49 31 73 55 79 14 29 93 71 40 67 53 88 30 03 49 13 36 65
          52 70 95 23 04 60 11 42 69 24 68 56 01 32 56 71 37 02 36 91
          22 31 16 71 51 67 63 89 41 92 36 54 22 40 40 28 66 33 13 80
          24 47 32 60 99 03 45 02 44 75 33 53 78 36 84 20 35 17 12 50
          32 98 81 28 64 23 67 10 26 38 40 67 59 54 70 66 18 38 64 70
          67 26 20 68 02 62 12 20 95 63 94 39 63 08 40 91 66 49 94 21
          24 55 58 05 66 73 99 26 97 17 78 78 96 83 14 88 34 89 63 72
          21 36 23 09 75 00 76 44 20 45 35 14 00 61 33 97 34 31 33 95
          78 17 53 28 22 75 31 67 15 94 03 80 04 62 16 14 09 53 56 92
          16 39 05 42 96 35 31 47 55 58 88 24 00 17 54 24 36 29 85 57
          86 56 00 48 35 71 89 07 05 44 44 37 44 60 21 58 51 54 17 58
          19 80 81 68 05 94 47 69 28 73 92 13 86 52 17 77 04 89 55 40
          04 52 08 83 97 35 99 16 07 97 57 32 16 26 26 79 33 27 98 66
          88 36 68 87 57 62 20 72 03 46 33 67 46 55 12 32 63 93 53 69
          04 42 16 73 38 25 39 11 24 94 72 18 08 46 29 32 40 62 76 36
          20 69 36 41 72 30 23 88 34 62 99 69 82 67 59 85 74 04 36 16
          20 73 35 29 78 31 90 01 74 31 49 71 48 86 81 16 23 57 05 54
          01 70 54 71 83 51 54 69 16 92 33 48 61 43 52 01 89 19 67 48"""

grid = grid.replace('\n', '').replace(' ', '')
# creat a matrix
matrix = [[0 for i in range(20)] for i in range(20)]
for i in range(20):
    for j in range(20):
        matrix[i][j] = int(grid[2*(20*i + j):2*(20*i + j) + 2])
        #print matrix[i][j]

# find the largest of rows
lor = 0 # the largest of rows
for i in range(20):
    for j in range(17):
        if matrix[i][j] * matrix[i][j + 1] * matrix[i][j + 2] * matrix[i][j + 3] > lor:
            lor = matrix[i][j] * matrix[i][j + 1] * matrix[i][j + 2] * matrix[i][j + 3]

# find the largest of cols
loc = 0 # the largest of cols
for i in range(17):
    for j in range(20):
        if matrix[i][j] * matrix[i + 1][j] * matrix[i + 2][j] * matrix[i + 3][j] > loc:
            loc = matrix[i][j] * matrix[i + 1][j] * matrix[i + 2][j] * matrix[i + 3][j]

# find the largest of diagonal lines
lod = 0 # the largest of diagonal lines
for i in range(17):
    for j in range(17):
        if matrix[i][j] * matrix[i + 1][j + 1] * matrix[i + 2][j + 2] * matrix[i + 3][j + 3] > lod:
            lod = matrix[i][j] * matrix[i + 1][j + 1] * matrix[i + 2][j + 2] * matrix[i + 3][j + 3]
# find the largest of anti diagonal lines
load = 0 # the largest of anti diagonal lines
for i in range(17):
    for j in range(3, 20):
        if matrix[i][j] * matrix[i + 1][j - 1] * matrix[i + 2][j - 2] * matrix[i + 3][j - 3] > load:
            load = matrix[i][j] * matrix[i + 1][j - 1] * matrix[i + 2][j - 2] * matrix[i + 3][j - 3]

print max(lor, loc, lod, load)

2014年5月8日星期四

Project Euler答题(6--10)

为防止帖子太长,故新开一贴,继续答题。
2014.5.8

第六题(python):

sum_of_squares = 0
sum_of_number = 0

for i in range(1, 101):
    sum_of_squares = sum_of_squares + i*i

for i in range(1, 101):
    sum_of_number = sum_of_number + i

print sum_of_number**2 - sum_of_squares

2014.5.8

第七题(python):

import math
import random
import sys   
sys.setrecursionlimit(1000000)

# def is_prime(num):
#     i = 2
#     isprime = True
#     while i <= int(math.sqrt(num)):
#         if num % i == 0:
#             isprime = False
#             break
#         i = i + 1
#     return isprime

# def is_prime(num):
#     if num % 2 == 0:
#         return False;
#     i = 3
#     isprime = True
#     while i <= int(math.sqrt(num)):
#         if num % i == 0:
#             isprime =  False
#             break;
#         i = i + 2
#     return isprime

# prime_list = [2, 3]
# def is_prime(num):
#     isprime = True
#     i = 0
#     while i < len(prime_list):
#         if num % prime_list[i] == 0:
#             isprime = False
#             break
#         i = i + 1
#     return isprime

prime_list = [2, 3]
def is_prime(num):
    isprime = True
    i = 0
    while i < len(prime_list) and prime_list[i] <= int(math.sqrt(num)):
        if num % prime_list[i] == 0:
            isprime = False
            break
        i = i + 1
    return isprime

# def is_prime(num):
#     f = 5
#     isprime = True
#     r = int(math.sqrt(num)) + 1
#     if num % 2 == 0 or num % 3 == 0:
#         return False
#     while f <= r:
#         if num % f == 0:
#             isprime = False
#         elif num % (f + 2) == 0:
#             isprime = False
#         f = f + 6
#     return isprime

# Wilson' theroem
# def factorial(n):
#     if (n == 1):
#         return n
#     else:
#         return n * factorial(n - 1)
# def is_prime(num):
#     if (factorial(num - 1) + 1) % num == 0:
#         return True
#     else:
#         return False

# Fermat's little theroem
# def FermatPrimalityTest(number):
#     for time in range(10):
#         randomNumber = random.randint(2, number - 1)
#         if ( pow(randomNumber, number-1, number) != 1 ):
#             return False
#     return True


number_of_prime = 2
num = 5
while True:
    if is_prime(num):
        prime_list.append(num)
        number_of_prime = number_of_prime + 1
        if number_of_prime == 40001:
            print num
            break
    # if is_prime(num) == False and FermatPrimalityTest(num) == True:
    #     print num
    num = num + 2

2014年5月7日星期三

复旦大学android客户端开发计划

        最近android的知识也学得差不多了,遂决定做个小应用玩一玩。转头一看发现市场上貌似什么应用都有了,不需要我再去做什么。于是想到了要做一个集成我旦一些服务的应用,这些服务包括但不仅限于复旦邮箱(重点打造),复旦图书馆,e-learning,urp服务等。未来也许还会添加其他功能。
        为什么要做这个东西?其实一开始我只是想把邮箱功能做好。虽然市场上已经有不少邮箱的app了,但是这些应用要么界面粗糙,要么不支持复旦邮箱,最重要的是它们都没有为复旦邮箱量身订做的功能!!!为什么呢?因为复旦邮箱有复旦“特色”,经常充斥着各种“求填问卷”,“社团广告”等。因为我的手机邮箱里设置了提醒功能,以便随时能接受重要邮件。一有新邮件,就会自动提醒。每天被这些“垃圾邮件”骚扰五遍以上,你知道我的感受吗!!!扫邮什么的最烦了,必须把它屏蔽掉!!!最近我突然发现了这些“问卷贴”的一个共同特点(虽然很明显),利用这个特点,我就可以轻松地把这些邮件屏蔽掉!!!这个特(ruo)点(dian)是什么暂时先保留,以免“道高一尺魔高一丈”。
        另一个想做好的功能是图书馆的各项服务。从此只要手机在手,你就可以轻松查询,预约,续借图书。虽说浏览器也能做到,但毕竟不适合在手机端浏览。
        目前基本框架已搭好,剩下的就是把各项功能添加进去了。我是极为注重用户体验的,所以花了很多时间来做UI。没有一个好的界面,用起来怎么可能爽!虽然我不是一个专业的UI/UX designer,但我会尽我所能把UI做好。我们并不缺技术(technology),缺少的是人文(liberal arts)关怀。所以我一直都很佩服apple,佩服乔帮主,apple是少数能把技术和人文光怀都做到极致的公司。我也一直很喜欢apple的各项产品,虽然我目前没有经济能力去拥有它。
科技与人文的十字路口
        目前项目托管在Bitbucket,项目地址:https://bitbucket.org/egrcc/fudanclient ,由于项目尚在开发中,所以不会放到github上去,等项目完成后会迁移到github上来。由于最近事情比较多,再过一段时间又要到学期末了,估计开发会很缓慢,希望在暑假我能把这个项目做完。