結果

問題 No.415 ぴょん
ユーザー uenoku
提出日時 2016-08-26 22:57:13
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 3,428 bytes
コンパイル時間 1,273 ms
コンパイル使用メモリ 101,672 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-08 09:42:14
合計ジャッジ時間 2,146 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 15 WA * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
#include <cstdio>
#include <queue>
#include "math.h"
#include <complex>
#include <iomanip>
#include <map>
#define ifor(i,a,b) for (int i=(a);i<(b);i++)
#define rfor(i,a,b) for (int i=(b)-1;i>=(a);i--)
#define rep(i,n) for (int i=0;i<(n);i++)
#define rrep(i,n) for (int i=(n)-1;i>=0;i--)
using namespace std;
typedef long double ld;
typedef long long int  lli;
typedef complex <double> P;
const double eps = 1e-11;
int vex[4]={1,0,-1,0};
int vey[4]={0,1,0,-1};
typedef vector<double> Vec;
typedef vector<int> vec;
typedef vector<Vec> MAT;
typedef vector<vec> mat;
lli MOD=1000000007;
//Ax=bをとくAは正方行列
//rankA<=min(m,n)ならば配列０のvecが帰る
Vec gauss_jordan(const MAT& A,const Vec&b){
	int n  = A.size();
	MAT B (n,Vec(n+1));
	rep(i,n)rep(j,n)B[i][j]=A[i][j];
	rep(i,n)B[i][n]=b[i];

	rep(i,n){
		int pivot = i;
		for(int j= i;j<n;j++){
			if(abs(B[j][i])>abs(B[pivot][i]))pivot = j;
		}

		swap(B[i],B[pivot]);
		if(abs(B[i][i])<eps)return Vec();
    //B_i_i成分が0であるつまり階数が
		for(int j =i +1 ;j<=n;j++)B[i][j]/=B[i][i];
		rep(j,n){
			if(i!=j)
				for(int k = i+1;k<=n;k++){
					B[j][k] -=B[j][i]*B[i][k];
				}
		}
	}
	Vec x(n);
	for(int i =0;i<n;i++){
		x[i]= B[i][n];
	}
	return x;


}
double det(const MAT& A){
	int n = A.size();
	MAT B (n,Vec(n));
	rep(i,n)rep(j,n)B[i][j]=A[i][j];
	double ans = 1;
	rep(i,n){
		int pivot = i;
		for(int j=i;j<n;j++){
			if(abs(B[j][i])>abs(B[pivot][i]))pivot = j;
		}
		if(i!=pivot)ans *= -1;
		swap(B[i],B[pivot]);
		if(abs(B[i][i])<eps)
			return 0;
		ans *= B[i][i];
		for(int j =i +1 ;j<n;j++)B[i][j]/=B[i][i];
		rep(j,n){
			if(i!=j)
				for(int k = i+1;k<n;k++){
					B[j][k] -=B[j][i]*B[i][k];
				}
		}
	}
	return ans;

}
int  rank(const MAT& A){
	int n = A.size();
	MAT B (n,Vec(n));
	rep(i,n)rep(j,n)B[i][j]=A[i][j];

	rep(i,n){
		int pivot = i;
		for(int j=i;j<n;j++){
			if(abs(B[j][i])>abs(B[pivot][i]))pivot = j;
		}
		swap(B[i],B[pivot]);
		if(abs(B[i][i])<eps)return i;
		for(int j =i +1 ;j<n;j++)B[i][j]/=B[i][i];
		rep(j,n){
			if(i!=j)
				for(int k = i+1;k<n;k++){
					B[j][k] -=B[j][i]*B[i][k];
				}
		}
	}
	return n;
}
lli prime(lli m){
  for(lli i =2;i*i<=m;i++){
    if(m%i==0)return i;
  }
  return m;
}
vector<lli> euler(lli m){
	vector<lli> p;
	lli M = m;
	for(int i = 2;i<=m;i++){
		if(m%i==0){
			p.push_back(i);
      return p;
			while(m%i==0)m/=i;
		}
		if(M<i*i&&p.size()==0){
			p.push_back(M);
			break;
		}
	}
	lli ans = M;
  return p;
}
lli powm(lli a,lli p,lli mod){
	lli ans =1;
	while(p>0){
		if(p&1)ans=(ans*a) % mod;
		a=(a*a)%mod;
		p >>=1;
	}
	return ans%mod;
}

lli inv(lli a){
	return  powm(a,MOD-2,MOD);
}
lli gcd(lli A,lli B){
  if(B>A)swap(A,B);
	if(A%B==0)return B;
	else return gcd(B,A%B);
}
mat mul_mat_mod(mat A,mat B,lli m){
	int n = A.size();
	mat C (n,vec(n));
	rep(i,n)rep(j,n)rep(k,n){
		C[i][j]+=A[i][k]*B[k][j] %m;
		C[i][j]%=m;
	}
	return C;
}
mat pow_mat(mat A,lli p,lli mod){
	int n= A.size();
	mat B = mat(n,vec(n));
	while(p>0){
		if(p&1){
			B=mul_mat_mod(A,B,mod);
		}
		A = mul_mat_mod(A,A,mod);
		p>>=1;
	}
	return B;
}
lli comb(lli a,lli b ){
	lli ans =1;
	rep(i,b){
		ans = ans %MOD* (a-i)%MOD*inv(b-i)%MOD;
	}
	return ans ;
}
int main(){
  lli N,D;
  cin >> N >> D;
  lli g = gcd(N,D);
  if(g==1)cout << N-1<<endl;
  else if(D%N==0){cout << 0<<endl;}
  else{
    cout << N/g<<endl;
  }

}
0