結果

問題 No.1871 divisXor
コンテスト
ユーザー pockyny
提出日時 2022-09-20 11:22:12
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 45 ms / 2,000 ms
コード長 1,278 bytes
コンパイル時間 903 ms
コンパイル使用メモリ 73,964 KB
最終ジャッジ日時 2025-02-07 13:00:57
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 29
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>

using namespace std;
typedef long long ll;
int a[8] = {2,3,5,7,13,17,19,31};
int sum[110],d[2010];
ll pw(ll a,ll x){
    ll ret = 1;
    while(x){
        if(x&1) (ret *= a);
        (a *= a); x /= 2;
    }
    return ret;
}
int main(){
    int i,j,k;
    for(i=0;i<(1<<8);i++){
        int s = 0;
        for(j=0;j<8;j++){
            if(i>>j&1) s += a[j];
        }
        sum[s] = i;
    }
    for(i=1;i<=2000;i++){
        for(j=i;j<=2000;j+=i) d[j] += i;
    }
    ll n;
    cin >> n;
    if(n==0){
        cout << -1 << "\n";
        return 0;
    }
    vector<ll> v;
    for(j=40;j>=0;j--){
        if(n>>j&1){
            if(sum[j]){
                ll pro = 1;
                for(k=0;k<8;k++){
                    if(sum[j]>>k&1) pro *= pw(2,a[k]) - 1;
                }
                n ^= (1LL<<j);
                v.push_back(pro);
            }else{
                ll x = 1LL<<j;
                for(k=1;k<=2000;k++){
                    if(d[k]>=x && d[k]<2*x){
                        v.push_back(k);
                        n ^= d[k];
                        break;
                    }
                }
            }
        }
    }
    cout << v.size() << "\n";
    for(ll x:v) cout << x << " ";
    cout << "\n";
}
0