CF1025D Recovering BST

题意：给定序列，问能否将其构成一颗BST，使得所有gcd(x, fa[x]) > 1

解：看起来是区间DP但是普通的f[l][r]表示不了根，f[l][r][root]又是n⁴的会超时，怎么办？

看了题解发现惊为天人......

f_l[l][r]表示[l, r]能否构成l-1的右子树，f_r[l][r]表示[l, r]能否构成r+1的左子树。

然后我们就发现这个神奇的东西变成n³了......

 #include <cstdio>

 const int N = ;

 int gcd(int ta, int tb) {
     if(!tb) {
         return ta;
     }
     return gcd(tb, ta % tb);
 }

 inline void read(int &x) {
     x = ;
     char c = getchar();
     while(c < '' || c > '') {
         c = getchar();
     }
     while(c >= '' && c <= '') {
         x = (x << ) + (x << ) + c - ;
         c = getchar();
     }
     return;
 }

 int a[N];
 bool G[N][N], val[N][N], var[N][N];

 int main() {
     int n;
     read(n);
     for(int i = ; i <= n; i++) {
         read(a[i]);
     }

     for(int i = ; i <= n; i++) {
         for(int j = i; j <= n; j++) {
             G[i][j] = G[j][i] = (gcd(a[i], a[j]) > );
         }
     }
     for(int i = ; i <= n; i++) {
         if(i > ) {
             val[i][i] = G[i][i - ];
         }
         if(i < n) {
             var[i][i] = G[i][i + ];
         }
     }

     for(int len = ; len < n; len++) {
         for(int l = ; l + len -  <= n; l++) {
             int r = l + len - ;
             for(int k = l; k <= r; k++) { // root
                 bool t = (k == l ?  : var[l][k - ]) & (k == r ?  : val[k + ][r]);
                 if(!t) {
                     continue;
                 }
                 if(G[k][r + ]) {
                     var[l][r] = ;
                 }
                 if(G[k][l - ]) {
                     val[l][r] = ;
                 }
             }
         }
     }

     for(int i = ; i <= n; i++) {
         bool t = (i ==  ?  : var[][i - ]) & (i == n ?  : val[i + ][n]);
         if(t) {
             printf("Yes");
             return ;
         }
     }
     printf("No");
     return ;
 }

AC代码