python学习应用笔记（一）

       之前一直用c++写程序  所以考虑程序一般都比较容易往数据结构的方向想 而自己设计数据结构往往要费很大事  昨天看了一下python  发现脚本语言 真是厉害    用来进行模拟运算确实不错  可以先想个解决问题的思路 用python 先测试一下  如果可行 如果为了运行效率 使用c++改写  确实是一个不错的选择

      昨天在看操作系统的时候看到这么一个题目  

这个题目是操纵系统中典型的为多个进程分配资源 最后达到安全状态的题目,可以理解为找一条安全路径

如果使用c++来写 那就用深度优先搜索算法  然后就要慢慢思考解决办法了 , 为了练习刚刚学习的python 我使用python的 numpy来解决这个问题

因为还没有学怎么写函数  所以写成四个循环了 不过算法复杂度应该差不多, 准备再学一段时间后 优化一下

 

import numpy

a=numpy.array(([1,1,2,4],

   [0,1,2,2],

   [1,1,1,0],

   [1,1,1,1]))

b=numpy.array(([1,2,3,6],

   [1,1,2,2],

   [1,2,1,1],

   [1,1,2,3]))

c=numpy.array([0,1,0,1])



d=b-a





for i in range(0,4):

	j = 0

	for j in range(0,4):

		if d[i][j]>c[j]:

			break

#	print 'e'

#	if j==3 and d[i][j]<=c[j]:

	else:

		e=numpy.array([0,0,0,0])

		f=c

		e[i]=1

		f=f+d[i]

		d[i]=numpy.array((0,0,0,0))

		

		for k in range(0,4):

			if e[k]==1:

				continue

			n = 0

			for n in range(0,4):

				if d[k][n]>f[j]:

					break;

#			if n==3 and  d[k][n]<=f[j]:

			else:

				

				

				e[k]=1

				f=f+d[k]

				d[k]=numpy.array((0,0,0,0))

				for m in range(0,4):

					if e[m]==1:continue

					p = 0

					for p in range(0,4):

						if d[m][p]>f[p]:break

#					print p

#					print m

#					if p==3 and d[m][p]<=f[p]:

					else:

						

						

						e[m]=1

						f=f+d[m]

						d[m]=numpy.array((0,0,0,0))

						for q in range(0,4):

							if e[q]==1:continue

							r = 0

							for r in range(0,4):

								if d[q][r]>f[r]:break

#							if r==3 and  d[q][r]<=f[r]:

							else:

								print i+1

								print k+1

								print m+1

								print q+1

#								print 'heihei'

								exit()

								

						d[m]=b[m]-a[m]

						e[m]=0

				d[k]=b[k]-a[k]

				e[k]=0

		d[i]=b[i]-a[i]

		e[i]=0