トラブル内容
こちらのサンプルでフェイスメッシュして遊ぼうとしたら、
エラー吐いて落ちる現象。
AttributeError: module 'mediapipe.python.solutions.face_mesh' has no attribute 'FACE_CONTOURS'
環境
Windows 10 64bit
Python 3.9.1
解決方法
海外サイトを漁っていたら見つけた情報。
「Update: FACEMESH_CONTOURS is the correct new name.」
サンプルのfacemesh.pyの３５行目あたり、
FACE_CONTOURSをFACEMESH_CONTOURSに書き換えたらすんなり動いた。
connections=mp_face_mesh.FACEMESH_CONTOURS ,
モチベーション
備忘録です。
for文を使って処理した値を空のリストに入れるというような操作をpytorchでやりたかったのですが、pytorchで空の配列を作ろうとすると処理後のtensorの次元にあったものを作って代入する必要があるっぽい。しかしパラメータによって処理後のtensorの次元が異なっていてそれに合わせた空のtensorを作るのはあまりにも面倒くさいと頭を抱えていたところようやくやり方を見つけました。
参考URL : Appending to a tensor
僕が頭を抱えていたところ
ここは僕個人の問題なので読み飛ばしていただいて構いません。
１次元畳み込みtorch.Conv1dをつかってチャネルごとの特徴量を抽出した後、それらを結合したかった。
やり方
一般的なやり方通り空の配列[]を作る
forループの中で処理したtensorを空の配列にappendする
forループ終了後はtorch.catでリストの中身のtensorをくっつけたい次元を指定してconcatenate.
model.py
features = []
for i in range(x.shape[1]):
x_conved = cnn1d(x[ :, i, :])
features.append(x_conved)
features = torch.cat(features, dim = 1)
なんでこんな簡単なこと思いつかなかったんだろう・・・
本記事はProjectEulerの「100番以下の問題の説明は記載可能」という規定に基づいて回答のヒントが書かれていますので、自分である程度考えてみてから読まれることをお勧めします。
問題 41. パンデジタル素数
原文 Problem 41: Pandigital prime
問題の要約：n桁のパンデジタル数で素数のものの最大値を求めよ
パンデジタル数はもう何度も出てきていますが、今回は素数。例として素数2143が4桁パンデジタル数として挙げられています。これも最大値を求めよなので大きい順に探していきます。
import itertools
import sympy
digits9 = "123456789"
found = False
for n in range(9,0,-1):
digits = digits9[n-1::-1] # "987654321", "87654321",,,"21","1"
for pstr in itertools.permutations(digits): # pandigital number in decending order
ps = ''.join(pstr)
if sympy.isprime(int(ps)):
found = True
break
if found: break
print(f"Answer: {ps}")
（開発環境：Google Colab）
本記事はProjectEulerの「100番以下の問題の説明は記載可能」という規定に基づいて回答のヒントが書かれていますので、自分である程度考えてみてから読まれることをお勧めします。
問題 40. チャンパーノウン定数
原文 Problem 40: Champernowne's constant
問題の要約：以下のように小数点以下のに自然数を 1 から小さい順に並べた少数点以下n桁目の数を$d_n$としたときの$d_1 × d_{10}× d_{100} × d_{1000} × d_{10000} × d_{100000} × d_{1000000} $を求めよ
0.12345678910112131415161718192021... \\
チャンパーノウン定数(Wikipedia) に詳しい説明があります。プログラム的にはかなりシンプルです。
import numpy as np
def Cham(maxlen): # generate Champernowne's constant as long as maxlen
ret = ""
for n in range(1,N+1):
ret += str(n)
if len(ret)>maxlen: break
return ret
print(f"Answer is {np.prod([int(Cham(10**6)[10**i-1]) for i in range(NP+1)])}")
（開発環境：Google Colab）
本記事はProjectEulerの「100番以下の問題の説明は記載可能」という規定に基づいて回答のヒントが書かれていますので、自分である程度考えてみてから読まれることをお勧めします。
問題 39. 整数の直角三角形
原文 Problem 39: Integer right triangles
問題の要約：整数の直角三角形の周囲の和pとするとp<1000で最も多いpを求めよ
整数の直角三角形すなわちピタゴラス数の生成は「Problem 9: 特別なピタゴラス数」で作ったpythを使います。生成された基本ピタゴラス数の周囲の和pの1000以下の倍数をすべて配列に入れ、これをCounterに入力して数えて最も多いものをmost_commonメソッドで取り出します。ただし要注意なのはpythは必ずしも周囲の和pの小さい順には生成しないのでp=1000で止めると、その後に1000より小さいものが出てくる可能性があるということです。
from collections import Counter
N, p, peris = 1000, 0, []
pgen = pyth()
while p > N*2: # need up to twice of the bound
p = sum(next(pgen))
for i in range(1,(N//p)+1): # Multiples up to N
peris.append(p*i)
print(Counter(peris).most_common()[0]) # extract most common number
そこでpをいくつまで生成すればいいか考えてみました。
pythでpを$10^6$まで生成して、p以前のの最大値pmaxとの比率p/pmaxの最小値rを求めます。数が大きくなっていくと0.5に収束していくようです。したがって1000以下のpをすべて生成するためには2倍の2000になるまでという条件で良さそうです。たぶん証明すれば出来るのでしょうが、誰か出来たら教えてください。
pgen = pyth()
p, pmax, r = 0, 1, 1
while p < 10**6:
p = sum(next(pgen))
if p>pmax: pmax=p
if (p/pmax) < r: r = (p/pmax)
print(f"Minimum s/smax ratio = {r}")
#Minimum s/smax ratio = 0.5035145516710452
（開発環境：Google Colab）
