結果

問題 No.590 Replacement
コンテスト
ユーザー ryoissy
提出日時 2017-11-06 15:21:27
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 88 ms / 2,000 ms
コード長 1,608 bytes
コンパイル時間 1,473 ms
コンパイル使用メモリ 179,848 KB
実行使用メモリ 17,544 KB
最終ジャッジ日時 2024-11-24 03:08:12
合計ジャッジ時間 5,843 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 47
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘void f(std::vector<int>&, std::vector<int>&, std::vector<int>&, std::vector<int>&)’:
main.cpp:20:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   20 |                 scanf("%d",&a[i]);
      |                 ~~~~~^~~~~~~~~~~~
main.cpp: In function ‘int main()’:
main.cpp:43:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   43 |         scanf("%d",&n);
      |         ~~~~~^~~~~~~~~

ソースコード

diff # 

#include <bits/stdc++.h>
#define MOD 1000000007LL
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
typedef pair<P,int> PP;
int n;

ll gcd(ll a,ll b){
	if(b==0LL)return a;
	return gcd(b,a%b);
}

ll lcm(ll a,ll b){
	return a/gcd(a,b)*b;
}

void f(vector<int>& a,vector<int>& c,vector<int>& d,vector<int>& s){
	for(int i=0;i<n;i++){
		scanf("%d",&a[i]);
		a[i]--;
	}
	int cnt=0;
	for(int i=0;i<n;i++){
		if(c[i]==-1){
			int j=i,num=0;
			do{
				c[j]=cnt;
				d[j]=num;
				num++;
				j=a[j];
			}while(j!=i);
			s.push_back(num);
			cnt++;
		}
	}
} 

vector<int> a,b,a_gr,a_ord,a_size,b_gr,b_ord,b_size;
map<PP,vector<P> > mp;

int main(void){
	scanf("%d",&n);
	a=b=a_ord=b_ord=vector<int>(n,0);
	a_gr=b_gr=vector<int>(n,-1);
	f(a,a_gr,a_ord,a_size);
	f(b,b_gr,b_ord,b_size);
	ll ans=0;
	for(int i=0;i<n;i++){
		int g=gcd(a_size[a_gr[i]],b_size[b_gr[i]]);
		mp[PP(P(a_gr[i],b_gr[i]),(a_ord[i]%g-b_ord[i]%g+g)%g)].push_back(P(a_ord[i],b_ord[i]));
	}
	map<PP,vector<P> >::iterator it=mp.begin();
	for(;it!=mp.end();it++){
		P tmp=it->first.first;
		vector<P> vp=it->second;
		int A=a_size[tmp.first],B=b_size[tmp.second],S=vp.size();
		ll l=lcm(A,B),g=gcd(A,B),P=vp[0].first,Q=vp[0].second;
		vector<ll> u(S+1),w(B/g);
		u[S]=l;
		int j=Q/g;
		for(int i=0;i<B/g;i++){
			(j+=A/g)%=B/g;
			w[j]=B/g-i-1;
		}
		for(int i=1;i<S;i++){
			ll X=vp[i].first,Y=vp[i].second;
			ll num=(P+A-X)%A;
			(Y+=num)%=B;
			num+=w[Y/g]*A;
			u[i]=num;
		}
		sort(u.begin(),u.end());
		for(int i=0;i<S;i++){
			ll sa=(u[i+1]-u[i])%MOD;
			(ans+=sa*(sa-1LL)/2LL)%=MOD;
		}
	}
	printf("%lld\n",ans);
	return 0;
}
0