poj3910欧拉函数的矩阵应用

最新推荐文章于 2021-08-04 15:15:41 发布

&小鹏鹏

最新推荐文章于 2021-08-04 15:15:41 发布

阅读量601

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/u012785169/article/details/44839049

本文介绍了一种基于素数筛的算法实现，该算法用于处理包含大量整数运算的问题，并通过预先筛选出一定范围内的所有素数来加速计算过程。文章详细展示了如何使用 C 语言实现这一算法，包括如何有效地标记合数、处理输入数据以及计算特定数学表达式的值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

//有点问题
#include <stdio.h> 
#include <math.h>
const int mo = 1000000007;
long long i, n, j, a[101000], bj[101000], k, tot, ans;
int main()
{
    //freopen("worry.in", "r", stdin);
    //freopen("worry.out", "w", stdout);
    while(1)
    {
        scanf("%I64d", &n);
    for (bj[a[tot = 1] = i = 2] = 1; i <= 100000; ++i)
    {
if (!bj[i]) bj[i] = 1, a[++tot] = i;
for (j = 1; j <= tot && a[j] * i <= 100000; ++j)
{
   bj[a[j] * i] = 1;
   if (i % a[j] == 0) break;
}
    }
    for (i = ans = 1; i <= n; ++i)
    {
scanf("%I64d", &k); int h = k;
for (j = 1; j <= tot && a[j] * a[j] <= h; ++j)
   for (int ok = 1; k % a[j] == 0;)
   {
if (ok) ans = ans * (a[j] - 1) % mo;
else ans = ans * a[j];
k /= a[j];
ok = 0;
   }
if (k > 1) ans = ans * (k - 1) % mo;
    }
    printf("%I64d\n", ans);
    }
    return 0;
}