結果

問題 No.2755 行列の共役類
ユーザー ttkkggwwttkkggww
提出日時 2024-05-24 01:55:45
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 265 ms / 3,500 ms
コード長 2,419 bytes
コンパイル時間 5,189 ms
コンパイル使用メモリ 263,568 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-20 19:12:04
合計ジャッジ時間 7,834 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 66
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <algorithm>
#include <bits/stdc++.h>
using namespace std;
#include <atcoder/all>
using namespace atcoder;
using ll = long long;
int B,C;
vector<pair<int,int>> units[1000];
using mint = modint;
/*
bool judge(int x,int n_1,int y,int n_2,bool test = false){
for(int z = 1;z<=B;z++){
if(gcd(z,B)!=1)
continue;
if(x==1){
if(mint(n_1-n_2*z).val()==0)
return true;
continue;
}
mint bunshi = mint(n_1 - n_2*z);
mint bunbo = mint(1-x);
if(test){
cerr<<n_2<<' '<<z<<' '<<bunshi.val()<<" "<<bunbo.val()<<endl;
}
if(bunshi.val()==0)
return true;
int g = gcd(bunshi.val(),bunbo.val());
int tmp = bunbo.val()/g;
if(gcd(tmp,B/g)==1){
return true;
}
}
return false;
}
*/
bool judge(int x,int n_1,int y,int n_2,bool test = false){
for(int z = 0;z<=B;z++){
if(gcd(z,B)!=1)
continue;
for(int n_3 = 0;n_3<B/C;n_3++){
mint tmp = mint(n_3*x+n_1- (n_2*z + n_3));
if(tmp.val()==0){
if(test){
cout<<x<<" "<<n_1<<" "<<y<<" "<<n_2<<" "<<z<<" "<<n_3<<endl;
}
return true;
}
}
}
return false;
}
void test(int x,int n_1){
for(auto &[y,n_2]:units[x]){
if(n_1==n_2)
continue;
if(judge(x,n_1,y,n_2,true)){
cout<<"error"<<endl;
return;
}
}
}
void solve(){
int cnt = 0;
for(int n_1 = 0;n_1<B/C;n_1++){
for(int x = 1;x<=B;x++){
if(gcd(x,B)!=1)
continue;
bool is_add = true;
for(auto &[y,n_2]:units[x]){
if(judge(x,n_1,y,n_2)){
is_add = false;
break;
}
if(judge(y,n_2,x,n_1)){
is_add = false;
break;
}
}
if(is_add){
//cerr<<x<<" "<<n_1<<endl;
units[x].emplace_back(x,n_1);
cnt++;
}
}
if(cnt>100){
cout<<"100+"<<endl;
return;
}
}
cout<<cnt<<endl;
}
signed main(){
ios::sync_with_stdio(false);
cin.tie(nullptr);
cin >> B >> C;
mint::set_mod(B/C);
solve();
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0