Skip to content

Instantly share code, notes, and snippets.

View junpeitsuji's full-sized avatar

Junpei Tsuji junpeitsuji

15 followers · 6 following
View GitHub Profile
@junpeitsuji
junpeitsuji / jonessatowadawiens.rb
Last active September 7, 2016 09:53
def jonessatowadawiens(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z)
res=(k+2)*(1-(w*z+h+j-q)**2 \
-((g*k+2*g+k+1)*(h+j)+h-z)**2 \
-(2*n+p+q+z-e)**2 \
-(16*((k+1)**3)*(k+2)*((n+1)**2)+1-f**2)**2 \
-((e**3)*(e+2)*((a+1)**2)+1-o**2)**2 \
-((a**2-1)*(y**2)+1-x**2)**2 \
-(16*(r**2)*(y**4)*(a**2-1)+1-u**2)**2 \
-(( (a+(u**2)*(u**2-a))**2 -1)*(n+4*d*y)**2 + 1-(x+c*u)**2)**2 \
@junpeitsuji
junpeitsuji / jonessatowadawiens.rb
Created September 7, 2016 08:54
def jonessatowadawiens(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z)
res=(k+2)*(1-(w*z+h+j-q)**2 \
-((g*k+2*g+k+1)*(h+j)+h-z)**2 \
-(2*n+p+q+z-e)**2 \
-(16*((k+1)**3)*(k+2)*((n+1)**2)+1-f**2)**2 \
-((e**3)*(e+2)*((a+1)**2)+1-o**2)**2 \
-((a**2-1)*(y**2)+1-x**2)**2 \
-(16*(r**2)*(y**4)*(a**2-1)+1-u**2)**2 \
-(((a+(u**2)*(u**2-a))**2-1)*(n+4*d*y)**2+1-(x+c*u)**2)**2 \
@junpeitsuji
junpeitsuji / parse.rb
Created August 7, 2016 08:49
# encoding: utf-8
array = []
File.open("leech.csv") do |io|
io.each do |line|
v = line.split(",")
v2 = v.map{|a| a.to_i}
array.push v2
@junpeitsuji
junpeitsuji / leech.py
Created August 7, 2016 08:40
import numpy as np
a = np.array([[8, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [4, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [4, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [4, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [4, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [4, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [4, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [2, 2, 2, 2, 2, 2, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [4, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [4, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [2, 2, 2, 2, 0, 0, 0, 0, 2, 2, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [2, 2, 0, 0, 2, 2, 0, 0, 2,
@junpeitsuji
junpeitsuji / leech.csv
Created August 7, 2016 08:23
8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0
@junpeitsuji
junpeitsuji / mathjax.inc.php
Created July 4, 2016 14:05
<?php
//
// pukiwiki用 数式プラグイン (mathjax.inc.php)
// Copywrite 2013 Osamu Takeuchi <[email protected]>
//
// [履歴]
// 2013.04.08 初期リリース
//
// [インストール]
// ソースファイルを (pukiwiki)/plugin/mathjax.inc.php として保存
@junpeitsuji
junpeitsuji / gaussian_prime.rb
Created May 27, 2016 12:02
# encoding: utf-8
require 'prime'
def fermat p, x, y
k = (x**2 + y**2)/p
b = (x % k)
@junpeitsuji
junpeitsuji / congruent.rb
Last active March 24, 2016 08:13
面積 n となる直角三角形の辺を発見するプログラム
# 面積 n となる直角三角形の辺を発見するプログラム
require 'rational'
# Rational class にメソッドを追加
class Rational
# 平方数かどうか判定する
@junpeitsuji
junpeitsuji / largest_dpn_check.cc
Last active February 22, 2016 06:51
periaさんの最大素数大富豪素数探索プログラムにおける結果 9999888877776666555544443333222213131313131312121212111111101101010111199998888777766665555444433332222131313131313121212121111111011010101111 の「最大性」をチェックするコード． ただし、最大性において以下の仮定をおいている．(1) JOKERはすべて "13" をもちいる（すなわち、13は6枚使用） (2) 相手の手元カードは "A" のみ（すなわち、1は3枚使用） 参考：Periaさんのコード https://gist.github.com/peria/00e294247aa422…
#include <gmpxx.h>
#include <algorithm>
#include <functional>
#include <iostream>
#include <string>
#include <vector>
using namespace std;
int main() {
@junpeitsuji
junpeitsuji / mpqs.rb
Last active May 11, 2018 04:00
複数多項式二次ふるい法（MPQS）
#複数多項式二次ふるい法（MPQS）
# の作りかけ
#
# 参考：
# https://en.wikipedia.org/wiki/Quadratic_sieve#Example_of_basic_sieve
# http://www.cs.t-kougei.ac.jp/nsim/lecture/2008/ws/QS.pdf
# http://www.asahi-net.or.jp/~KC2H-MSM/mathland/math12/math1207.htm
# http://inaz2.hatenablog.com/entry/2016/01/09/032521
# http://d.hatena.ne.jp/lemniscus/20130226/1361874593#special
#