結果
問題 | No.890 移調の限られた旋法 |
ユーザー | ttttan2 |
提出日時 | 2019-09-20 22:02:28 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 155 ms / 2,000 ms |
コード長 | 2,736 bytes |
コンパイル時間 | 1,180 ms |
コンパイル使用メモリ | 149,012 KB |
実行使用メモリ | 20,780 KB |
最終ジャッジ日時 | 2023-10-12 19:10:39 |
合計ジャッジ時間 | 3,911 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 11 ms
7,044 KB |
testcase_01 | AC | 10 ms
7,192 KB |
testcase_02 | AC | 12 ms
7,192 KB |
testcase_03 | AC | 11 ms
7,100 KB |
testcase_04 | AC | 11 ms
7,192 KB |
testcase_05 | AC | 10 ms
7,064 KB |
testcase_06 | AC | 10 ms
7,288 KB |
testcase_07 | AC | 10 ms
7,252 KB |
testcase_08 | AC | 10 ms
7,060 KB |
testcase_09 | AC | 11 ms
7,056 KB |
testcase_10 | AC | 10 ms
7,096 KB |
testcase_11 | AC | 10 ms
7,064 KB |
testcase_12 | AC | 11 ms
7,056 KB |
testcase_13 | AC | 62 ms
20,780 KB |
testcase_14 | AC | 60 ms
20,764 KB |
testcase_15 | AC | 62 ms
20,672 KB |
testcase_16 | AC | 61 ms
20,768 KB |
testcase_17 | AC | 57 ms
19,844 KB |
testcase_18 | AC | 59 ms
19,892 KB |
testcase_19 | AC | 37 ms
17,420 KB |
testcase_20 | AC | 28 ms
15,368 KB |
testcase_21 | AC | 14 ms
9,696 KB |
testcase_22 | AC | 48 ms
19,464 KB |
testcase_23 | AC | 56 ms
20,216 KB |
testcase_24 | AC | 38 ms
17,272 KB |
testcase_25 | AC | 17 ms
10,076 KB |
testcase_26 | AC | 59 ms
20,380 KB |
testcase_27 | AC | 59 ms
20,500 KB |
testcase_28 | AC | 42 ms
19,324 KB |
testcase_29 | AC | 32 ms
15,372 KB |
testcase_30 | AC | 105 ms
20,048 KB |
testcase_31 | AC | 85 ms
15,296 KB |
testcase_32 | AC | 76 ms
19,392 KB |
testcase_33 | AC | 155 ms
19,840 KB |
testcase_34 | AC | 105 ms
19,896 KB |
ソースコード
#include<bits/stdc++.h> //ios::sync_with_stdio(false); //cin.tie(0); using namespace std; typedef long long ll; typedef unsigned long long ull; typedef pair<int,int> pii; typedef pair<pii,int> ppii; typedef pair<int,pii> pipi; typedef pair<ll,ll> pll; typedef pair<pll,ll> ppll; typedef pair<ll,pll> plpl; typedef tuple<ll,ll,ll> tl; ll mod=1000000007; ll mod2=998244353; ll mod3=1000003; ll mod4=998244853; ll inf=1000000000000000000; double pi=2*acos(0); #define rep(i,m,n) for(ll i=m;i<n;i++) #define rrep(i,n,m) for(ll i=n;i>=m;i--) int dh[4]={1,-1,0,0}; int dw[4]={0,0,1,-1}; int ddh[8]={-1,-1,-1,0,0,1,1,1}; int ddw[8]={-1,0,1,-1,1,-1,0,1}; ll lmax(ll a,ll b){ if(a<b)return b; else return a; } ll lmin(ll a,ll b){ if(a<b)return a; else return b; } ll gcd(ll a,ll b){ if(a<b)swap(a,b); if(a%b==0)return b; return gcd(b,a%b); } ll Pow(ll n,ll k){ ll ret=1; ll now=n; while(k>0){ if(k&1)ret*=now; now*=now; k/=2; } return ret; } ll gya[1000010]; ll kai[1000010]; ll beki(ll n,ll k){ ll ret=1; ll now=n; while(k>0){ if(k%2==1){ ret*=now; ret%=mod; } now*=now; now%=mod; k/=2; } return ret; } ll gyaku(ll n){ return beki(n,mod-2); } void nckinit(ll n){ kai[0]=1; kai[1]=1; for(int i=2;i<=n;i++){ kai[i]=kai[i-1]*i; kai[i]%=mod; } gya[n]=gyaku(kai[n]); for(int i=n-1;i>=1;i--){ gya[i]=gya[i+1]*(i+1); gya[i]%=mod; } gya[0]=1; } ll nck(ll n,ll k){ if(k<0)return 0; if(k==0||n==k)return 1; ll ret=kai[n]; ret*=gya[n-k]; ret%=mod; ret*=gya[k]; ret%=mod; return ret; } ll npk(ll n,ll k){ if(k<0)return 0; if(k==0)return 1; ll ret=kai[n]; ret*=gya[n-k]; ret%=mod; return ret; } int main(){ ios::sync_with_stdio(false); cin.tie(0); ll n,k;cin>>n>>k; vector<ll> v; rep(i,2,n+1){ if(n%i==0)v.push_back(i); } nckinit(n+10); ll ans=0; vector<ll> pri; bool used[1000001]; fill(used,used+1000001,false); rep(i,2,1000001){ if(used[i])continue; pri.push_back(i); for(ll j=i;j<=1000000;j+=i)used[j]=true; } rep(i,0,v.size()){ ll now=v[i]; if(k%now!=0)continue; ll u=now; ll cnt=0; rep(j,0,pri.size()){ ll t=pri[j]; if(u%t!=0)continue; if(u%(t*t)==0){ cnt=-1; break; } cnt++; } if(cnt==-1)continue; ll r=nck(n/now,k/now); if(cnt%2==1)ans+=r; else ans-=r; ans=(ans+mod)%mod; } cout<<ans<<endl; }