Submission #856615

#	Time	Username	Problem	Language	Result	Execution time	Memory
856615		nnhzzz	Longest beautiful sequence (IZhO17_subsequence)	C++14	100 / 100	3218 ms	51544 KiB

subsequence

// cre: Nguyen Ngoc Hung - Train VOI 2024 :>

#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <exception>
#include <fstream>
#include <functional>
#include <iomanip>
#include <ios>
#include <iosfwd>
#include <iostream>
#include <istream>
#include <iterator>
#include <limits>
#include <list>
#include <locale>
#include <map>
#include <memory>
#include <new>
#include <numeric>
#include <ostream>
#include <queue>
#include <set>
#include <unordered_set>
#include <sstream>
#include <stack>
#include <stdexcept>
#include <streambuf>
#include <string>
#include <typeinfo>
#include <utility>
#include <valarray>
#include <vector>
#include <cstring>
#include <unordered_map>
#include <cmath>
#include <array>
#include <cassert>
#include <random>

using namespace std;

#define        __nnhzzz__  signed main()
#define          BIT(i,j)  ((i>>j)&1LL)
#define           MASK(i)  (1LL<<i)
#define            ALL(x)  (x).begin(),(x).end()
#define             SZ(x)  (int)(x).size()
#define                fi  first
#define                se  second
#define                ll  long long
#define                vi  vector<int>
#define               vvi  vector<vi>
#define               pii  pair<int,int>
#define              vpii  vector<pii>
#define REPDIS(i,be,en,j)  for(int i = (be); i<=(en); i+=j)
#define     REPD(i,be,en)  for(int i = (be); i>=(en); i--)
#define      REP(i,be,en)  for(int i = (be); i<=(en); i++)
//-------------------------------------------------------------//
const int oo = 1e9,LOG = 20,MAXN = 1e5+7,N = 1e2+3;
const int MOD = 1e9+7,MOD1 = 1e9+207,MOD2 = 1e9+407,MOD3 = 998244353;
//-------------------------------------------------------------//
template<typename T1, typename T2> bool mini(T1 &a, T2 b){if(a>b){a=b;return true;}return false;}
template<typename T1, typename T2> bool maxi(T1 &a, T2 b){if(a<b){a=b;return true;}return false;}
template<typename T> T gcd(T a, T b) { while(b) { a %= b; swap(a,b); } return a; }
template<typename T> T lcm(T a, T b) { return a/gcd(a,b)*b; }

/*
----------------------------------------------------------------
    END OF TEMPLATE
----------------------------------------------------------------
    Nguyen Ngoc Hung - nnhzzz
    Training for VOI24 gold medal
----------------------------------------------------------------
*/

int a[MAXN],k[MAXN],dp[MASK(10)|1][MASK(10)|1][11],cnt[MASK(20)|1];
int n;

void sub2(){
    vi dp(n+3),par(n+3),cnt(2100000,0);
    REP(i,0,SZ(cnt)-1) cnt[i] = __builtin_popcount(i);
    REP(i,1,n){
        dp[i] = 1;
        par[i] = 0;
    }
    REP(i,1,n){
        REP(j,1,i-1){
            int val = a[i]&a[j];
            if(cnt[val]==k[i]){
                if(maxi(dp[i],dp[j]+1)){
                    par[i] = j;
                }
            }
        }
    }
    int res = 0,pos = 0;
    REP(i,1,n){
        if(maxi(res,dp[i])) pos = i;
    }
    vi ans;
    while(pos!=0){
        ans.push_back(pos);
        pos = par[pos];
    }
    reverse(ALL(ans));
    cout << res << "\n";
    for(const auto &i:ans){
        cout << i << " ";
    }
}

void sub3(){
    vi pos(MASK(8)|1,-1),par(n+3,0),dp(n+3,1);
    REP(i,0,MASK(8)) cnt[i] = __builtin_popcount(i);
    REP(i,1,n){
        REP(j,0,MASK(8)){
            int val = a[i]&j;
            if(cnt[val]!=k[i] || pos[j]==-1) continue;
            if(maxi(dp[i],dp[pos[j]]+1)) par[i] = pos[j];
        }
        if(dp[pos[a[i]]]<dp[i]) pos[a[i]] = i;
    }
    vi ans;
    int res = 0,p = 0;
    REP(i,1,n){
        if(maxi(res,dp[i])) p = i;
    }
    while(p!=0){
        ans.push_back(p);
        p = par[p];
    }
    reverse(ALL(ans));
    cout << res << "\n";
    for(const auto &i:ans){
        cout << i << " ";
    }
}

void sub4(){
    // vector<vvi> dp(MASK(10)|1,vvi(MASK(10)|1,vi(11,-1)));
    REP(i,1,MASK(20)) cnt[i] = __builtin_popcount(i);
    memset(dp,-1,sizeof dp);
    vpii f(n+3);
    pii res = {1,-1};
    REP(i,1,n){
        int v1 = a[i]&(MASK(10)-1);
        int v2 = a[i]>>10;
        f[i] = {1,-1};
        REP(pre_v1,0,MASK(10)-1){
            int diff_v1 = __builtin_popcount(v1&pre_v1);
            int diff_v2 = k[i]-diff_v1;
            if(diff_v2<0 || diff_v2>10) continue;
            int cur = dp[pre_v1][v2][diff_v2];
            if(cur==-1) continue;
            maxi(f[i],make_pair(f[cur].fi+1,cur));
        }
        maxi(res,make_pair(f[i].fi,i));
        REP(nxt_v2,0,MASK(10)-1){
            int diff_v2 = cnt[v2&nxt_v2];
            int &val = dp[v1][nxt_v2][diff_v2];
            if(val!=-1 && f[i].fi<f[val].fi) continue;
            val = i;
        }
    }
    cout << res.fi << "\n";
    vi ans(1,res.se);
    while(f[ans.back()].se!=-1) ans.push_back(f[ans.back()].se);
    reverse(ALL(ans));
    for(const auto &i:ans){
        cout << i << " ";
    }
}

void solve(){
    cin >> n;
    REP(i,1,n) cin >> a[i];
    REP(i,1,n) cin >> k[i];
    if(n<=5000){
        sub2(); return ;
    }
    int ma = 0;
    REP(i,1,n) maxi(ma,a[i]);
    if(ma<=MASK(8)){
        sub3(); return ;
    }
    sub4();
}

__nnhzzz__{
    ios::sync_with_stdio(0);
    cin.tie(0); cout.tie(0);
    #define name "test"
    if(fopen(name".inp","r")){
        freopen(name".inp","r",stdin);
        freopen(name".out","w",stdout);
    }
    #define name1 "subsequence"
    if(fopen(name1".in","r")){
        freopen(name1".in","r",stdin);
        freopen(name1".out","w",stdout);
    }

    int test = 1;

    while(test--){
      solve();
    }
    cerr << "\nTime elapsed: " << 1000*clock()/CLOCKS_PER_SEC << "ms\n";
    return 0;
}

Compilation message (stderr)

subsequence.cpp: In function 'int main()':
subsequence.cpp:196:16: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
  196 |         freopen(name".inp","r",stdin);
      |         ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~
subsequence.cpp:197:16: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
  197 |         freopen(name".out","w",stdout);
      |         ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~
subsequence.cpp:201:16: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
  201 |         freopen(name1".in","r",stdin);
      |         ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~
subsequence.cpp:202:16: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
  202 |         freopen(name1".out","w",stdout);
      |         ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~

• Subtask 1: 7 / 7
• Subtask 2: 16 / 16
• Subtask 3: 17 / 17
• Subtask 4: 60 / 60