Python 100例——第五章---不定方程式の解

2043 ワード

# -*- coding: cp936 -*-
#   ----      
#《c    》36-41 
#21:45 2006-11-4
def z36():
    '''     
x+y+z=100
5*x+3*y+z/3=100
x,y,z        0'''
    #       7*x+4*y=100
    for i in range(1,25):
        t1=(100-4*i)
        t=int(t1/7)
        if t*7==t1:
            print [t,i,100-t-i]
    return
def findx(r2,r1):
    #  x  x*r2%r1=1
    y=r1
    if r2%r1==1:return 1
    q2,r1=divmod(r1,r2)
    q3,r2=divmod(r2,r1)
    q3=q3*q2+1
    t=False
    while  r1!=1:
        if r2==1 :r1-=1
        q1,r3=divmod(r1,r2)
        t,q2,q3,r1,r2=not t,q3,q1*q3+q2,r2,r3
    if t:q3=y-q3
    return q3
def get23(a,b):
    m=reduce(lambda x, y: x*y, a)
    t=map(lambda x,y=m: m/x, a)
    k=sum(map(lambda x,y,z:findx(x,y)*x*z,t,a,b))
    return k%m
def z37():
    '''     
x=1(mod 2)
x=2(mod 3)
x=4(mod 5)
x=0(mod 7)'''
    a,b=[2,3,5,7],[1,2,4,0]
    print get23(a,b)
def z38():
    #   x*5+y*2+z=100
    num=1
    for i in range(0,100,5):
        for j in range(0,100,2):
            for m in range(0,100):
                if i+j+m==100:
                    #print num,[i/5,j/2,m]
                    num+=1
    return
def z39():
    #       ,       26,  880,              
    for i in range(1,4):
        j=26-6*i
        j,yu=divmod(j,4)
        if yu==0 and j *(j + i)*(j + 2 * i)*(j + 3 * i) == 880:
            print "             :",i,j
            print "          :",[j + n * i for n in range(0,20)]
    return
def z40():
    #    12   ,  3 ,3 ,6 ,              
    num=1
    i=range(0,4)
    for j in i:
        for m in i:
            n=8-j-m
            if n<=6:
                print num,[j,m,n]
                num+=1
    return
def z41():
    '''   
x+y+z=30
3*x+2*y+z=50'''
    #       2*x+y=20
    for i in range(0,11):
        t1=(20-2*i)
        print [i,t1,30-t1-i]
    return
if __name__ == '__main__':
    s=""
    for i in range(36,42):
        s+='z'+str(i)+'()
' exec(s)