結果
問題 | No.1262 グラフを作ろう! |
ユーザー |
|
提出日時 | 2020-10-16 23:22:40 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 765 ms / 3,000 ms |
コード長 | 1,888 bytes |
コンパイル時間 | 2,369 ms |
コンパイル使用メモリ | 210,196 KB |
最終ジャッジ日時 | 2025-01-15 09:46:35 |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 96 |
コンパイルメッセージ
main.cpp: In function ‘int main()’: main.cpp:86:23: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result] 86 | int n,m; scanf("%d%d",&n,&m); | ~~~~~^~~~~~~~~~~~~~ main.cpp:88:23: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result] 88 | rep(i,m) scanf("%d",&a[i]); | ~~~~~^~~~~~~~~~~~
ソースコード
#include <bits/stdc++.h> #define rep(i,n) for(int i=0;i<(n);i++) using namespace std; using lint=long long; class linear_sieve{ vector<int> lpf,p; public: linear_sieve(int n):lpf(n+1){ for(int i=2;i<=n;i++){ if(lpf[i]==0){ lpf[i]=i; p.emplace_back(i); } for(int j=0;j<p.size()&&p[j]<=lpf[i]&&i*p[j]<=n;j++) lpf[i*p[j]]=p[j]; } } const vector<int>& primes()const{ return p; } bool is_prime(int a)const{ assert(a<=(int)lpf.size()-1); return a>0 && lpf[a]==a; } map<int,int> prime_factorize(int a)const{ assert(a<=(int)lpf.size()-1); map<int,int> pf; for(;a>1;a/=lpf[a]) ++pf[lpf[a]]; return pf; } int number_of_divisors(int a)const{ assert(a<=(int)lpf.size()-1); int res=1,cnt=0,pre=-1; for(;a>1;a/=lpf[a]){ if(pre==-1 || pre==lpf[a]){ cnt++; } else{ res*=cnt+1; cnt=1; } pre=lpf[a]; } return res*(cnt+1); } }; template<class T> vector<T> divisors(const map<T,int>& pf){ vector<T> res={T(1)}; for(const auto& q:pf){ int m=res.size(); T pp=1; rep(i,q.second){ pp*=q.first; rep(i,m) res.emplace_back(res[i]*pp); } } sort(res.begin(),res.end()); return res; } template<class T> vector<pair<T,int>> run_length_encoding(const vector<T>& a){ vector<pair<T,int>> res; int n=a.size(),pre=0; rep(i,n) if(i==n-1 || a[i]!=a[i+1]) res.emplace_back(a[i],i-pre+1), pre=i+1; return res; } int phi[1000001]; void build_phi(){ rep(i,1000001) phi[i]=i; for(int i=2;i<=1000000;i++) if(phi[i]==i) { for(int j=i;j<=1000000;j+=i) phi[j]=phi[j]/i*(i-1); } } int main(){ build_phi(); int n,m; scanf("%d%d",&n,&m); vector<int> a(m); rep(i,m) scanf("%d",&a[i]); sort(a.begin(),a.end()); lint ans=0; linear_sieve LS(1e6); for(auto [x,c]:run_length_encoding(a)){ for(int d:divisors(LS.prime_factorize(x))){ ans+=lint(c)*phi[x/d]*((x-1)%d); } } printf("%lld\n",ans); return 0; }