add gcd,mp,readme/ upd mint,vector

Author: atcoder-live

README.md

@@ -0,0 +1,38 @@
+AtCoder解説放送ライブラリ集
+----
+## これは何？
+[解説放送](https://www.youtube.com/channel/UCtG3StnbhxHxXfE6Q4cPZwQ)で作ったライブラリを公開しています。
+
+## 目次
+
+### ユーティリティ
+|名前|コード|説明|
+|:--|:--|:--|
+|テンプレート|[template.cpp](template.cpp)||
+|ModInt|[mint.cpp](mint.cpp)|自動でmodを取ってくれる整数型|
+
+### データ構造
+|名前|コード|説明|
+|:--|:--|:--|
+|BIT|[bit.cpp](bit.cpp)|Binary Indexed Tree (Fenwick Tree)|
+
+### 数学
+|名前|コード|説明|
+|:--|:--|:--|
+|GCD/LCM|[gcd.cpp](gcd.cpp)|最大公約数と最小公倍数|
+|Combination|[comb.cpp](comb.cpp)|nCkをmod素数で求める|
+
+### グラフ
+|名前|コード|説明|
+|:--|:--|:--|
+|LCA|[lca.cpp](lca.cpp)|最小共通祖先|
+
+### 文字列
+|名前|コード|説明|
+|:--|:--|:--|
+|KMP|[mp.cpp](mp.cpp)|文字列検索アルゴリズム（正確にはMP）|
+
+### 幾何
+|名前|コード|説明|
+|:--|:--|:--|
+|Vector|[vector.cpp](vector.cpp)|ベクトル（点を扱う際にも使う）|

gcd.cpp

@@ -0,0 +1,6 @@
+// GCD, LCM
+// https://youtu.be/8lm8o8L9Bmw?t=2285
+// https://youtu.be/XI8exXVxZ-Q?t=3595
+// https://youtu.be/F2p_e6iKxnk?t=843
+ll gcd(ll a, ll b) { return b?gcd(b,a%b):a;}
+ll lcm(ll a, ll b) { return a/gcd(a,b)*b;}

mint.cpp

@@ -6,6 +6,7 @@ const int mod = 1000000007;
 struct mint {
   ll x; // typedef long long ll;
   mint(ll x=0):x((x%mod+mod)%mod){}
+  mint operator-() const { return mint(-x);}
   mint& operator+=(const mint a) {
     if ((x += a.x) >= mod) x -= mod;
     return *this;

mp.cpp

@@ -0,0 +1,34 @@
+// Morris-Pratt
+// https://youtu.be/9MphwmIsO7Q?t=7283
+template<typename T>
+struct MP {
+  int n;
+  T t;
+  vector<int> a;
+  MP() {}
+  MP(const T& t): t(t) {
+    n = t.size();
+    a = vector<int>(n+1);
+    a[0] = -1;
+    int j = -1;
+    for (int i = 0; i < n; ++i) {
+      while (j != -1 && t[j] != t[i]) j = a[j];
+      j++;
+      a[i+1] = j;
+    }
+  }
+  int operator[](int i) { return a[i];}
+  vector<int> findAll(const T& s) {
+    vector<int> res;
+    int j = 0;
+    for (int i = 0; i < s.size(); ++i) {
+      while (j != -1 && t[j] != s[i]) j = a[j];
+      j++;
+      if (j == n) {
+        res.push_back(i-j+1);
+        j = a[j];
+      }
+    }
+    return res;
+  }
+};

geom.cpp renamed to vector.cpp

@@ -1,4 +1,4 @@
-// geometry
+// vector
 // https://youtu.be/UWbGRhF3Ozw?t=9564
 const double eps = 1e-9;
 struct V {