2020-03-01から1日間の記事一覧
# -*- coding: utf-8 -*- from pulp import LpProblem, LpVariable, LpMaximize, LpStatus problem = LpProblem('sample', LpMaximize) x1 = LpVariable('x1', lowBound=0,upBound=None,cat='Continuous')x2 = LpVariable('x2', lowBound=0,upBound=None,cat…
# -*- coding: utf-8 -*- from pulp import LpProblem, LpVariable, LpMaximize, LpStatus problem = LpProblem('sample', LpMaximize) x1 = LpVariable('x1', lowBound=0, upBound=2,cat='Continuous')x2 = LpVariable('x2', lowBound=0, upBound=3,cat='Co…
# -*- coding: utf-8 -*- from pulp import LpProblem, LpVariable, LpMaximize, LpStatus,lpDot cost = [700,500,400,800,600,200,300]profit = [2700,1650,1200,2600,1800,1000,1200] upCost = 2000r = range(len(cost)) problem = LpProblem('sample',LpM…
import randoma = 0b = 1h = 1n = 100000count = 0for i in range(n): x = random.uniform(a,b) y = random.uniform(0,h) if x**2 + y**2 < 1: count += 1print("領域内の点の数 :" + str(count) + '/' + str(n))print("積分値の近似値 :" + str*1 / n))prin…
import math def numOfPrime(n): m = 0 def isPrime(s): for i in range(2,int(math.sqrt(s) + 1)): if s % i == 0: return False return True for j in range(2,n+1): if isPrime(j): m += 1 return m print(numOfPrime(100))
感覚を身に付けよう # -*- coding: utf-8 -*- import networkx as nximport matplotlib.pyplot as plt G = nx.DiGraph()#空の有向グラフのインスタンスの作成 G.add_node('a')#ノードの追加G.add_node('b')G.add_node('c')G.add_node('d')G.add_node('e')G.a…
# -*- coding: utf-8 -*- import networkx as nximport matplotlib.pyplot as plt G = nx.Graph()#空の無向グラフのインスタンスの作成 G.add_nodes_from([1,2,3,4,5])#ノードG.add_edges_from([(1,2),(3,4)])#辺 nx.draw_circular(G,with_labels=True,node_…