博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
hdu4059The Boss on Mars 容斥原理
阅读量:6072 次
发布时间:2019-06-20

本文共 1193 字,大约阅读时间需要 3 分钟。

//求1到n之间与n互质的数的四次方的和
//segma(n^4) = (6n^5+15n^4+10n^3-n)/30
//对于(a/b)%mod能够转化为(a*inv(b))%mod
//inv(b)为b的逆元
//由费马小定理a^(p-1) = 1(modp) a , p 互质可得
//30的逆元为30^(mod-2)
//由容斥原理非常easy得到与n不互质的数之和为
//对于全部的n的素数因子
//一个素数因子的全部数的四次方之和-有两个素数因子的全部数的四次方之和+有三个。

。。。

//然后用总的减去不互质的即为互质的
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std ;
const int maxn = 1010 ;
typedef __int64 ll ;
const __int64 mod = 1e9+7 ;
int p[maxn] ;
int len ;
ll  inv ;
void get_prime(ll n)
{
    len = 0 ;
    for(int i = 2;i*i <= n;i++)
    {
        if(n%i == 0)p[++len] = i;
        while(n%i==0)n/=i ;
    }
    if(n>1)p[++len] = n ;
}
ll Pow(ll a , ll b)
{
    ll c = 1 ;
    while(b)
    {
        if(b&1)c = (c*a)%mod;
        a = (a*a)%mod;
        b >>= 1;
    }
    return c ;
}
ll dfs(int pos , int n)
{
    ll ans = 0 ;
    for(int i = pos;i <= len ;i++)
    {
        ll t = n/p[i] ;
        ll t_1 = ((((6*Pow(t,5) + 15*Pow(t,4) + 10*Pow(t,3) - t))%mod)*inv)%mod;
        ans = (ans + (Pow(p[i] , 4)*(t_1- dfs(i+1 , n/p[i]))))%mod;
    }
    return ans ;
}
int main()
{
    int T ;
    scanf("%d" , &T) ;
    inv = Pow(30 , mod - 2) ;
    while(T--)
    {
        ll n ;
        scanf("%I64d" , &n) ;
        get_prime(n) ;
        ll ans = (((((6*Pow(n,5) + 15*Pow(n,4) + 10*Pow(n,3) - n))%mod)*inv)%mod - dfs(1 , n))%mod ;
        printf("%I64d\n" , (ans+mod)%mod);
    }
    return  0 ;
}

转载地址:http://lgngx.baihongyu.com/

你可能感兴趣的文章
yum mysql on centos 7
查看>>
Jquery动态绑定事件处理函数 bind / on / delegate
查看>>
node 下查看安装插件的最新版本号的方法
查看>>
雨课堂知识点总结(十四)
查看>>
16个时髦的扁平化设计+视觉滚差的网站模板
查看>>
tesseract text2image windows
查看>>
Linux系统下如何配置JDK环境变量
查看>>
CodeForces 767C Garland
查看>>
postman/jmeter 请求https接口配置问题
查看>>
事务的学习,从jdbc开始:jdbc对事务的支持与实现
查看>>
[人工智能]追逐与躲闪
查看>>
CentOS上搭建Yii2 --2017
查看>>
防范DDOS攻击脚本
查看>>
在阿里云上挂在/data脚本
查看>>
解决文件大小上传限制
查看>>
Anconda 3.7安装以及使用详细教程
查看>>
three.js之让物体动起来方式(二)移动物体
查看>>
ARM总线架构
查看>>
Nexus构建
查看>>
php读入mysql数据并以表格形式显示(表单实现无刷新提交)
查看>>