結果
問題 | No.890 移調の限られた旋法 |
ユーザー | ttttan2 |
提出日時 | 2019-09-20 22:02:28 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 156 ms / 2,000 ms |
コード長 | 2,736 bytes |
コンパイル時間 | 1,449 ms |
コンパイル使用メモリ | 164,232 KB |
実行使用メモリ | 20,920 KB |
最終ジャッジ日時 | 2024-09-14 17:35:01 |
合計ジャッジ時間 | 3,997 ms |
ジャッジサーバーID (参考情報) |
judge6 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 10 ms
5,248 KB |
testcase_01 | AC | 9 ms
5,376 KB |
testcase_02 | AC | 12 ms
5,376 KB |
testcase_03 | AC | 11 ms
5,376 KB |
testcase_04 | AC | 10 ms
5,376 KB |
testcase_05 | AC | 10 ms
5,376 KB |
testcase_06 | AC | 10 ms
5,376 KB |
testcase_07 | AC | 10 ms
5,376 KB |
testcase_08 | AC | 10 ms
5,376 KB |
testcase_09 | AC | 11 ms
5,376 KB |
testcase_10 | AC | 10 ms
5,376 KB |
testcase_11 | AC | 10 ms
5,376 KB |
testcase_12 | AC | 11 ms
5,428 KB |
testcase_13 | AC | 62 ms
20,920 KB |
testcase_14 | AC | 60 ms
20,920 KB |
testcase_15 | AC | 62 ms
20,920 KB |
testcase_16 | AC | 62 ms
20,796 KB |
testcase_17 | AC | 57 ms
19,236 KB |
testcase_18 | AC | 60 ms
19,364 KB |
testcase_19 | AC | 37 ms
13,496 KB |
testcase_20 | AC | 27 ms
10,528 KB |
testcase_21 | AC | 14 ms
6,404 KB |
testcase_22 | AC | 51 ms
17,388 KB |
testcase_23 | AC | 59 ms
19,972 KB |
testcase_24 | AC | 37 ms
13,688 KB |
testcase_25 | AC | 16 ms
7,264 KB |
testcase_26 | AC | 61 ms
20,280 KB |
testcase_27 | AC | 59 ms
20,536 KB |
testcase_28 | AC | 42 ms
15,076 KB |
testcase_29 | AC | 32 ms
11,564 KB |
testcase_30 | AC | 105 ms
19,372 KB |
testcase_31 | AC | 87 ms
13,272 KB |
testcase_32 | AC | 76 ms
18,292 KB |
testcase_33 | AC | 156 ms
19,204 KB |
testcase_34 | AC | 105 ms
19,208 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; }