VBA 使えるソースコード

つたないソースコードを載せます。これは、他人の書いたソースコードを読む練習に最適です。初心者の方は、どうしたらきれいになるかなど考えながら活用してください。

2020-03-01から1日間の記事一覧

Python 線形計画問題 その2

# -*- 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…

Python 線形計画問題 その1

# -*- 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…

Python 0-1ナップサック問題を解く

# -*- 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…

Python モンテカルロ法で円周率を求める

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…

Python ある整数以下の素数の数を数える

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))

Python 有向グラフとダイクストラ法

感覚を身に付けよう # -*- 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…

Python 無向グラフの生成 ノードと枝の確認

# -*- 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_…