# |
Submission time |
Handle |
Problem |
Language |
Result |
Execution time |
Memory |
19816 |
2016-02-25T05:54:34 Z |
tonyjjw |
팔찌 (kriii4_V) |
C++14 |
|
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;
}
//*/
# |
Verdict |
Execution time |
Memory |
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 |
# |
Verdict |
Execution time |
Memory |
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 |
- |