結果

問題 No.1411 Hundreds of Conditions Sequences
ユーザー 沙耶花沙耶花
提出日時 2021-02-27 00:12:33
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 239 ms / 2,000 ms
コード長 1,274 bytes
コンパイル時間 4,745 ms
コンパイル使用メモリ 269,760 KB
実行使用メモリ 16,412 KB
最終ジャッジ日時 2024-10-02 16:36:24
合計ジャッジ時間 19,410 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 25 ms
15,512 KB
testcase_01 AC 27 ms
15,500 KB
testcase_02 AC 25 ms
15,664 KB
testcase_03 AC 26 ms
15,628 KB
testcase_04 AC 27 ms
15,368 KB
testcase_05 AC 25 ms
15,500 KB
testcase_06 AC 27 ms
15,488 KB
testcase_07 AC 26 ms
15,456 KB
testcase_08 AC 26 ms
15,496 KB
testcase_09 AC 27 ms
15,584 KB
testcase_10 AC 25 ms
15,568 KB
testcase_11 AC 26 ms
15,452 KB
testcase_12 AC 206 ms
16,144 KB
testcase_13 AC 160 ms
15,992 KB
testcase_14 AC 95 ms
15,628 KB
testcase_15 AC 155 ms
15,928 KB
testcase_16 AC 35 ms
15,500 KB
testcase_17 AC 151 ms
16,008 KB
testcase_18 AC 205 ms
16,268 KB
testcase_19 AC 224 ms
16,272 KB
testcase_20 AC 97 ms
15,632 KB
testcase_21 AC 234 ms
16,396 KB
testcase_22 AC 102 ms
15,632 KB
testcase_23 AC 205 ms
16,116 KB
testcase_24 AC 25 ms
15,456 KB
testcase_25 AC 204 ms
16,268 KB
testcase_26 AC 152 ms
16,008 KB
testcase_27 AC 98 ms
15,756 KB
testcase_28 AC 178 ms
15,992 KB
testcase_29 AC 45 ms
15,516 KB
testcase_30 AC 134 ms
15,884 KB
testcase_31 AC 126 ms
15,756 KB
testcase_32 AC 236 ms
16,400 KB
testcase_33 AC 235 ms
16,268 KB
testcase_34 AC 233 ms
16,268 KB
testcase_35 AC 235 ms
16,392 KB
testcase_36 AC 236 ms
16,268 KB
testcase_37 AC 236 ms
16,264 KB
testcase_38 AC 233 ms
16,236 KB
testcase_39 AC 236 ms
16,268 KB
testcase_40 AC 231 ms
16,396 KB
testcase_41 AC 231 ms
16,344 KB
testcase_42 AC 235 ms
16,300 KB
testcase_43 AC 233 ms
16,268 KB
testcase_44 AC 239 ms
16,268 KB
testcase_45 AC 232 ms
16,232 KB
testcase_46 AC 231 ms
16,412 KB
testcase_47 AC 237 ms
16,272 KB
testcase_48 AC 236 ms
16,268 KB
testcase_49 AC 231 ms
16,140 KB
testcase_50 AC 236 ms
16,264 KB
testcase_51 AC 236 ms
16,264 KB
testcase_52 AC 220 ms
16,400 KB
testcase_53 AC 222 ms
16,268 KB
testcase_54 AC 222 ms
16,268 KB
testcase_55 AC 221 ms
16,228 KB
testcase_56 AC 222 ms
16,272 KB
testcase_57 AC 221 ms
16,400 KB
testcase_58 AC 220 ms
16,228 KB
testcase_59 AC 224 ms
16,396 KB
testcase_60 AC 223 ms
16,396 KB
testcase_61 AC 219 ms
16,352 KB
testcase_62 AC 201 ms
16,272 KB
testcase_63 AC 198 ms
16,392 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <atcoder/all>
#include <bits/stdc++.h>
using namespace std;
using namespace atcoder;
using mint = modint1000000007;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf 1000000001


int main(){
	
	vector<int> isp(1000001,-1);
	vector<int> p;
	for(int i=2;i<isp.size();i++){
		if(isp[i]==-1){
			p.push_back(i);
			for(int j=i;j<isp.size();j+=i)isp[j] = i;
		}
	}
	int N;
	cin>>N;
	vector<int> A(N);
	mint sum = 1;
	vector<int> cnt(1000001,0);
	rep(i,N){
		scanf("%d",&A[i]);
		cnt[A[i]]++;
		sum *= A[i];
	}
	vector<int> cnt2(1000001,0);
	rep(i,p.size()){
		long long cur = p[i];
		while(cur <= 1000000){
			
			for(int j=cur;j<=1000000;j+=cur){
				cnt2[cur] += cnt[j];
			}
			if(cnt2[cur]>=2){
				sum /= mint(p[i]).pow(cnt2[cur]-1);
			}
			cur *= p[i];
		}
	}
	vector<mint> ans(N,0);
	//cout<<sum.val()<<endl;
	rep(i,N){
		mint temp = sum / A[i];
		map<int,int> mp;
		int t = A[i];
		while(t!=1){
			mp[isp[t]]++;
			t/=isp[t];
		}
		for(auto a:mp){
			int x = 1;
			rep(j,a.second){
				x *= a.first;
				if(cnt2[x]>=2){
					temp *= a.first;
				}
			}
		}
		ans[i] = temp;
	}
	
	sum = 1;
	rep(i,N){
		sum *= A[i];
	}
	
	rep(i,N){
		ans[i] = sum / A[i] - ans[i];
	}
	
	rep(i,N){
		cout<<ans[i].val()<<endl;
	}
	
	
		
    return 0;
}
0