1~10000までの素数をグラフで表す(Python)

# 概要
　アルゴリズムを使ってある数までの素数を見つけ出し、配列に詰め込んでいくコードです。
　ついでに素数をグラフにして可視化しました。

```Python:sosuu.py

from numpy import linspace
import numpy as np
import matplotlib.pyplot as plt
import random

N = 1000

def factor(N):
    x = 2
    insuu = []
    sosuu = 0
    while x <= N ** 0.5:
        if N % x == 0:
            insuu.append(x)
        else:
            pass
        x += 1
        # print(insuu)
        
    if insuu == []:
        sosuu = N

    if sosuu != 0:
        return sosuu
    else:
        return 0

x = range(0, N)
y = []

for i in x:
    if factor(i) != 0:
        y.append(factor(i))
        #print(y)
    #else:
        #y.append(0)
#random.shuffle(y)
print(y)

for n in range(0, len(y)):
    color1 = [n/N, 0.0, 1-n/N]
    color2 = [1-n/N, n/N, 0.0]
    #print(n)
    #print(y[n])
    #plt.scatter(n, y[n], s=y[n], c=y[n], vmin=0, vmax=N, cmap='cool', marker="*", edgecolors="orange")
    p1 = plt.bar(n,y[n], bottom=n, color=color2)
    p2 = plt.bar(n,n, color=color1)
    plt.legend((p1[0],p2[0]),("sosuu","order"))
    
plt.show()
#plt.colorbar()

print('----------------------------------------------')

def sosuuFibonacci(array):
    num = 0
    sosuuFibonacci = []

    for i in range(0, len(array)-1):
        #sosuuFibonacci.append(array[i]+array[i+1])
        num += array[i]
        sosuuFibonacci.append(num)

    return sosuuFibonacci


y2 = sosuuFibonacci(y)
print(y2)

for n in range(0, len(y2)):
    color1 = [n/N, 0.0, 1-n/N]
    color2 = [1-n/N, n/N, 0.0]
    #print(n)
    #print(y[n])
    #plt.scatter(n, y[n], s=y[n], c=y[n], vmin=0, vmax=N, cmap='cool', marker="*", edgecolors="orange")
    p1 = plt.bar(n,y2[n], bottom=n, color=color2)
    p2 = plt.bar(n,n, color=color1)
    plt.legend((p1[0],p2[0]),("sosuuFibonacci","order"))
    
plt.show()
#plt.colorbar()

```

　コメントアウトしている部分は、デバッグ用の出力や、グラフなどの設定を変えたいとき用のものです。
　素数を使って色々なグラフが書けましたが、素数はしょせん素数で、あまり面白みのないグラフになりました。素数がどのくらいの割合で分布しているかなどは分かりやすいと思います。

　ついでに、フィボナッチ数列のように、各素数について、次の素数と足した数を配列にしたり、それまでの素数の和を配列にしたりしてみました。何か法則性を見つけたら、数学者を目指してください。あと私に教えてください。ちなみに、各素数とその次の素数の和は、5以外偶数なような気がします。あと、何故か6の倍数が67%程入っていました。

![sosuuNumber3.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/b09e4740-a7d3-f1ee-3389-72ec79ec937e.png)



　10000までの素数です。
![sosuuNumber1.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/ee20a146-0fa4-9917-0e02-543fec3a4dcb.png)


グラフはこんな感じになりました
二つ目は素数と次の素数を足したもののグラフ、三つめは、それまでの素数の合計のグラフです。
※4つ目のグラフは、同じように素数を使っていますが、全く別のコードからつくりました。他のグラフも、多少式や値を変えて作っています。

![sosuuNotFibonacci1.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/3ca41979-863e-3ae1-7a9e-0df3df28e2a7.png)


![sosuuFibonacci1.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/89d5eeb9-210c-e940-7495-b092347afe9a.png)

![sosuuFibonacci2.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/6faabbde-3667-f012-a371-9531aebc2b52.png)

![sosuu1.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/8becb44d-d687-4556-5ef4-317f242e72d8.png)

![sosuuGraph8.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/1b7f840b-8764-476f-fdc4-cf18b1d34379.png)

![sosuuGraph1.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/8a2f5bfb-e38d-9984-02d8-49681468bb56.png)

![sosuuGraph9.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/5c781992-b0e3-d43a-1249-29a2b83a9bf9.png)

![sosuuGraph10.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/6b0a9c27-6472-346c-700b-b2ced41a7d0a.png)

![sosuuGraph11.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/33662c96-b6d4-4019-ef0a-c8d98d60954f.png)

![sosuuGraph12.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/b0d851e5-7f93-2821-7e9a-856ac9bd0b7e.png)

![sosuuGraph15.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/8ee78453-2d6e-1879-bd4e-717a9a46fdfd.png)

![sosuuGraph17.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/e12ba3ce-6e83-f3a3-757f-5d42660b725d.png)

![sosuuGraph16.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/a04428a5-eab8-68cf-0e6f-efc8df8bf6e2.png)

![sosuuGraph13.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/3fbce6ff-5619-b16e-5e5b-9811a48a5e9b.png)


![sosuuGraph14.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2475646/8526cc12-eea9-1546-8e32-9f8d8eec9cde.png)


# 結論
　素数は面白い