답안 #20110

#	제출 시각	아이디	문제	언어	결과	실행 시간	메모리
20110	2016-02-25T16:07:53 Z	gs14004	팔찌 (kriii4_V)	C++14	6 / 100	1000 ms	13444 KB

V

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <limits.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <algorithm>
#include <string>
#include <functional>
#include <vector>
#include <numeric>
#include <deque>
#include <utility>
#include <bitset>
#include <iostream>
using namespace std;
typedef long long lint;
typedef long double llf;
typedef pair<int, int> pi;
const int mod = 1e9 + 7;

lint ipow(lint x, lint p){
    lint ret = 1, piv = x % mod;
    while(p){
        if(p&1) ret *= piv;
        piv *= piv;
        ret %= mod;
        piv %= mod;
        p >>= 1;
    }
    return ret;
}

int dp[1000005];

int phi(int x){
	if(~dp[x]) return dp[x];
	int tmp = x;
	lint ret = x;
	set<int> s;
	int cnt = 0;
	for(int i=2; i*i<=x; i++){
		if(x % i == 0) s.insert(i);
		while(x % i == 0){
			x /= i;
		}
	}
	if(x > 1) s.insert(x);
	for(auto &i : s){
		ret *= (i-1);
		ret /= i;
	}
	return dp[tmp] = ret;
}

int n, k;

int solve(int n, int k){
	lint ret = 0;
	for(int i=1; i*i<=n; i++){
		if(n%i == 0){
			ret += ipow(k, i) * phi(n/i) % mod;
			if(i*i!=n) ret += ipow(k, n/i) * phi(i) % mod;
			ret %= mod;
		}
	}
	ret *= ipow(n, mod - 2);
	ret %= mod;
	return ret;
}

lint kpow[1000005];

int main(){
	memset(dp, -1, sizeof(dp));
	cin >> n >> k;
	kpow[0] = 1;
	for(int i=1; i<=n; i++){
		kpow[i] = kpow[i-1] * k % mod;
	}
	lint ret = 2;
	for(int i=1; i<=n; i++){
		ret += solve(i, k);
		ret %= mod;
		if(i%2 == 0) ret += (1ll * (k+1) * kpow[i/2] % mod) * ((mod + 1) / 2) % mod;
		else ret += kpow[i/2+1];
		ret %= mod;
	}
	ret *= (mod + 1) / 2;
	ret %= mod;
	cout << ret;
}

Subtask #1

6.0 / 6.0

#	결과	실행 시간	메모리	Grader output
1	Correct	0 ms	13444 KB	Output is correct
2	Correct	0 ms	13444 KB	Output is correct
3	Correct	0 ms	13444 KB	Output is correct
4	Correct	0 ms	13444 KB	Output is correct
5	Correct	0 ms	13444 KB	Output is correct
6	Correct	0 ms	13444 KB	Output is correct
7	Correct	1 ms	13444 KB	Output is correct
8	Correct	0 ms	13444 KB	Output is correct
9	Correct	0 ms	13444 KB	Output is correct
10	Correct	0 ms	13444 KB	Output is correct
11	Correct	0 ms	13444 KB	Output is correct
12	Correct	1 ms	13444 KB	Output is correct
13	Correct	0 ms	13444 KB	Output is correct
14	Correct	0 ms	13444 KB	Output is correct
15	Correct	0 ms	13444 KB	Output is correct
16	Correct	0 ms	13444 KB	Output is correct
17	Correct	0 ms	13444 KB	Output is correct
18	Correct	0 ms	13444 KB	Output is correct
19	Correct	1 ms	13444 KB	Output is correct
20	Correct	0 ms	13444 KB	Output is correct

Subtask #2

0.0 / 94.0

#	결과	실행 시간	메모리	Grader output
1	Correct	0 ms	13444 KB	Output is correct
2	Correct	0 ms	13444 KB	Output is correct
3	Correct	0 ms	13444 KB	Output is correct
4	Correct	172 ms	13444 KB	Output is correct
5	Execution timed out	1000 ms	13440 KB	Program timed out
6	Halted	0 ms	0 KB	-