結果

問題 No.590 Replacement
ユーザー Lepton_sLepton_s
提出日時 2017-11-04 00:53:58
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 2,462 bytes
コンパイル時間 998 ms
コンパイル使用メモリ 95,828 KB
実行使用メモリ 11,392 KB
最終ジャッジ日時 2024-05-02 21:04:53
合計ジャッジ時間 6,503 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 3 ms
5,376 KB
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 TLE -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
testcase_43 -- -
testcase_44 -- -
testcase_45 -- -
testcase_46 -- -
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:61:17: warning: format ‘%d’ expects argument of type ‘int*’, but argument 2 has type ‘ll*’ {aka ‘long long int*’} [-Wformat=]
   61 |         scanf("%d", &n);
      |                ~^   ~~
      |                 |   |
      |                 |   ll* {aka long long int*}
      |                 int*
      |                %lld
main.cpp:61:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   61 |         scanf("%d", &n);
      |         ~~~~~^~~~~~~~~~
main.cpp:63:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   63 |                 scanf("%lld", a[i]+j);
      |                 ~~~~~^~~~~~~~~~~~~~~~

ソースコード

diff # 

// kinen submit
#include <algorithm>
#include <climits>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <ctime>
#include <iostream>
#include <sstream>
#include <functional>
#include <map>
#include <string>
#include <cstring>
#include <vector>
#include <queue>
#include <stack>
#include <deque>
#include <set>
#include <list>
#include <numeric>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll,ll> P;
typedef pair<P,ll> PP;
const double PI = 3.14159265358979323846;
const double EPS = 1e-12;
const ll INF = 1LL<<29;
const ll mod = 1e9+7;
#define rep(i,n) for(int (i)=0;(i)<(ll)(n);++(i))
#define repd(i,n,d) for(ll (i)=0;(i)<(ll)(n);(i)+=(d))
#define all(v) (v).begin(), (v).end()
#define pb(x) push_back(x)
#define mp(x,y) make_pair((x),(y))
#define mset(m,v) memset((m),(v),sizeof(m))
#define chmin(X,Y) ((X)>(Y)?X=(Y),true:false)
#define chmax(X,Y) ((X)<(Y)?X=(Y),true:false)
#define fst first
#define snd second
#define UNIQUE(x) (x).erase(unique(all(x)),(x).end())
template<class T> ostream &operator<<(ostream &os, const vector<T> &v){int n=v.size();rep(i,n)os<<v[i]<<(i==n-1?"":" ");return os;}

#define N 100100
ll n;
ll a[2][N], c[2][N], p[2][N], pr[2][N], used[N], ix[N], lp[N], dd[N], rr[N];
map<PP, vector<ll>> s;

void extgcd(ll x, ll y, ll &a, ll &b, ll &g){
	if(y == 0){
		a = 1;
		b = 0;
		g = x;
	} else {
		extgcd(y, x%y, a, b, g);
		swap(a, b);
		b -= a*x/y;
	}
}

int main(){
	scanf("%d", &n);
	rep(i, 2) rep(j, n){
		scanf("%lld", a[i]+j);
		a[i][j]--;
	}
	mset(p, -1);
	rep(i, 2) rep(j, n){
		if(p[i][j] != -1) continue;
		ll k = a[i][j];
		p[i][j] = 1;
		while(k != j){
			c[i][k] = p[i][j]++;
			pr[i][k] = j;
			k = a[i][k];
		}
		pr[i][j] = j;
	}
	ll res = 0;
	rep(j, n){
		ll b[2], d;
		extgcd(p[0][j], p[1][j], b[0], b[1], d);
		ll r[2], q[2];
		rep(i, 2) r[i] = c[i][j]%d;
		if(r[1]<r[0]) r[1] += d;
		rep(i, 2) q[i] = c[i][j]/d;
		ll qq = q[1]-q[0];
		b[0] *= qq; b[0] %= p[1][j]/d; if(b[0]<0) b[0] += p[1][j]/d;
		b[1] * -qq; b[1] %= p[0][j]/d; if(b[1]<0) b[1] += p[0][j]/d;
		ix[j] = b[0]*p[0][j]*d+r[0];
		lp[j] = p[0][j]/d*p[1][j]/d;
		dd[j] = d;
		rr[j] = r[1]-r[0];
		s[PP(P(pr[0][j], pr[1][j]), r[j])].push_back(j);
	}
	for(auto &&x: s){
		auto &&v = x.snd;
		int m = v.size();
		vector<ll> vv(m+1);
		rep(i, m) vv[i] = ix[v[i]];
		v[m] = vv[0]+lp[v[0]];
		rep(i, m){
			ll y = vv[i+1]-vv[i]-1;
			res += y*(y+1)/2;
			res %= mod;
		}
	}
	cout<<res<<endl;
	return 0;
}
0