#471192 (C++17) No.382 シャイな人たち (2)

• 問題: No.382 シャイな人たち (2)
• ユーザー: Ryuhei Mori
• 提出日時: 2020-04-26 10:14:51
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 1,376 ms / 8,000 ms
• コード長: 1,427 bytes
• コンパイル時間: 333 ms
• コンパイル使用メモリ: 34,176 KB
• 実行使用メモリ: 5,248 KB
• 最終ジャッジ日時: 2024-11-16 00:40:52
• 合計ジャッジ時間: 25,599 ms
• ジャッジサーバーID （参考情報）: judge2 / judge4

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	1,376 ms; 5,248 KB
testcase_01	AC	1,367 ms; 5,248 KB
testcase_02	AC	1,160 ms; 5,248 KB
testcase_03	AC	1,138 ms; 5,248 KB
testcase_04	AC	1,096 ms; 5,248 KB
testcase_05	AC	1,222 ms; 5,248 KB
testcase_06	AC	920 ms; 5,248 KB
testcase_07	AC	1,261 ms; 5,248 KB
testcase_08	AC	947 ms; 5,248 KB
testcase_09	AC	1,076 ms; 5,248 KB
testcase_10	AC	1,058 ms; 5,248 KB
testcase_11	AC	921 ms; 5,248 KB
testcase_12	AC	1,146 ms; 5,248 KB
testcase_13	AC	1,169 ms; 5,248 KB
testcase_14	AC	1,281 ms; 5,248 KB
testcase_15	AC	1,330 ms; 5,248 KB
testcase_16	AC	1,092 ms; 5,248 KB
testcase_17	AC	1,028 ms; 5,248 KB
testcase_18	AC	1,250 ms; 5,248 KB
testcase_19	AC	1,076 ms; 5,248 KB
testcase_20	AC	2 ms; 5,248 KB

ソースコード

#include <cstdio>

using u32 = unsigned;
using u64 = unsigned long long;
using u128 = __uint128_t;

u32 step(u32 s){
  return (u64) s * 12345 % 1000003;
}

u128 edges[121];

int ctz(const u128 &n){
  if((u64) n){
    return __builtin_ctzll((u64) n);
  }
  return 64 + __builtin_ctzll((u64) (n >> 64));
}

int popc(const u128 &n){
  return __builtin_popcountll((u64) n) + __builtin_popcountll((u64) (n >> 64));
}

u128 mis(const u128 &g){
  if(g == 0) return 0;
  int y = 127;
  int d = 0;
  for(u128 h = g; h; h &= h - 1){
    int z = ctz(h);
    int k = popc(g & edges[z]);
    if(k <= 2) return mis(g & ~edges[z]) | ((u128) 1 << z);
    if(k > d){
      y = z;
      d = k;
    }
  }

  u128 g0 = mis(g & ~edges[y]) | ((u128) 1 << y);

  if(d == 3) return g0;

  u128 g1 = mis(g & ~((u128) 1 << y));

  return popc(g0) > popc(g1) ? g0 : g1;
}



int main(){
  int s;
  scanf("%d", &s);
  s = step(s);
  int n = (s % 120) + 2;
  s = step(s);
  int p = s;
  for(int i = 0; i < n; i++){
    edges[i] |= ((u128) 1 << i);
    for(int j = i+1; j < n; j++){
      s = step(s);
      edges[i] |= ((u128) (s >= p) << j);
      edges[j] |= ((u128) (s >= p) << i);
    }
  }
  u128 g = mis(((u128) 1 << n) - 1);
  int smis = popc(g);
  if(smis == n){
    puts("-1");
    return 0;
  }
  printf("%d\n", popc(g)+1);
  for(u128 h = g; h; ){
    printf("%d", ctz(h) + 1);
    h &= h - 1;
    printf("%c", h ? ' ' : '\n');
  }
  return 0;
}