結果

問題 No.1189 Sum is XOR
ユーザー Rho
提出日時 2020-08-02 00:44:44
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 143 ms / 2,000 ms
コード長 1,298 bytes
コンパイル時間 1,613 ms
コンパイル使用メモリ 169,372 KB
実行使用メモリ 12,940 KB
最終ジャッジ日時 2024-07-16 01:00:04
合計ジャッジ時間 5,309 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
#define int long long
#define rep(i,n) for(int i=0;i<n;i++)
#define vi vector<int>
#define all(a) a.begin(),a.end()
typedef pair<int,int> P;
const long long mod=998244353;
const long long inf=1ll<<61;
const int maxN=100;
int cnt[1111];
int cmb[1111];
int dp[12][1024];int kj[600006],kji[600006];
int modpow(int a,int x){
  int res=1;
  while(x){
    if(x&1)res=res*a%mod;
    a=a*a%mod;
    x>>=1;
  }
  return res;
}
void setkj(){
  kj[0]=kji[0]=1;
  for(int i=1;i<=600002;i++){
    kj[i]=kj[i-1]*i%mod;
    kji[i]=modpow(kj[i],mod-2);
  }
}
int comb(int r,int c){
  if(r<c)return 0;
  return kj[r]*kji[c]%mod*kji[r-c]%mod;
}
signed main(){
    setkj();
    int n,k;cin>>n>>k;
    rep(i,n){
        int a;cin>>a;cnt[a]++;
    }
    dp[0][0]=1;
    rep(j,11){
        rep(l,1<<10){
            if(!dp[j][l])continue;
            int m=1023-l;
            for(int t=m;;t=(t-1)&m){
                if(t==0)break;
                dp[j+1][t+l]+=dp[j][l]*cnt[t]%mod;
                dp[j+1][t+l]%=mod;
            }
        }
    }

    int sum=0;
    rep(i,1<<10){
        rep(j,min(k+1,11ll)){
            dp[j][i]*=kji[j];dp[j][i]%=mod;
            sum+=dp[j][i]*comb(cnt[0],k-j)%mod;
            sum%=mod;
        }
    }
    
    cout<<sum%mod<<endl;
}
0