結果
問題 | No.1358 [Zelkova 2nd Tune *] 語るなら枚数を... |
ユーザー |
![]() |
提出日時 | 2021-01-23 00:38:25 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 543 ms / 2,000 ms |
コード長 | 1,956 bytes |
コンパイル時間 | 2,057 ms |
コンパイル使用メモリ | 201,828 KB |
最終ジャッジ日時 | 2025-01-18 06:57:41 |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 17 |
ソースコード
#include<bits/stdc++.h>using namespace std;#define ALL(x) begin(x),end(x)#define rep(i,n) for(int i=0;i<(n);i++)#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;#define mod 1000000007using ll=long long;const int INF=1000000000;const ll LINF=1001002003004005006ll;int dx[]={1,0,-1,0},dy[]={0,1,0,-1};template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}#pragma GCC target("avx2")#pragma GCC optimize("O3")#pragma GCC optimize("unroll-loops")struct IOSetup{IOSetup(){cin.tie(0);ios::sync_with_stdio(0);cout<<fixed<<setprecision(12);}} iosetup;template<typename T>ostream &operator<<(ostream &os,const vector<T>&v){for(int i=0;i<(int)v.size();i++) os<<v[i]<<(i+1==(int)v.size()?"":" ");return os;}template<typename T>istream &operator>>(istream &is,vector<T>&v){for(T &x:v)is>>x;return is;}template< typename T >inline T extgcd(T a, T b, T &x, T &y) {T d = a;if(b != 0) {d = extgcd(b, a % b, y, x);y -= (a / b) * x;} else {x = 1;y = 0;}return d;}vector<ll> ve(3);inline void solve(){cin>>ve;ll y;cin>>y;sort(ALL(ve));ll a,b,c;a=ve[2],b=ve[1],c=ve[0];ll res=0;ll g=gcd(b,c);ll j,k;ll u=b/g,v=c/g;ll gg=extgcd(u,v,j,k);for(ll i=0;i<=y;i+=a){ll t=y-i;if(t%g) continue;t/=g;ll w=t/gg; // t=g*w// g=j*u+k*v;ll jj=j*w,kk=k*w;// j=j' + x * v// k=k' - x * ujj*=-1;ll dw=(jj<0ll?jj/v:(jj+v-1)/v);ll up=(kk>0ll?kk/u:(kk-u+1)/u);// cout<<up<<" ~ "<<dw<<endl;res+=max(0ll,(up-dw+1));if(res>=mod) res%=mod;}cout<<res<<"\n";}signed main(){int q;cin>>q;while(q--) solve();return 0;}