結果
問題 |
No.453 製薬会社
|
ユーザー |
![]() |
提出日時 | 2018-03-25 18:36:55 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 31 ms / 2,000 ms |
コード長 | 1,277 bytes |
コンパイル時間 | 93 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 11,008 KB |
最終ジャッジ日時 | 2024-06-25 06:07:35 |
合計ジャッジ時間 | 1,074 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 9 |
ソースコード
c,d=[int(i) for i in input().split()] gyou=5 retsu=7 matrix=[[0 for i in range(retsu)] for j in range(gyou)] #評価関数 matrix[0][0], matrix[0][1]= -1000, -2000 #x, yは正 matrix[3][0], matrix[4][1]= -1, -1 #薬品の使用量 matrix[1][0], matrix[1][1]= 3/4, 2/7 matrix[2][0], matrix[2][1]= 1/4, 5/7 #定数項 matrix[1][6], matrix[2][6]=c, d #スラック変数 matrix[1][2]=1 matrix[2][3]=1 matrix[3][4]=1 matrix[4][5]=1 for i in range(10): if min(matrix[0])>=0: break else: p_column=min(enumerate(matrix[0]), key=lambda x: x[1])[0] const_per_p=[10**6 for i in range(gyou)] for j in range(gyou): if matrix[j][p_column]>0: const_per_p[j]=matrix[j][6]/matrix[j][p_column] k_row=min(enumerate(const_per_p), key=lambda x: x[1])[0] pivot=matrix[k_row][p_column] new_matrix=[[0 for i in range(retsu)] for j in range(gyou)] for l in range(gyou): for j in range(retsu): if l==k_row: new_matrix[l][j]/=pivot else: new_matrix[l][j]=matrix[l][j]-matrix[l][p_column]*matrix[k_row][j]/pivot matrix=new_matrix[:][:] # print(k_row, p_column) #for i in range(5): # print(matrix[i]) print(matrix[0][6])