TDPC I イウィ

##問題
https://atcoder.jp/contests/tdpc/tasks/tdpc_iwi
##見た感じ
初めての区間dpと呼ばれる問題を解いた。
問題を見た感じから、
区間[i,j]をに絞った小問題についての解をdpしていけばよさそうだと思った。
##解法
```math
dp[i][j]=i番目からj番目までの文字列からとれるiwiの個数とあまりのpair\\
dp[i][j]=\max_{i\leq l<j}(dp[i][l]+dp[l+1][j]+(文字列を結合してiwiができれば+1))
```
この漸化式を文字数nを1からnまで増やしながら計算した。
##ソースコード
```c++:
//https://atcoder.jp/contests/tdpc/tasks/tdpc_iwi
//小さくできるものは小さくすべき！！！！
#include "bits/stdc++.h"

using namespace std;

#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define rep1(i, n) for (int i = 1; i <= (int)(n); i++)
typedef long long ll;
typedef pair<int,int> pint;
typedef pair<long,long> plong;

const double pi=M_PI;
const int INF10 = 1000000001, ID_MAX = 20;
const long long INF18 = 1e18 + 1;
const int mod = 1000000007;
//const int mod = 998244353;

struct edge
{
    int to;
    long cost;
};

void sort2array(double *a, long *b, long n)
{
    for (int i = 0; i < n; i++)
    {
        b[i] = i;
    }
    sort(b, b + n, [a](long i, long j) -> bool { return a[i] < a[j]; });
}
//bを何回足せばaを超えるか(O(a/b))
int wtover(int a, int b)
{
    if (a % b > 0)
        return (a / b) + 1;
    else
        return a / b;
}
int bi_e[ID_MAX] = {0};
//2進数表示したときの最高桁(O(log n))
int bi_max(long n)
{
    int m = 0;
    for (m; (1 << m) <= n; m++)
        ;
    m = m - 1;
    return m;
}
//bi_eに二進数表示したやつを代入(O(log^2 n))
void bi_exs(long n)
{
    memset(bi_e, 0, sizeof(bi_e));
    if (n < (1 << ID_MAX))
    {
        for (int i = 0; n > 0; n = (n >> 1), i++)
            bi_e[i] = n & 1;
    }
}
//x^n mod m (nが負の時は０)(O(log n))
long myPow(long x, long n, long m)
{
    if (n < 0)
        return 0;
    if (n == 0)
        return 1;
    if (n % 2 == 0)
        return myPow(x * x % m, n / 2, m);
    else
        return x * myPow(x, n - 1, m) % m;
}
// auto mod int
// https://youtu.be/L8grWxBlIZ4?t=9858
// https://youtu.be/ERZuLAxZffQ?t=4807 : optimize
// https://youtu.be/8uowVvQ_-Mo?t=1329 : division
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;
    }
    mint &operator-=(const mint a)
    {
        if ((x += mod - a.x) >= mod)
            x -= mod;
        return *this;
    }
    mint &operator*=(const mint a)
    {
        (x *= a.x) %= mod;
        return *this;
    }
    mint operator+(const mint a) const { return mint(*this) += a; }
    mint operator-(const mint a) const { return mint(*this) -= a; }
    mint operator*(const mint a) const { return mint(*this) *= a; }
    mint pow(ll t) const
    {
        if (!t)
            return 1;
        mint a = pow(t >> 1);
        a *= a;
        if (t & 1)
            a *= *this;
        return a;
    }

    // for prime mod
    mint inv() const { return pow(mod - 2); }
    mint &operator/=(const mint a) { return *this *= a.inv(); }
    mint operator/(const mint a) const { return mint(*this) /= a; }
};
istream &operator>>(istream &is, const mint &a) { return is >> a.x; }
ostream &operator<<(ostream &os, const mint &a) { return os << a.x; }
// combination mod prime
// https://www.youtube.com/watch?v=8uowVvQ_-Mo&feature=youtu.be&t=1619
struct combination
{
    vector<mint> fact, ifact;
    combination(int n) : fact(n + 1), ifact(n + 1)
    {
        assert(n < mod);
        fact[0] = 1;
        for (int i = 1; i <= n; ++i)
            fact[i] = fact[i - 1] * i;
        ifact[n] = fact[n].inv();
        for (int i = n; i >= 1; --i)
            ifact[i - 1] = ifact[i] * i;
    }
    mint operator()(int n, int k)
    {
        if (k < 0 || k > n)
            return 0;
        return fact[n] * ifact[k] * ifact[n - k];
    }
};

//分割数のmodの表生成(tof=falseで重複あり)d[i][j]=jのi分割
struct divNum
{
    vector<vector<mint>> d;
    int N, M;
    divNum(int n, int m, bool tof) : d(n + 1, vector<mint>(m + 1, 0))
    {
        N = n;
        M = m;
        if (tof)
        {
            rep(j, n + 1)
            {
                d[j][0] = 1;
            }
            rep(j, m)
            {
                d[0][j + 1] = 0;
            }
            for (int j = 1; j <= n; j++)
            {
                for (int k = 1; k <= m; k++)
                {
                    if (j > k)
                        d[j][k] = d[j - 1][k];
                    else
                        d[j][k] = d[j - 1][k] + d[j][k - j];
                }
            }
        }
        else
        {
            d[0][0] = 1;
            rep(i, n)
            {
                rep(j, m + 1)
                {
                    rep(k, m - j + 1)
                    {
                        d[i + 1][j + k] = d[i + 1][j + k] + d[i][j];
                    }
                }
            }
        }
    }
    mint operator()(int n, int m)
    {
        if (n < 0 || m < 0 || n > N || m > M)
            return 0;
        else
        {
            return d[n][m];
        }
    }
};
const int N_MAX=300;
string s;
pair<string,long> dp[N_MAX][N_MAX];

int matchstr(string a,string tem){
    bool tof=false;
    int res;
    rep(i,a.size()-tem.size()+1){
        if(a[i]==tem[0]){
            string x=a.substr(i,tem.size());
            if(x==tem){
                return i;
            }
        }
    }
    return -1;
}

long solve(){
    long n=s.size(),cou=0;
    string temp="iwi",x;
    pair<string,long> res[2];
    for(long i=1;i<=n;i++){
        if(i<=3){
            rep(j,n-i+1){
                string x=s.substr(j,i);
                if(x==temp){
                    dp[j][j+i-1]=pair<string,long>("",1);
                }else{
                    dp[j][j+i-1]=pair<string,long>(x,0);
                }
            }
            continue;
        }
        rep(j,n-i+1){
            pair<string,long> ans=pair<string,long>("",-1),a;
            for(long l=j;l<j+i-1;l++){
                x=dp[j][l].first+dp[l+1][j+i-1].first;
                a.second=dp[j][l].second+dp[l+1][j+i-1].second;
                int d=x.find(temp);
                if(d!=string::npos){
                    x.erase(d,3);
                    a.second++;
                }
                a.first=x;
                if(ans.second<a.second){
                    ans=a;
                }
            }
            dp[j][i+j-1]=ans;
        }
    }
    return dp[0][n-1].second;
}

int main()
{
    long i;
    //入力
    cin>>s;
    //処理
    if(s.size()<3){
        cout<<0<<endl;
        return 0;
    }
    long ans=solve();
    //出力
    cout<<ans<<endl;
}
```
##感想
区間について小問題が構成できただけで、dpの基本は変わらないみたい。なるべく早く解けるようになりたい