結果
問題 | No.1463 Hungry Kanten |
ユーザー | snow39 |
提出日時 | 2021-04-02 21:58:44 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 955 ms / 2,000 ms |
コード長 | 1,800 bytes |
コンパイル時間 | 2,277 ms |
コンパイル使用メモリ | 125,584 KB |
最終ジャッジ日時 | 2025-01-20 09:13:39 |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 20 |
ソースコード
#include <iostream> #include <algorithm> #include <string> #include <vector> #include <cmath> #include <map> #include <queue> #include <iomanip> #include <set> #include <tuple> #define mkp make_pair #define mkt make_tuple #define rep(i,n) for(int i = 0; i < (n); ++i) #define all(v) v.begin(),v.end() using namespace std; typedef long long ll; const ll MOD=1e9+7; template<class T> void chmin(T &a,const T &b){if(a>b) a=b;} template<class T> void chmax(T &a,const T &b){if(a<b) a=b;} const int L=1e3; const ll INF=1e9; int main(){ cin.tie(0); ios::sync_with_stdio(false); int N,K; cin>>N>>K; vector<int> A(N); rep(i,N) cin>>A[i]; sort(all(A)); vector<int> isPrime(L+1,1); isPrime[0]=0;isPrime[1]=0; for(int i=2;i<=L;i++){ if(isPrime[i]==0) continue; for(int j=2;i*j<=L;j++) isPrime[i*j]=0; } vector<int> primes; for(int i=2;i<=L;i++) if(isPrime[i]) primes.push_back(i); int V=primes.size(); vector<vector<int>> pcnt(N); rep(i,N){ rep(j,V){ int value=A[i]; while(value%primes[j]==0){ pcnt[i].push_back(primes[j]); value/=primes[j]; } } } map<int,int> mp; map<vector<int>,int> mulmp; for(int bit=0;bit<(1<<N);bit++){ int sz=__builtin_popcount(bit); if(sz<K) continue; vector<int> v; rep(i,N) if(bit>>i&1) v.push_back(i); int sum=0; for(auto i:v) sum+=A[i]; mp[sum]=1; ll value=1; bool over=false; for(auto i:v){ value*=A[i]; if(value>INF){ over=true; break; } } if(over==false){ mp[value]=1; continue; } vector<int> ps; for(auto i:v){ for(auto n:pcnt[i]) ps.push_back(n); } sort(all(ps)); mulmp[ps]=1; } int ans=mp.size()+mulmp.size(); cout<<ans<<endl; return 0; }