「行列プログラマー」の問題と課題を解いていくスレ [無断転載禁止]©2ch.net ->画像>1枚

https://www.amazon.co.jp/行列プログラマー-_Pythonプログラムで学ぶ線形代数-Philip-N-Klein/dp/4873117771

行列のできるプログラミング相談所はどこですか？

P10 問題 0.3.21
g:U→V、f:V→W、とする。 命題0.3.12により
(g^(-1)◦f^(-1))◦(f◦g)=g^(-1)◦(f^(-1)◦(f◦g))=g^(-1)◦((f^(-1)◦f)◦g)=g^(-1)◦(id_V◦g)=g^(-1)◦g=id_U
(f◦g)◦(g^(-1)◦f^(-1))=f◦(g◦(g^(-1)◦f^(-1)))=f◦((g◦g^(-1))◦f^(-1))=f◦(id_V◦f^(-1))=f◦f^(-1)=id_W
よってg^(-1)◦f^(-1)=(f◦g)^(-1)

買ってください！！！

文字化けしちゃう、、、いきなりの企画倒れ、、、orz

行列プログラマー Pythonプログラムで学ぶ線形代数
・著者：Philip N. Klein／松田 晃一、弓林 司、脇本 佑紀、中田 洋、齋藤 大吾
・定価：5,832 円(本体5,400 円＋税)
・B5　672頁
・ISBN 978-4-87311-777-5
・発売日 2016/10


買って作者のサポートをお願いします！
よろしくお願いします！

原著のタイトルは、『Coding the Matrix: Linear Algebra through Computer Science Applications』である。

その訳が『行列プログラマー ―Pythonプログラムで学ぶ線形代数』である。

訳者らが、どういうセンスの持ち主なのかがこのタイトルの訳だけで分かる。

P10 問題 0.3.22
g=e^x : R→(0,∞)
f=x^2 : R→R
f(g(x))は可逆だけどfは単射じゃない。

P19 課題 0.5.1
In [1]: 60*30*24*7
Out[1]: 604800

P19 課題 0.5.2
In [2]: 2304811-(2304811//47)*47
Out[2]: 25

P21 課題 0.5.4
In [4]
x=-9
y=1/2
2**(y+1/2) if x+10<0 else 2**(y-1/2)
Out[4]
1.0

P25 課題0.5.5
{x**2 for x in {1,2,3,4,5}}

P25 課題0.5.6
{2**x for x in {0,1,2,3,4}}

P25 課題0.5.7
{x*y for x in {1,2,3} for y in {5,7,11}}

P25 課題0.5.8
{x*y for x in {2,6,18} for y in {5,15,45}}

P26 課題0.5.9
S={1,2,3,4}
T={3,4,5,6}
{x for x in S if x in T}

P27 課題0.5.10
sum([20,10,15,75])/len([20,10,15,75])

P28 課題 0.5.11
[[x,y] for x in ['A','B','C'] for y in [1,2,3]]

P28 課題 0.5.12
sum(sum([[0.25,0.75,0.1],[-1,0],[4,4,4,4]],[]))

P30 課題 0.5.13
In [27]: [x,y]=[1,2,3]
---------------------------------------------------------------------------
ValueError Traceback (most recent call last)
<ipython-input-109-093f1be0fe93> in <module>()
----> 1 [x,y]=[1,2,3]

ValueError: too many values to unpack (expected 2)


[28]: [x,y,z]=[1,2]
---------------------------------------------------------------------------
ValueError Traceback (most recent call last)
<ipython-input-110-02e4cd0205e8> in <module>()
----> 1 [x,y,z]=[1,2]

ValueError: not enough values to unpack (expected 3, got 2)

P31 課題 0.5.14
S={-4,-2,1,2,5,0}
[(i,j,k) for i in S for j in S for k in S if i+j+k==0]

P31 課題 0.5.15
S={-4,-2,1,2,5,0}
[(i,j,k) for i in S for j in S for k in S if (i+j+k==0) & ((i,j,k)!=(0,0,0))]

P31 課題 0.5.16
（問題の意味が分からないからパス）

P32 課題 0.5.17
L=[0,1,1}

P32 課題 0.5.18
{2*x+1 for x in range(50)}

P33 課題 0.5.19
L=['A','B','C','D','E']
list(zip(range(5),L))

P33 課題 0.5.20
[x+y for (x,y) in zip([10,25,40],[1,15,20])]

早く課題を説いてください。
提出期限に間に合いません。

アメリカのSEは1,000万円以上だけど
日本のSEは1,000万円以下の低収入！

【アメリカ】
スーパープログラマ
時給10万円だったり、ストックオプションで数億〜数十億円稼ぎだす
Javaフレームワーク
〜3600万円
PHP
〜2400万円
COBOL
〜1700万円
技術サポート
〜1200万円

年収中央値：１１７５万円（アメリカ労働統計局調査、サンプルは１１１万人）
７０歳でも仕事があり、年収は下がらない
数年単位で転職する（一つの会社に長くいるのは危険）
管理系の職種は雇用が不安定で、報酬も高くない

【日本】
平均年収：４３０万円（情報処理推進機構調査）

Web/ゲーム業界（昔ほど報酬は高くない）以外は人月単位のため、報酬には上限あり

年功賃金を採用する企業では２０代後半までの給料は一部の例外を除き低い
間接雇用が基本（大手のSIerでも客先常駐派遣が少なくない）
４０歳以降になるとリストラ候補となり、一旦リストラされると低賃金職か、長期間無職となる
大企業の場合は管理職トラックに進むためコーディングはしなくなり、
プログラミング経験が昔あっても３５歳以降の転職は難しい
転職回数が３回超えるだけで大手には書類で落とす。

SEの経済損失

SI業界は善悪損益が世間と正反対！
プログラム生産するほど貧困非婚！

プログラムの巨額利益を客先に提供
プログラムの巨額報酬を人売に提供
会社員なのに短勤続年数
人手不足なのに無職意識
人手不足なのに低収入
高生産なのに低収入
高利益なのに低収入
高需要なのに低収入
勤勉なのに低収入
高稼働残業して非婚
高稼働残業して離婚
多重派遣なのに稼働
偽装請負なのに稼働

ピンハネ商流
http://freelance-programmer.com/shouryu

この本どんな感じ？数学苦手でも読めそう？

>>1は30ページ読んで挫折したのか？

>>21　まだ0章なので何とも言えない。
>>22　一応まだ挫折してないけど、問題の多さにかなり萎えてる。

P35 #0.5.21
dlist = [{'James':'Sean','director':'Terence'},{'James':'Rojer','dierctor':'Lewis'},{'James':'Pierce','director':'Rojer'}]
k = 'James'
[dic[k] for dic in dlist]

P35 #0.5.22
dlist = [{'Bilbo':'Ian','Frodo':'Elijah'},{'Bilbo':'Martin','Thorin':'Richard'}]
k = 'Bilbo'
print("k = Bilbo : ",[dic[k] if k in dic else 'NOT PRESENT' for dic in dlist])
k = 'Frodo'
print("k = Frodo : ",[dic[k] if k in dic else 'NOT PRESENT' for dic in dlist])

P35 #0.5.23
dic={x:x**2 for x in range(100) }
print(dic)

P36 #0.5.24
D={'red','white','blue'}
{x:x for x in D}

P36 #0.5.25
base=2
digits=set(range(base))
print("base=2 \n",{a*base**2+b*base+c :[a,b,c] for a in digits for b in digits for c in digits },'\n')

base=10
digits=set(range(base))
print("base=10 \n",{a*base**2+b*base+c :[a,b,c] for a in digits for b in digits for c in digits })

P37 d={0:1000.0, 3:990, 1:1200.50}
L=['Larry','Curly','','Moe']
{L[x]:d[x] for x in d.keys()}

[問題番号を追記]

P37 #0.5.26
d={0:1000.0, 3:990, 1:1200.50}
L=['Larry','Curly','','Moe']
{L[x]:d[x] for x in d.keys()}

# P37 0.5.27
def twice(z): return 2*z
print("twice(1)=",twice(1))
print("twice('x')=",twice('x'))
print("twice(['a','b','c'])=",twice(['a','b','c']))

# P37 0.5.28
def nextInts(L): return [x+1 for x in L]
a=[1,5,7]
nextInts(a)

# P38 0.5.29
def cubes(L): return [x**3 for x in L]
a=[1,2,3]
cubes(a)

# P38 0.5.30
def dict2list(dct,keylist): return [dct[keylist[x]] for x in range(len(keylist))]
x={'a':'A','b':'B','c':'C'}
y=['b','c','a']
dict2list(x,y)

# P38 0.5.31
def list2dict(L,keylist): return {keylist[x]:L[x] for x in range(len(keylist))}
x=['A','B','C']
y=['a','b','c']
list2dict(x,y)

# P38 0.5.31 zipを使って
def list2dict(L,keylist): return {x:y for x, y in zip(keylist,L)}
x=['A','B','C']
y=['a','b','c']
list2dict(x,y)

# P38 0.5.31　失敗例
def list2dict(L,keylist): return {x:y for x in keylist for y in L}
x=['A','B','C']
y=['a','b','c']
list2dict(x,y)

# P38 0.5.32
def all_3_digits_numbers(base,digits): return { a*base**2+b*base+c for a in digits for b in digits for c in digits}
b=10
d={x for x in range(b)}
print(all_3_digits_numbers(b,d))

# P39 0.6.1
import math
help(math)
print(math.sqrt(3)**2,math.cos(math.pi),math.log(math.e))
math.sqrt(-1) #これはエラーになる

# P40 0.6.2
from random import randint as ri
rev=["See it!","A gem","Ideological claptrap"]
def movie_review(name):
return name+" : " + rev[ri(0,2)]
movie_review("劇場版 探偵オペラ ミルキィホームズ 〜逆襲のミルキィホームズ〜")

# P40 0.6.3

# dictutil.py で保存
def dict2list(dct, keylist): return [dct[keylist[x]] for x in range(len(keylist))]
def list2dict(L, keylist): return {keylist[x]:L[x] for x in range(len(keylist))}

#同じディレクトリで下記を実行

import dictutil

x={'a':'A','b':'B','c':'C'}
y=['b','c','a']
print("dict2list",dictutil.dict2list(x,y))

x=['A','B','C']
y=['a','b','c']
print("list2dict",dictutil.list2dict(x,y))

# P41 0.6.4

# dictutil.py で保存
def dict2list(dct, keylist): return [dct[keylist[x]] for x in range(len(keylist))]
def list2dict(L, keylist): return {keylist[x]:L[x] for x in range(len(keylist))}
def listrange2dict(L): return {x:L[x] for x in range(len(L))}

#同じディレクトリで下記を実行
from imp import reload

reload(dictutil)

x=['A','B','C']
print("listrange2dict",dictutil.listrange2dict(x))

# P42 0.6.5
for x in [1,2,3]:
y=x*x
print(y)

# P44 0.6.6　全然pythonぽくない。誰かに書き直して欲しい。
def makeInverseIndex(strlist):

sl=strlist.replace('.',' ').replace(',',' ').replace(':',' ').replace(';',' ')
word_list=sl.split()
word_dic={wd:set() for wd in word_list}

f=open('stories_small.txt')
line_list=dict(enumerate(list(f)))

for num in line_list:
line_set = set(line_list[num].split())
for w in word_dic:
if w in line_set:
word_dic[w]=word_dic[w]|{num}
return word_dic

g=open('stories_small.txt')
linest=dict(enumerate(list(g)))

w_dic=makeInverseIndex(lines[0])
print(w_dic)

# P43 0.6.7
def orSearch(inverseIndex, query):
setObj=set()
for s in inverseIndex.keys:
for word in query:
word == s
setObj=setObj|inverseIndex[s]
return setObj

ss=makeInverseIndex('stories_small.txt')

print(orSearch(ss,["word","found","went"]))

# P44 0.6.8

def andSearch(inverseIndex,query):
check_word = "__z?z!"
keys=list(inverseIndex.keys())
while check_word in inverseIndex.keys():
check_word = check_word*2
and_set={check_word}
for k in inverseIndex.keys():
and_set=and_set|inverseIndex[k]
for w in query:
if keys.count(w)==0:
and_set=set()
break
and_set=and_set&inverseIndex[w]

return and_set

g=open('stories_small.txt')
lstr=""
for line in list(g):
lstr=lstr+' '+line

w_dic=makeInverseIndex(lstr)

w_list=['an','then']

orObj=andSearch(w_dic,w_list)
print(orObj)

インデントが消えてる、、、orz

【犯罪者追放のお願い】

違法派遣（偽装請負・多重派遣・偽装出向・事前面接等）についての刑事罰
【告訴権者＝業務委託、準委任、共同受注、業務請負契約および特定派遣（契約・正規）、一般派遣、正規社員】

�@職業安定法第４４条の労働者供給事業の禁止規定に違反（1年以下の懲役または20万円以下の罰金）
　■偽装請負・多重派遣・偽装出向・多重出向
　■事前面接（顔合わせ・面談・職場見学等）と履歴書・職務経歴書・スキルシート等提出による労働者の特定(※）
（音声録音で立証可能）
�A労働基準法第６条（中間搾取の禁止） (1年以下の懲役又は50万円以下の罰金)
　■多重派遣・多重出向

※違法派遣（派遣労働者の特定）→派遣法で認められた派遣労働者ではない→労働者供給事業→職業安定法44条違反というの
が前提となる法解釈となります。派遣法における罰則が軽微なのは法律の不備や労働者軽視などが原因ではありません。
違法派遣は全て職業安定法44条で裁くことが可能なため、刑罰の重複を避けるために派遣法には軽微な罰則（主に裁量行政による）しかないのです。

使用者に有利な民事訴訟や労働関係諸局への通報等の対極にあるのが書面（告訴状）による刑事告訴（※告訴先は検察の直告班）です。
労働関係諸局への通報・斡旋による軽微な「適正化」や監督・指導に対して、法律に定められた刑事罰を問うことになり、
違法派遣業者にとって有罪は考えられる限り最大の処罰となります。同時に刑事罰を受けた
担当者が取引先に与える悪印象を考慮すれば、通常会社側は告訴が受理された時点で告訴取り下げに
動くのが妥当でしょう。懲役、前科がつく刑罰が下される可能性から、告訴取り下げの和解金は高額となることが多いのです。

告訴の流れとしては、
刑事告訴⇒告訴受理⇒告訴取下げ要請⇒取下げ和解金入金⇒告訴取下げ
となります。告訴の懲役刑適応は犯罪者個人に対してのみですので、告訴する対象は
派遣先・派遣元 社員
派遣先・派遣元 社長
派遣先・派遣元 担当者・責任者・管理役員・取締役
派遣先・派遣元 人事管理担当者・人事管理役員・取締役
が妥当です。
刑事告訴取り下げの和解金額は犯罪者個人と交渉するとよいでしょう。（告訴状は人数分提出する必要あり）

rubyに行列ライブラリあるし
python使わなくたっていいじゃん

>>34
そもそも問題を読み間違えているのでは
grading.codingthematrix.com/edition1/inverse_index_lab.py

年収1,000万円以下はロジック障害者
迷惑だからSE・PG辞めろ

主な受注SEの結婚障害

無能残業して共働き妨害するな！
巨額搾取させて生活妨害するな！

・異性にモテない
・コミュニケーションが苦手
・ファッションセンスがない
・コンピューターが趣味
・プログラムの巨額利益を客先に提供
・プログラムの巨額報酬を人売に提供
・会社員なのに短勤続年数
・人手不足なのに無職意識
・人手不足なのに低収入
・高生産なのに低収入
・高利益なのに低収入
・高需要なのに低収入
・勤勉なのに低収入
・多重派遣なのに稼働
・偽装請負なのに稼働
・運動不足で不健康
・高稼働で不健康
・高稼働で家事困難
・低収入で生活困難

【SE】結婚障害【PG】
片働き共働き両方とも結婚困難
http://hanabi.2ch.net/test/read.cgi/infosys/1477540734/

買って今p7まで読み始めたけどp6の0.3.3が何言ってるかよく分からなくて飛ばした…

逃げちゃダメだダメだダメだ

【偽装請負・多重派遣の横奪被害】

巨額の料金を奪われてんだぞ！
横取り相手を訴えないのかよ！

裁判で横取り料金を取り戻せ！
盗難被害額分の作業を減らせ！

システム開発報酬盗難被害の事件例

加害者
発注者 支払 140万円/人月 1億円/人月の大儲け

被害者
1次受注者 報酬 120万円/人月 20万円/人月の盗難被害
2次受注者 報酬 80万円/人月 60万円/人月の盗難被害
3次受注者 報酬 60万円/人月 80万円/人月の盗難被害

>>6-7
各地で行列グルメを楽しみつつ、行列中にノーパソでさくさくプログラムを書いている、放浪プログラマーのことかと思ってここに来てみた

無能実態派遣残業して時間報酬相場下げるな！
【知的財産と契約料金の搾取促進者ばかり】
［SI生涯損害助長SEを追放すべき］
偽装請負従犯SEの動機
コミュニケーション障害
コンピュータ趣味
人格障害
文系大卒
低偏差値大卒
情報処理資格非保有者

偽装請負従犯SEの損害
無償プログラム提供
事前面接
契約外期限遵守
客先指示遵守
知的財産譲渡
中間搾取促進
時間外労働違反
低予備工数見積
残業見積
無料追加
学習不足
裁判苦手
対人障害
健康障害
孤独死

偽装請負従犯SEの代償
低収入低技術
非婚離婚
鬱病早死

無能実態派遣残業して時間報酬相場下げるな！
【知的財産と契約料金の搾取促進者ばかり】
［SI生涯損害助長SEを追放すべき］
偽装請負従犯SEの動機
コミュニケーション障害
コンピュータ趣味
人格障害
文系大卒
低偏差値大卒
情報処理資格非保有者

偽装請負従犯SEの損害
無償プログラム提供
事前面接
契約外期限遵守
客先指示遵守
知的財産譲渡
中間搾取促進
時間外労働違反
低予備工数見積
残業見積
無料追加
学習不足
裁判苦手
対人障害
健康障害
孤独死

偽装請負従犯SEの代償
低収入低技術
非婚離婚
鬱病早死

P252課題4.15.2 ってどう書き始めれば良いんだ？
引数をとらずにどう操作する行列のデータを持ってくれば良いのか分からん

【主な偽装請負従犯結婚障害者の作業】
[文系多数の使い捨て貧困非婚スキル]
コマンド
スクリプト
データ > ロジック
簡単ロジック
大量データ
SE適性不要
IT資格不要
大卒資格不要
文科系対象
体育系対象
商業系業種
業務系処理

[技術不要の使い捨て貧困非婚ソフト]
ノンプログラミングツール
フレームワーク
Web
COBOL
VＢ
.net
Java
DB
ERP
SAP

ちょっと質問
この本興味あってまだ買ってないんだけど
課題って答えは書いてないの？

主なSEの結婚障害

無能残業して共働き妨害するな！
巨額搾取させて生活妨害するな！

・常識がない
・モラルがない
・異性にモテない
・ファッションセンスがない
・コミュニケーションが苦手
・コンピューターが趣味
・プログラムの巨額利益を客先に提供
・プログラムの巨額報酬を人売に提供
・会社員なのに短勤続年数
・人手不足なのに無職意識
・人手不足なのに低収入
・高生産なのに低収入
・高利益なのに低収入
・高需要なのに低収入
・勤勉なのに低収入
・多重派遣なのに稼働
・偽装請負なのに稼働
・PC使用過多で不健康
・運動不足で不健康
・高稼働で不健康
・高稼働で家事困難
・低収入で生活困難

【SE】結婚障害【PG】
http://hanabi.2ch.net/test/read.cgi/infosys/1477540734/

>>53
https://github.com/franzip/coursera/tree/master/coding-the-matrix

公式のは無いけどこうやって解いて公開してくれてる人はいる
課題4.5.12のもweek4>geometry labにあったぞ

このコピペ荒しなんとかならんのか

>>56
貧困非婚なんとかならんのか

【犯罪者追放のお願い】

偽装請負とは、実際は人材を派遣して利益を得ているにも関わらず、業務請負など別の契約形態で労働者を働かせることを言い以下に該当します。
●労働者派遣法違反
●職業安定法第４４条違反
●労働基準法第６条違反

偽装請負の相談窓口
東京労働局需給調整事業第二課（海岸庁舎）
03-3452-1474
神奈川労働局需給調整事業課
045-650-2810
045-650-2880

来所の際に下記証拠品をご持参ください
・現場入館証のコピー
・顧客とやり取りしたメール
・雇用契約書
・委託契約書
・時間指定
・場所指定
・方法指定
・就業規則
・給与明細

>>55
githubで公開してる人がちょくちょくいるんだね
ありがとう

>>17

-4 + 2 + 2 == 0
2 - 4 + 2 == 0
2 + 2 - 4 == 0

のように解の組が重複しています。

重複解のうち一番最初に現れたもののみを表示せよということです。

>>51

def identity():
■return Mat(({'x','y','u'},{'x','y','u'}),{('x','x'):1,('y','y'):1,('u','u'):1})

>>45

集合 D から 集合 F への関数全体の集合を F^D と表記するという話です。

要するに約束です。

で、なぜこのように約束するかという理由が書いてあります。

D = {a, b}
F = {1, 2, 3}

とします。

F^D = {f1, f2, f3, f4, f5, f6, f7, f8, f9}

f1(a) = 1
f1(b) = 1

f2(a) = 1
f2(b) = 2

f3(a) = 1
f3(b) = 3

f4(a) = 2
f4(b) = 1

f5(a) = 2
f5(b) = 2

f6(a) = 2
f6(b) = 3

f7(a) = 3
f7(b) = 1

f8(a) = 3
f8(b) = 2

f9(a) = 3
f9(b) = 3

です。

|F| = 3
|D| = 2

です。

ですので、

|F|^|D| = 3^2 = 9 = |F^D|

が成り立ちます。

とても見やすい式になっています。

これが、集合 D から 集合 F への関数全体の集合を F^D と表記するように
約束した理由です。

>>34-37

def makeInverseIndex(strlist):
■dic = {}
■for doc in strlist:
■■for s in doc.split():
■■■if s not in dic:
■■■■dic[s] = set()
■for s in dic.keys():
■■for (i, doc) in enumerate(strlist):
■■■if s in doc.split():
■■■■dic[s].add(i)
■return dic

def orSearch(inverseIndex, query):
■s = set()
■for word in query:
■■s |= inverseIndex[word]
■return s

def andSearch(inverseIndex, query):
■s = inverseIndex[query[0]]
■for word in query[1:]:
■■s &= inverseIndex[word]
■return s

Gilbert Strang の本とどっちがいいですか？

WandboxやIdeoneやpaiza.IO使えばいいのに

p68 課題1.5.1
鍵が同じだから総当たりすればいいやと思ったんだけど
肝心な平文がわかんない、どれかな？

VEVLZDLVEZ
WFWM EMWF A
XGXNAFNXGAB
YHYOBGOYHBC
ZIZPCHPZICD
J QDIQ JDE
AKAREJRAKEF
BLBSFKSBLFG
CMCTGLTCMGH
DNDUHMUDNHI
EOEVINVEOIJ
FPFWJOWFPJK
GQGXKPXGQKL
HRHYLQYHRLM
ISIZMRZISMN
JTJ NS JTNO
KUKAOTAKUOP
LVLBPUBLVPQ
MWMCQVCMWQR
NXNDRWDNXRS
OYOESXEOYST
PZPFTYFPZTU
Q QGUZGQ UV
RARHV HRAVW
SBSIWAISBWX
TCTJXBJTCXY
UDUKYCKUDYZ

何か間違ってない？
自分も総当りで調べたけど下記になったよ

key : 10001
decode : ['00100', '10101', '00100', '11010', '01000', '10010', '11010', '00100', '10101', '01000', '01011']
EVE IS EVIL

>>69
返事ありがとうございます
前のページのGF(2)勘違いして
暗号化された文に0-26を足していって、それを%27してました

>>45
訳者が何もわかってない馬鹿ばっかりだから
ウソばっかり書いてあるのさ。
つか著者も馬鹿だろ？

こんな簡単な線形代数を間違いだらけ。
どんだけ馬鹿なんだよ。

だが、悪い本ではない。
作者の間違いを超えていくことで、
間違いだらけの仕様書で作っていくということの訓練になる。
仕事で疲労困憊して脳内変換しているのに、
勉強でも疲労困憊するわけだ。

楽しそうだな。

9章まで精読した。

簡単だった。

p.138の問題2.14.10でvec.pyの中身を書いたんだけど、意味不明のエラーが発生しました。

下のコメントアウトしたコードだとdoctestでエラーが発生する。
上のコードだとエラーが発生しない。

上のコードと下のコードは等価だと思うんだけどなんででしょうか？

def neg(v):
"""
Returns the negation of a vector.

>>> u = Vec({2, 4, 6, 8}, {2:1, 4:2, 6:3, 8:4})
>>> neg(u)
Vec({8, 2, 4, 6}, {8:-4, 2:-1, 4:-2, 6:-3})
>>> u == Vec({2, 4, 6, 8}, {2:1, 4:2, 6:3, 8:4})
True
"""
return scalar_mul(v, -1)


#def neg(v):
# """
# Returns the negation of a vector.

# >>> u = Vec({2,4,6,8},{2:1,4:2,6:3,8:4})
# >>> neg(u)
# Vec({8, 2, 4, 6},{8: -4, 2: -1, 4: -2, 6: -3})
# >>> u == Vec({2,4,6,8},{2:1,4:2,6:3,8:4})
# True
# """
# return scalar_mul(v, -1)

逆でした。

下のコメントアウトしたコードだとdoctestでエラーが発生しない。
上のコードだとエラーが発生する。

エラーメッセージは以下です：

**********************************************************************
File "C:\Users\***\Documents\Python\module\行列プログラマー\vec.py", line 171, in vec.neg
Failed example:
neg(u)
Expected:
Vec({8, 2, 4, 6}, {8:-4, 2:-1, 4:-2, 6:-3})
Got:
Vec({8, 2, 4, 6},{8: -4, 2: -1, 4: -2, 6: -3})
**********************************************************************
1 items had failures:
1 of 3 in vec.neg
***Test Failed*** 1 failures.
Out[15]:
TestResults(failed=1, attempted=48)

上下で異なるのは、テストコードの見た目（半角スペース）だけです。

もしかして、doctestでの結果の比較って出力の文字列を比較しているんですか？

考えてみれば文字列を比較しているはずですね。

クラインの描くグラフが分かり肉すぎる。

Jupyter Notebook で解答集を作ったんだけど、どうやって公開すればいいの？

p.197に

（１）
出力：　Ａ＊ｕ　＝　v を許容誤差の範囲内で満たすあるベクトル u が存在するとき、
そのベクトルを返す。

と書いてあります。その下に

（２）
solve(A, b) から得られた解答 u に対し、 A*u と b を比べて、それが解であるかどうか
チェックする必要がある。

と書いてあります。

（１）から、チェックはプロシージャ内で行われていると解釈できると思います。

これはどういうことでしょうか？

p.205 問題4.6.12ってグラフアルゴリズム知らなくても余裕じゃん。

プロットすると

SVGタグで直接プロットデータが出力されるますが、これを画像ファイルとして
出力するにはどうすればいいんですか？

https://github.com/GoroSatake/coding-the-matrix

4章までの解答があるね。

>>83
完璧だな。

行列ができるプログラマー相談所

誰でも簡単にパソコン1台で稼げる方法など
参考までに、
⇒　『宮本のゴウリエセレレ』 というブログで見ることができるらしいです。

グーグル検索⇒『宮本のゴウリエセレレ』

VNZ7TTITFR

とても簡単な自宅で稼げる方法
参考までに書いておきます
グーグルで検索するといいかも『ネットで稼ぐ方法　モニアレフヌノ』

APM5G

>>7
訳者つか編集だな
技術も知識もない文系

>>8
f(g) は単射でないことから、ある実数　a, bが存在して
f( g(a) ) = f( g(b) ) = c

よって( f(g) )^-1( c ) = { a, b, .... }

以上より f(g) は可逆じゃねーぞカス