답안 #19816

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
19816 2016-02-25T05:54:34 Z tonyjjw 팔찌 (kriii4_V) C++14
6 / 100
1000 ms 101452 KB
//*
#include <stdio.h>
#include <string.h>
#include <ctype.h>
#include <time.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <numeric>
#include <functional>
#define MOD 1000000007
#define MAX 0x3f3f3f3f
#define MAX2 0x3f3f3f3f3f3f3f3fll
#define ERR 1e-10
#define mp make_pair
#define all(x) (x).begin(), (x).end()
#pragma warning(disable:4996)
using namespace std;

typedef long long ll;
typedef long double ldb;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef pair<double, double> pdd;

int n, m;
vector<int> d[1040000];
ll pi[1040000];
ll powm[1040000];
int comp[1040000];

ll mul_inv(ll a, ll b=MOD)
{
	ll b0=b, t, q;
	ll x0=0, x1=1;
	if(b == 1) return 1;
	while(a > 1)
	{
		q=a/b;
		t=b, b=a%b, a=t;
		t=x0, x0=x1-q*x0, x1=t;
	}
	if(x1 < 0) x1+=b0;
	return x1;
}


ll modpow(ll x, ll y)
{
	ll ans=1;
	while(y)
	{
		if(y%2) ans*=x, ans%=MOD;
		x*=x, x%=MOD, y/=2;
	}
	return ans;
}

int main()
{
	int i, j, k;
	cin>>n>>m;
	for(i=1;i<=n;i++) pi[i]=i;
	powm[0]=1;
	for(i=1;i<=n;i++) powm[i]=powm[i-1]*m, powm[i]%=MOD;
	for(i=1;i<=n;i++)
	{
		for(j=i;j<=n;j+=i)
		{
			if(j != i && i >= 2) comp[j]=1;
			if(i >= 2 && !comp[i])
			{
				pi[j]/=i, pi[j]*=i-1;
			}
			d[j].push_back(i);
		}
	}
	ll ans=0;
	for(i=1;i<=n;i++)
	{
		ll neck=0;
		for(auto e : d[i])
		{
			neck+=pi[e]*powm[i/e]%MOD;
			neck%=MOD;
		}
		neck*=mul_inv(i*2, MOD);
		neck%=MOD;
		if(i%2 == 0)
		{
			neck+=mul_inv(4, MOD)*(m+1)%MOD*modpow(m, i/2)%MOD;
			neck%=MOD;
		}
		else
		{
			neck+=mul_inv(2, MOD)*modpow(m, (i+1)/2)%MOD;
			neck%=MOD;
		}
		ans+=neck, ans%=MOD;
	}
	ans++;
	ans%=MOD;
	cout<<ans;
	return 0;
}
//*/
# 결과 실행 시간 메모리 Grader output
1 Correct 4 ms 46408 KB Output is correct
2 Correct 4 ms 46408 KB Output is correct
3 Correct 9 ms 46408 KB Output is correct
4 Correct 0 ms 46408 KB Output is correct
5 Correct 0 ms 46408 KB Output is correct
6 Correct 6 ms 46408 KB Output is correct
7 Correct 4 ms 46408 KB Output is correct
8 Correct 3 ms 46408 KB Output is correct
9 Correct 4 ms 46408 KB Output is correct
10 Correct 3 ms 46408 KB Output is correct
11 Correct 6 ms 46408 KB Output is correct
12 Correct 9 ms 46408 KB Output is correct
13 Correct 9 ms 46408 KB Output is correct
14 Correct 4 ms 46408 KB Output is correct
15 Correct 4 ms 46408 KB Output is correct
16 Correct 0 ms 46408 KB Output is correct
17 Correct 4 ms 46408 KB Output is correct
18 Correct 4 ms 46408 KB Output is correct
19 Correct 0 ms 46408 KB Output is correct
20 Correct 6 ms 46408 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 6 ms 46408 KB Output is correct
2 Correct 6 ms 46408 KB Output is correct
3 Correct 3 ms 46408 KB Output is correct
4 Correct 76 ms 50632 KB Output is correct
5 Execution timed out 1000 ms 101452 KB Program timed out
6 Halted 0 ms 0 KB -