結果
問題 |
No.990 N×Mマス計算(Kの倍数)
|
ユーザー |
![]() |
提出日時 | 2020-02-15 01:24:45 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,763 ms / 2,000 ms |
コード長 | 1,832 bytes |
コンパイル時間 | 1,600 ms |
コンパイル使用メモリ | 146,896 KB |
最終ジャッジ日時 | 2025-01-09 00:38:49 |
ジャッジサーバーID (参考情報) |
judge3 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 19 |
ソースコード
#include<iostream> #include<string> #include<vector> #include<algorithm> #include<map> #include<set> #include<cstdio> #include<cmath> #include<deque> #include<numeric> #include<queue> #include<stack> #include<cstring> #include<limits> #include<functional> #include<unordered_set> #include<iomanip> #include<cassert> #include<regex> #include<bitset> #include<complex> #include<chrono> #define rep(i,a) for(int i=(int)0;i<(int)a;++i) #define pb push_back #define eb emplace_back using ll=long long; constexpr ll mod = 1e9 + 7; constexpr ll INF = 1LL << 60; template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; } using namespace std; //ここから #define int64_t ll ll gcd(ll n, ll m) { ll tmp; while (m!=0) { tmp = n % m; n = m; m = tmp; } return n; } ll lcm(ll n, ll m) { return abs(n * m) / gcd(n, m);//gl=xy } vector< int64_t > divisor(int64_t n) { vector< int64_t > ret; for(int64_t i = 1; i * i <= n; i++) { if(n % i == 0) { ret.push_back(i); if(i * i != n) ret.push_back(n / i); } } sort(begin(ret), end(ret)); return (ret); } void solve(){ ll n,m,k; cin>>n>>m>>k; char op; cin>>op; vector<ll>b(m),a(n); rep(i,m)cin>>b[i]; rep(i,n)cin>>a[i]; map<ll,ll>mp; rep(i,m)mp[b[i]%k]++; ll ans=0; if(op=='+'){ rep(i,n)ans+=mp[(k-(a[i]%k))%k]; } else { map<ll,ll>mp2; rep(j,m){ ll x=gcd(b[j],k); for(auto y:divisor(x))mp2[y]++; } rep(j,n){ ll x=gcd(a[j],k); ans+=mp2[k/x]; } } cout<<ans<<"\n"; } int main(){ ios::sync_with_stdio(false); cin.tie(0); cout<<fixed<<setprecision(15); solve(); return 0; }