GTW likes function
Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others)
Total Submission(s): 243 Accepted Submission(s): 134
Problem Description
Now you are given two definitions as follows.
f(x)=∑xk=0(−1)k22x−2kCk2x−k+1,f0(x)=f(x),fn(x)=f(fn−1(x))(n≥1)
Note that φ(n) means Euler’s totient function.( φ(n) is an arithmetic function that counts the positive integers less than or equal to n that are relatively prime to n.)
For each test case, GTW has two positive integers — n and x , and he wants to know the value of the function φ(fn(x)) .
f(x)=∑xk=0(−1)k22x−2kCk2x−k+1,f0(x)=f(x),fn(x)=f(fn−1(x))(n≥1)
Note that φ(n) means Euler’s totient function.( φ(n) is an arithmetic function that counts the positive integers less than or equal to n that are relatively prime to n.)
For each test case, GTW has two positive integers — n and x , and he wants to know the value of the function φ(fn(x)) .
Input
There is more than one case in the input file. The number of test cases is no more than 100. Process to the end of the file.
Each line of the input file indicates a test case, containing two integers, n and x , whose meanings are given above. (1≤n,x≤1012)
Each line of the input file indicates a test case, containing two integers, n and x , whose meanings are given above. (1≤n,x≤1012)
Output
In each line of the output file, there should be exactly one number, indicating the value of the function
φ(fn(x))
of the test case respectively.
Sample Input
1 1 2 1 3 2
Sample Output
2 2 2
问题描述
现在给出下列两个定义: f(x)=f0(x)=∑k=0x(−1)k22x−2kC2x−k+1k,fn(x)=f(fn−1(x))(n≥1) φ(n)为欧拉函数。指的是不超过n的与n互质的正整数个数。 对于每组数据,GTW有两个正整数n,x,现在他想知道函数φ(fn(x))的值。
输入描述
输入有多组数据,不超过100组。每数据输入一行包含2个整数组n和x。(1≤n,x≤1012)
输出描述
对于每组数据输出一行,表示函数φ(fn(x))的值。
打表代码:
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <vector>
#define INF 0x3f3f3f3f
#define eps 1e-8
#define MAXN (50000+1)
#define MAXM (5000000000)
#define Ri(a) scanf("%d", &a)
#define Rl(a) scanf("%lld", &a)
#define Rf(a) scanf("%lf", &a)
#define Rs(a) scanf("%s", a)
#define Pi(a) printf("%d\n", (a))
#define Pf(a) printf("%.2lf\n", (a))
#define Pl(a) printf("%lld\n", (a))
#define Ps(a) printf("%s\n", (a))
#define W(a) while(a--)
#define CLR(a, b) memset(a, (b), sizeof(a))
#define MOD 1000000007
#define LL long long
#define lson o<<1, l, mid
#define rson o<<1|1, mid+1, r
#define ll o<<1
#define rr o<<1|1
using namespace std;
LL C[110][110];
void getC()
{
for(int i = 0; i <= 50; i++)
C[i][0] = 1, C[i][i] = 1;
for(int i = 1; i <= 50; i++)
for(int j = 1; j < i; j++)
C[i][j] = C[i-1][j-1] + C[i-1][j];
}
LL Pow(LL a, LL n)
{
LL ans = 1;
while(n)
{
ans *= a;
n--;
}
return ans;
}
LL f(LL x)
{
LL ans = 0;
for(LL k = 0; k <= x; k++)
{
if(k & 1LL)
ans -= Pow(2, 2*x-2*k) * C[2*x-k+1][k];
else
ans += Pow(2, 2*x-2*k) * C[2*x-k+1][k];
}
return ans;
}
int main()
{
getC();
for(LL i = 1; i <= 20LL; i++)
Pl(f(i));
return 0;
}
规律明显为f(x) = x + 1,这样f n(x) = x + n + 1。下面就是求解欧拉函数了。
AC代码:数据很大,先除后乘。
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <vector>
#define INF 0x3f3f3f3f
#define eps 1e-8
#define MAXN (50000+1)
#define MAXM (5000000000)
#define Ri(a) scanf("%d", &a)
#define Rl(a) scanf("%lld", &a)
#define Rf(a) scanf("%lf", &a)
#define Rs(a) scanf("%s", a)
#define Pi(a) printf("%d\n", (a))
#define Pf(a) printf("%.2lf\n", (a))
#define Pl(a) printf("%lld\n", (a))
#define Ps(a) printf("%s\n", (a))
#define W(a) while(a--)
#define CLR(a, b) memset(a, (b), sizeof(a))
#define MOD 1000000007
#define LL long long
#define lson o<<1, l, mid
#define rson o<<1|1, mid+1, r
#define ll o<<1
#define rr o<<1|1
using namespace std;
LL euler(LL n)
{
LL ans = n;
for(LL i = 2; i * i <= n; i++)
{
if(n % i == 0)
{
ans = ans / i * (i-1);
while(n % i == 0)
n /= i;
}
}
if(n > 1) ans = ans / n * (n-1);
return ans;
}
int main()
{
LL n, x;
while(scanf("%lld%lld", &n, &x) != EOF){
printf("%I64d\n", euler(n+x+1));
}
return 0;
}