题目id: 1 just print a+b
give you two var a and b, print the value of a+b, just do it!!
print a+b
题目id: 2 list排序
给你一个list L, 如 L=[2,8,3,50], 对L进行升序排序并输出
print sorted(L)
题目id： 3 字符串逆序
给你一个字符串 a， 如a=‘12345’，对a进行逆序输出a。
print a[::-1]
题目id：4 输出字典key
给你一字典a，如a={1:1,2:2,3:3}，输出字典a的key，以','链接，如‘1,2,3'。
print ','.join(a.keys())
题目id: 5 输出字符奇数位置的字符串
给你一个字符串 a， 输出字符奇数位置的字符串。如a=‘12345’，则输出135。
str1=''
for i in range(len(a)):
if i%2==0:
str1 += str(a[i])
print str1
题目id：6 求解100以内的所有素数
输出100以内的所有素数，素数之间以一个空格区分
import math
def isPrime(n):
if n<=1:
return False
for i in range(2,int(math.sqrt(n))+1):
if n%i == 0:
return False
return True
list = []
for i in range(1,100):
if isPrime(i):
list.append(str(i))
print ' '.join(list)
题目id：7 求矩形面积
已知矩形长a,宽b,输出其面积和周长，面积和周长以一个空格隔开
print "%d %d"%(a*b,2*(a+b))
题目id：8 求中位数
给你一个list L, 如 L=[0,1,2,3,4], 输出L的中位数（若结果为小数，则保留一位小数）。
n = len(L)
L.sort() //注意sorted()函数和sort()函数的区别，sorted()不影响L的结构，而sort()影响L的结构
print (L[n/2] if n&1 else round(float(L[n/2-1]+L[n/2])/2,1))
题目id：9 最大公约数
给你两个正整数a和b， 输出它们的最大公约数。
def gcd(a,b):
if a < b:
a,b = b,a
while b:
a,b=b,a%b
return a
print gcd(a,b)
题目id：10 最小公倍数
给你两个正整数a和b， 输出它们的最小公倍数
def gcd(a,b):
if a < b:
a,b = b,a
while b:
a,b=b,a%b
return a
print a*b/gcd(a,b)
题目id：11 结尾0的个数
给你一个正整数列表 L, 如 L=[2,8,3,50], 输出L内所有数字的乘积末尾0的个数,
def getFactorNum(num,factor):
cnt = 0
while num%factor == 0:
cnt+=1
num /=factor
return cnt
fiveNum = 0
twoNum = 0
for num in L:
fiveNum += getFactorNum(num,5)
twoNum += getFactorNum(num,2)
print min(fiveNum,twoNum)
题目id：12 结尾非零数的奇偶性
给你一个正整数列表 L, 如 L=[2,8,3,50], 判断列表内所有数字乘积的最后一个非零数字的奇偶性,
def getFactorNum(num,factor):
cnt = 0
while num%factor == 0:
cnt+=1
num /=factor
return cnt
fiveNum = 0
twoNum = 0
for num in L:
fiveNum += getFactorNum(num,5)
twoNum += getFactorNum(num,2)
print ( 1 if fiveNum >= twoNum else 0)
题目id：13 光棍的悲伤
给你一个整数a，数出a在二进制表示下1的个数，并输出。
cnt = 0
while a:
cnt += a&1
a >>= 1
print cnt
题目id：14 Python之美
输出Python之禅
注意：输出python之禅的源码即可，不用转换为英文。
import this
print this.s
题目id：15 大小写转换
给定一个字符串a, 将a中的大写字母 转换成小写，其它字符不变，并输出。
print a.lower()
题目id：16 人民币金额打印
银行在打印票据的时候，常常需要将阿拉伯数字表示的人民币金额转换为大写表示，现在请你来完成这样一个程序。
在中文大写方式中，0到10以及100、1000、10000被依次表示为：
零壹贰叁肆伍陆柒捌玖拾佰仟万
以下的例子示范了阿拉伯数字到人民币大写的转换规则：
1 壹圆
11 壹拾壹圆
111 壹佰壹拾壹圆
101 壹佰零壹圆
-1000 负壹仟圆
1234567 壹佰贰拾叁万肆仟伍佰陆拾柒圆
现在给你一个整数a(|a|<100000000), 打印出人民币大写表示
d={0:'零',1:'壹',2:'贰',3:'叁',4:'肆',5:'伍',6:'陆',7:'柒',8:'捌',9:'玖'}
unit = {0:'',1:'拾',2:'佰',3:'仟',4:'万'}
b = abs(a)
def transferMoney(b,s,outOfRange):
s.append('万') if outOfRange else s.append('圆')
cnt = 0
flag = False
if b == 0:
s.append(d[0])
else:
while b:
if b%10:
flag = True
s.append(unit[cnt])
if flag:
s.append(d[b%10])
b/=10
cnt+=1
s=[]
transferMoney(b%10000,s,False)
if b >= 10000:
transferMoney(b/10000,s,True)
if a < 0:
s.append('负')
print "".join(s[::-1])
题目id：17 公约数的个数
给你两个正整数a,b, 输出它们公约数的个数。
def gcd(a,b):
if a < b:
a,b = b,a
while b:
a,b = b,a%b
return a
cnt = 1
for i in range(1,gcd(a,b)):
if a%i == 0 and b%i ==0:
cnt +=1
print cnt
题目id：18 逆解最大公约数与最小公倍数
我们经常遇到的问题是给你两个数，要你求最大公约数和最小公倍数。
今天我们反其道而行之，给你两个数a和b，计算出它们分别是哪两个数的最大公约数和最小公倍数。
输出这两个数，小的在前，大的在后，以空格隔开。若有多组解，输出它们之和最小的那组。
import sys
def gcd(a,b):
if a < b:
a,b = b,a
while b:
a,b = b,a%b
return a c= b/a+1
minv,mina,minb= sys.maxint,0,0
for i in range(1,c):
for j in range(i,c):
if gcd(i,j) == 1 and a*i*j == b:
if minv > (a*i+a*j):
minv,mina,minb = (a*i+a*j),a*i,a*j
print '%d %d'%(mina,minb)
题目id：19 单身情歌
抓不住爱情的我
总是眼睁睁看它溜走
...
现在来练习一下发现爱的能力，给你一个字符串a,如果其中包含"LOVE"（love不区分大小写)则输出LOVE，否则输出SINGLE。
print ("LOVE" if "love" in a.lower() else "SINGLE" )
题目id：20 信息加密
给你个小写英文字符串a和一个非负数b(0<=b<26), 将a中的每个小写字符替换成字母表中比它大b的字母。
这里将字母表的z和a相连，如果超过了z就回到了a。例如a="cagy",b=3, 则输出 fdjb
c=[]
for letter in a:
if letter.islower():
c.append(chr((ord(letter)+b-97)%26+97))
else:
c.append(letter)
print "".join(c)
题目id：21 回文子串
给你一个字符串a和一个正整数n,判断a中是否存在长度为n的回文子串。如果存在，则输出YES，否则输出NO。
回文串的定义：记串str逆序之后的字符串是str1，若str=str1,则称str是回文串，如"abcba".
def is_palindrome(a):
return (True if a == a[::-1] else False) for index in range(len(a)-n+1):
if is_palindrome(a[index:index+n]):
print 'YES'
break
else:
print 'NO'
题目id：22 时间就是金钱
给你两个时间st和et(00:00:00<=st <= et<=23:59:59), 请你给出这两个时间间隔的秒数。
如：st="00:00:00", et="00:00:10", 则输出10.
stTime = st.split(':')
etTime = et.split(':')
hour,minute,second = int(etTime[0])-int(stTime[0]),int(etTime[1])-int(stTime[1]),int(etTime[2])-int(stTime[2])
print hour*3600+minute*60+second
题目id：23 365 Or 366？
一年有多少天，这是个大问题，很值得思考。现在给你一个年份year(year为四位数字的字符串，如"2008","0012"),
你输出这一年的天数。如year="2013", 则输出365。
def isLeapYear(year):
return (year%4 == 0 and year%100 )or year%400 == 0
print (366 if isLeapYear(int(year)) else 365)
题目id：24 一马当先
下过象棋的人都知道，马只能走'日'字形（包括旋转90°的日），现在想象一下，给你一个n行m列网格棋盘，
棋盘的左下角有一匹马，请你计算至少需要几步可以将它移动到棋盘的右上角，若无法走到，则输出-1.
如n=1，m=2,则至少需要1步；若n=1，m=3,则输出-1
主要用深度优先搜索算法去做
n +=1
m +=1
visit = [([False]*n) for i in range(m)]
step = ((2,1),(1,2),(-1,2),(-2,1),(-2,-1),(-1,-2),(1,-2),(2,-1))
flag = False
cnt = 0
minCnt = n*m
def dfs(row,col):
global minCnt
global cnt
global flag
global step
global visit
if row == m-1 and col == n-1:
flag = True
if minCnt > cnt:
minCnt = cnt
return
for item in step:
newRow,newCol= row+item[0],col+item[1]
if newRow>=0 and newRow<m and newCol>=0 and newCol < n:
if visit[newRow][newCol] == False:
visit[newRow][newCol] =True
cnt +=1
dfs(newRow,newCol)
cnt -=1
visit[newRow][newCol] =False
visit[0][0] =True
dfs(0,0)
print (minCnt if flag else -1)
题目id：25 格式化时间
给你一个时间t（t是一个字典，共有六个字符串key(year,month,day,hour,minute,second),值为每个值为数字组成的字符串，
如t={'year':'2013','month':'9','day':'30','hour':'16','minute':'45','second':'2'}
请将其按照以下格式输出， 格式:XXXX-XX-XX XX:XX:XX。如上例应该输出： 2013-09-30 16:45:02。
d =[t['year'].zfill(4),t['month'].zfill(2),t['day'].zfill(2),t['hour'].zfill(2),t['minute'].zfill(2),t['second'].zfill(2)]
print '-'.join(d[:3]),':'.join(d[3:])
题目id：26 序列判断
给你一个整数组成的列表L，按照下列条件输出：
若L是升序排列的,则输出"UP";
若L是降序排列的,则输出"DOWN";
若L无序，则输出"WRONG"。
if L == sorted(L):
print 'UP'
elif L == sorted(L,reverse=True):
print 'DOWN'
else:
print 'WRONG'
题目id：27 加油站
每个加油站加油都有一个上限，保存在列表limit中，即limit[i]为第i个加油站加油的上限，
而从第i个加油站开车开到第(i+1)%n个加油站需要cost[i]升油,cost为一个列表。
现在有一辆开始时没有油的车，要从一个加油站出发绕这个公路跑一圈回到起点。
给你整数n，列表limit和列表cost,你来判断能否完成任务。
如果能够完成任务，输出起始的加油站编号，如果有多个,输出编号最小的。
如果不能完成任务，输出-1。
for i in range(n):
sumLimit,sumCost = 0,0
for k in range(n):
sumLimit +=limit[(k+i)%n]
sumCost += cost[(k+i)%n]
if sumCost > sumLimit:
break
else:
print i
break
else:
print -1
题目id：28 相同数字
给你一个整数列表L,判断L中是否存在相同的数字，
若存在，输出YES，否则输出NO。
print ('NO' if len(set(L)) == len(L) else 'YES')
题目id：29 判断三角形
给你三个整数a,b,c, 判断能否以它们为三个边长构成三角形。
若能，输出YES，否则输出NO。
d=[a,b,c]
d.sort()
if d[0]+d[1]>d[2] and abs(d[0]-d[1])<d[2]:
print 'YES'
else:
print 'NO'
题目id：30 National Day
马上国庆节了，用一个英文单词描述你此时此刻的心情。
print 'Happy'
题目id：31 山峰的个数
十一假期,小P出去爬山,爬山的过程中每隔10米他都会记录当前点的海拔高度(以一个浮点数表示),
这些值序列保存在一个由浮点数组成的列表h中。回到家中，小P想研究一下自己经过了几个山峰，请你帮他计算一下，输出结果。
例如：h=[0.9,1.2,1.22,1.1,1.6,0.99], 将这些高度顺序连线，会发现有两个山峰，故输出一个2(序列两端不算山峰)
cnt = 0
for i in xrange(1,len(h)-1):
if h[i]>h[i-1] and h[i]>h[i+1]:
cnt += 1
print cnt
题目id：32 三角形形状
给以一个三角形的三边长a,b和c(边长是浮点数),请你判断三角形的形状。
若是锐角三角形，输出R,
若是直角三角形，输出Z，
若是钝角三角形，输出D，
若三边长不能构成三角形，输出W.
d=[a,b,c]
d.sort()
if d[0]+d[1] > d[2] and abs(d[0]-d[1]) < d[2]:
if d[0]**2+d[1]**2 == d[2]**2:
print 'Z'
elif d[0]**2+d[1]**2 >d[2]**2:
print 'R'
else:
print 'D'
else:
print 'W'
题目id：33 大幂次运算
给你两个正整数a(0 < a < 100000)和n(0 <= n <=100000000000)，计算(a^n) % 20132013并输出结果
print pow(a,n,20132013)
题目id：34 密码生成
生活在当代社会，我们要记住很多密码，银行卡，qq，人人，微博，邮箱等等。小P经过一番思索之后，发明了下面这种生成密码方法：
给定两个正整数a和b, 利用a / b我们会到的一个长度无限的小数(若a / b不是无限小数，
比如1/2=0.5,我们认为0.5是0.5000000...，同样将其看做无限长的小数），小P将该小数点后第x位到第y位的数字
当做密码，这样，无论密码有多长，小P只要记住a,b,x,y四个数字就可以了，牢记密码再也不是那么困难的事情了。
现在告诉你a,b,x,y（0 < a,b <= 20132013, 0 < x <= y < 100000000000),请你输出密码。
例如：a = 1, b = 2, x = 1, y = 4, 则 a / b = 0.5000000..., 输出小数点后第1到4位数字，即5000
decimals = []
if a > b :
a =a%b
cnt = 0
flag = False
x-=1
y-=1
while cnt <= y:
decimalPart,a = divmod(a*10**1024,b)
if flag:
decimals.append(str(decimalPart).zfill(1024))
if x >= cnt and x < cnt+1024:
flag = True
decimals.append((str(decimalPart).zfill(1024))[x-cnt:])
if y>=cnt and y< cnt+1024:
break
cnt +=1024
if y/1024 !=x/1024:
decimals[-1]=(str(decimals[-1]))[:y-cnt+1]
else:
decimals[-1]=(str(decimals[-1]))[:y-x+1]
print ''.join(decimals)
题目id：35 最大连续子序列
给你一个整数list L, 如 L=[2,-3,3,50], 求L的一个连续子序列，使其和最大，输出最大子序列的和。
例如，对于L=[2,-3,3,50]， 输出53（分析：很明显，该列表最大连续子序列为[3,50]).
def maxSubSum(L):
maxSum,thisSum= 0,0
for item in L:
thisSum += item;
if thisSum > maxSum:
maxSum = thisSum
elif thisSum < 0:
thisSum = 0
return maxSum print maxSubSum(L)
题目id：36 最大非连续子序列
给你一个整数list L, 如 L=[2,-3,3,50], 求L的一个非连续子序列，使其和最大，输出最大子序列的和。
这里非连续子序列的定义是，子序列中任意相邻的两个数在原序列里都不相邻。
例如，对于L=[2,-3,3,50]， 输出52（分析：很明显，该列表最大非连续子序列为[2,50]).
注意类似这种情况L=[-1,-1,-1,-1,-1,0]
if len(L) >= 2:
dp=[L[0],max(L[0],L[1])]
for i in range(2,len(L)):
dp.append(max(dp[i-1],max(L[i]+dp[i-2],L[i])))
print dp[len(L)-1]
else:
print L[0]
题目id：37 简单题之勾股定理
给你直角三角形的两个直角边的边长a,b,请你求出其斜边边长，结果保留小数点后三位小数。
如a=3, b =4, 则输出5.000。
import math
print "%.3f"%math.sqrt(a**2+b**2)
题目id：38 简单题之列表转换
给你一个字符串列表L，请用一行代码将列表所有元素拼接成一个字符串并输出。
如L=['abc','d','efg'], 则输出abcdefg。
print ''.join(L)
题目id：39 简单题之输出格式练习
给你一个字符串列表L，用一行代码顺序输出L中的元素，元素之间以一个空格隔开，注意行尾不要有空格，输出单独占一行。
如L=['abc','d','efg'], 则输出abc d efg。
print ' '.join(L)
题目id：40 整数解
给你两个整数a和b（-10000<a,b<10000），请你判断是否存在两个整数，他们的和为a，乘积为b。
若存在，输出Yes,否则输出No
例如：a=9,b=15, 此时不存在两个整数满足上述条件，所以应该输出No。
import math
x1 = int(float(a+math.sqrt(a**2-4*b))/2)
x2 = int(float(a-math.sqrt(a**2-4*b))/2)
if x1+x2 == a and x1*x2 == b:
print 'Yes'
else:
print 'No'
题目id：41 Py数
Py从小喜欢奇特的东西，而且天生对数字特别敏感，一次偶然的机会，他发现了一个有趣的四位数2992，
这个数，它的十进制数表示，其四位数字之和为2+9+9+2=22，它的十六进制数BB0，其四位数字之和也为22，
同时它的十二进制数表示1894，其四位数字之和也为22，啊哈，真是巧啊。
Py非常喜欢这种四位数，由于他的发现，所以这里我们命名其为Py数。
现在给你一个十进制4位数n，你来判断n是不是Py数，若是，则输出Yes，否则输出No。
如n=2992，则输出Yes； n = 9999，则输出No。
def getSumByValue(n,value):
res = 0
while n:
res += n%value
n /=value
return res if getSumByValue(n,10) == getSumByValue(n,16) == getSumByValue(n,12):
print 'Yes'
else:
print 'No'
题目id：42 分拆素数和
把一个偶数拆成两个不同素数的和，有几种拆法呢？
现在来考虑考虑这个问题，给你一个不超过10000的正的偶数n，
计算将该数拆成两个不同的素数之和的方法数，并输出。
如n=10，可以拆成3+7，只有这一种方法，因此输出1.
import math
def is_prime(num):
for i in range(2,int(math.sqrt(num)+1)):
if num%i == 0:
return False
return True cnt = 0
for i in range(2,n/2): #注意1不是素数
if is_prime(i) and is_prime(n-i):
cnt +=1
print cnt
题目id：43 斐波那契数列
斐波那契数列为1,1,2,3,5,8...。数列从第三项起满足，该项的数是其前面两个数之和。
现在给你一个正整数n（n < 10000), 请你求出第n个斐波那契数取模20132013的值（斐波那契数列的编号从1开始）。
a,b=1,1
for i in range(3,n+1):
a,b=b%20132013,(a+b)%20132013
print b
题目id：44 超级楼梯
有一楼梯共n级，刚开始时你在第一级，若每次只能跨上一级或二级，要走上第n级，共有多少种走法？
现在给你一个正整数n（0<n<40),请你输出不同的走法数。
如n=2,则输出1（你只有一种走法，走一步，从第一级到第二级）
a,b=1,1
for i in range(3,n+1):
a,b=b,(a+b)
print b
题目id：45 砝码问题
有一组砝码，重量互不相等，分别为m1、m2、m3……mn；它们可取的最大数量分别为x1、x2、x3……xn。
现要用这些砝码去称物体的重量,问能称出多少种不同的重量。
现在给你两个正整数列表w和n， 列表w中的第i个元素w[i]表示第i个砝码的重量，列表n的第
i个元素n[i]表示砝码i的最大数量。i从0开始，请你输出不同重量的种数。
如：w=[1,2], n=[2,1], 则输出5（分析：共有五种重量：0,1,2,3,4）
母函数的问题
m = set()
for i in range(n[0]+1):
m.add(w[0]*i)
for i in range(1,len(n)):
tmp = m.copy()
m.clear()
for j in range(n[i]+1):
for k in tmp:
m.add(w[i]*j+k)
print len(m)
题目id：45 测试贡献题目
给你一个list L, 如 L=[2,8,3,50], 对L进行降序排序并输出, 如样例L的结果为[50,8,3,2]
print sorted(L,reverse = True)
题目id：46 取石子游戏
有两堆石子，数量任意，可以不同。游戏开始由两个人轮流取石子。游戏规定，每次有两种不同的取法，
一是可以在任意的一堆中取走任意多的石子；二是可以在两堆中同时取走相同数量的石子。最后把石子全部取完者为胜者。
现在给出初始的两堆石子的数目a和b，如果轮到你先取，假设双方都采取最好的策略，问最后你是胜者还是败者。
如果你是胜者，输出Win,否则输出Loose。
例如，a=3,b=1, 则输出Win(你先在a中取一个，此时a=2，b=1,此时无论对方怎么取，你都能将所有石子都拿走). 博弈论问题
import math
p = (math.sqrt(float(5))+1)/float(2)
c = abs(a-b)
if a > b:
a = b
if a == int(p*c):
print "Loose"
else: print "Win"
题目id：47 杨辉三角
还记得中学时候学过的杨辉三角吗？具体的定义这里不再描述，你可以参考以下的图形：
1
1 1
1 2 1
1 3 3 1
1 4 6 4 1
1 5 10 10 5 1
..............
先在给你一个正整数n，请你输出杨辉三角的前n层
注意：层数从1开始计数,每层数字之间用一个空格隔开，行尾不要有空格。
如n=2,则输出：
1
1 1
trangle = [[''],['','']] def print_trangle(n):
for i in range(n):
print ' '.join(trangle[i])
for i in range(3,n+1):
tmp=['']
newRow= trangle[i-2]
for j in range(len(newRow)-1):
tmp.append(str(int(newRow[j])+int(newRow[j+1])))
tmp.append('')
trangle.append(tmp) print_trangle(n)
题目id：48 砝码问题II
有一组砝码，重量互不相等，分别为m1、m2、m3……mn；每种砝码的数量有无限个。
现要用这些砝码去称物体的重量,给你一个重量n,请你判断有给定的砝码能否称出重量n。
现在给你一个正整数列表w和一个正整数n，列表w中的第i个元素w[i]表示第i种砝码的重量，
n表示要你判断的重量。如果给定砝码能称出重量n，输出Yes，否则输出No。
例如，w=[2,5,11], n=9,则输出Yes（取两个2，一个5）。
此类问题用母函数做比较简单
def solve():
num=[]
for i in range(len(w)):
num.append(n/w[i])
m = set()
for i in range(num[0]+1):
m.add(w[0]*i)
for i in range(1,len(num)):
tmp = m.copy()
m.clear()
for j in range(num[i]+1):
for k in tmp:
if w[i]*j+k < n:
m.add(w[i]*j+k)
elif w[i]*j+k == n:
return "Yes"
return "No" print solve()
题目id：49 进制转换
给你一个十进制数a，将它转换成b进制数,如果b>10,用大写字母表示（10用A表示，等等）
a为32位整数，2 <= b <= 16
如a=3,b = 2, 则输出11 不要忘记a为负数的情况
base = ('','','','','','','','','','','A','B','C','D','E','F')
def transfer(a,b):
res = []
aa = abs(a)
while aa:
res.append(base[aa%b])
aa/=b
if a<0:
res.append('-')
return ''.join(res[::-1])
print transfer(a,b)
题目id：50 Py扔铅球
Py不但是编程大牛，而且是运动健将。比如说扔铅球，1000m，现在Py参加校园扔铅球比赛，
给你Py的身高a（双精度数），球落地点与Py头部的连线与水平线的夹角 b（弧度），
要你编写一个程序计算Py扔铅球的水平距离。
a，b都是浮点数，注意b是弧度，其中， 140 < a < 200, 0 < b < 1.5.
输出你求出的水平距离，保留到小数点后三位。
如，a = 165.5, b=1.1, 则输出84.234
import math
print "%.3f" % (a/math.tan(b))
题目id：53 神の安♂排
记得有一次全班去唱K, 其中有个活动是情歌对唱. 具体操作流程是这样的:
准备好 21 个阄(我们班 15 男 6 女), 其中只有两个是有标记的, 每人随意抓取一个, 最后取到有标记的阄的两个人去点首情歌对唱.
旁边一哥们儿幽幽地对我说, 看来搅基真是神的安排啊, 你看我们班的男女人数, 搅基的几率 C(15,2)/C(21,2) 刚好是 1/2.
给跪了, 这哥们儿对数字太敏感了, 简直是拉马努金转世啊. 不过我随之想到一个问题: (21, 15) 真的是神的唯一安排吗? 其实不是的,
神还有很多类似的安排. 比如 (4, 3), 显然 C(4,2)/C(3,2) 也等于 1/2, 当然还有 (120, 85) 等等等等.
神的安排太多太多了, 如果我们定义 (n, m) 是一个安排(其中 1 < m < n), 而如果 C(m,2)/C(n,2) = 1/2, 它就是神的安排.
现在的问题是, 给你一个不大于 10^9 的正整数 N, 有多少组神的安排 (n, m) 满足 n <= N 呢? 具体参考 关于求解不定方程的n(n-1)=2m(m-1)的解法的总结
cnt = 0
x,y=1,1
while x<= 2*N-1:
x,y=3*x+4*y,2*x+3*y
if x > 2*N-1:
break
if x%2 and y%2:
cnt +=1
print cnt
题目id：57 那些年我们集过的卡片
不知道大家的童年有没有过和我相似的经历。我记得小时候经常买干脆面，不为别的，只是因为里面有一张人物卡片。
其实有很多这样的活动都有一个相同的模式：N 种人物卡片，每次买一包干脆面随机得到一张。当你集齐这 N 种人物时，就会有相应的奖励。
那时候还不懂怎么计算概率，白白给人家送了好多钱，吃了好多干脆面。
现在的任务是，给你一个正整数 N (1 <= N <= 10^4)，请你帮我从期望的角度计算平均需要买多少包干脆面才能集齐这 N 种人物。
提醒：由于结果可能不是整数，所以结果只保留到小数点后两位。
ans = 0
for i in range(N,0,-1):
ans+=N*1.0/i
print "%.2f" % ans