【数论】bzoj3601一个人的数论
来源:互联网 发布:基础网络的构成 编辑:程序博客网 时间:2024/05/27 14:15
题目链接
如果这题叫“一只蛤的神论”…
反正窝只会30分乱搞…
test1:暴力水过 test2:欧拉函数 test3:n/2*欧拉函数…
蒟蒻之前从未见过莫比乌斯反演与高斯消元如此奇妙的结合在一起…
所以来瞎bb一下…
明确目标:
利用莫比乌斯反演最基本的运用将
利用和式的性质将柿子变一下性。
此处
然后是窝觉得此题最令人惊叹的地方。
令
然后我们就可以用高斯消元来求这个多项式的系数…
当然也可以用拉格朗日插值法,理论复杂度比高斯消元优(然而窝写不来…逃…)
就是暴力求出
然后就知道辣
观察与
加之题目中的
根据莫比乌斯函数的定义,在这些柿子中只有
至此,问题解决完毕。
大概重新理一下解题过程。
1.用待定系数+高斯消元(拉格朗日插值法)求出某多项式的系数。
2.直接上计算式子
#include <iostream>#include <cstdio>#define mod 1000000007#define LL long long intusing namespace std;int n, d, p[1005][2];int sum[105], mat[155][155], A[155];int power(int a,LL pos){ int ans=1; if(pos<0)pos+=mod-1; for(;pos;pos>>=1, a=1ll*a*a%mod) if(pos&1)ans=1ll*ans*a%mod; return ans;}void init(){ for(int i=1;i<=d+2;++i) { sum[i]=(power(i,d)+sum[i-1])%mod; mat[i-1][d+2]=sum[i], mat[i-1][0]=1; for(int j=1;j<=d+1;++j)mat[i-1][j]=1ll*mat[i-1][j-1]*i%mod; } for(int i=0, j, k, tmp;i<=d+1;++i) { for(j=i;j<=d+1;++j)if(mat[j][i])break; if(i!=j)for(k=0;k<=d+2;++k)swap(mat[i][k],mat[j][k]); for(j=0;j<=d+1;++j) if(j!=i&&mat[j][i]) { tmp=1ll*mat[j][i]*power(mat[i][i],mod-2)%mod; for(k=0;k<=d+2;++k)mat[j][k]=(mat[j][k]-1ll*tmp*mat[i][k]%mod+mod)%mod; } } for(int i=0;i<=d+1;++i)A[i]=1ll*mat[i][d+2]*power(mat[i][i],mod-2)%mod;}int main(){ scanf("%d%d",&d,&n); init(); for(int i=1;i<=n;++i)scanf("%d%d",&p[i][0],&p[i][1]); int ans=0; for(int i=1, tmp;i<=d+1;++i) { tmp=1; for(int j=1;j<=n;++j) tmp=1ll*power(p[j][0],1ll*p[j][1]*i%(mod-1))*(1-power(p[j][0],d-i)+mod)%mod*tmp%mod; ans=(ans+1ll*A[i]*tmp%mod)%mod; } printf("%d\n",ans); return 0;}
1 0
- [BZOJ3601] 一个人的数论
- 【数论】bzoj3601一个人的数论
- 【BZOJ3601】一个人的数论,莫比乌斯反演+高斯消元
- bzoj3601 一个人的数论 莫比乌斯反演+高斯消元
- 数论
- 数论
- 数论
- 数论
- 数论
- 数论
- 数论
- 数论
- 数论
- 数论
- 数论
- 数论 ?
- 数论
- 数论
- 第四课讲C语言程序设计之体验---项目一《计算长方形的周长和面积》
- ubuntu卡在登陆界面,无法登入桌面。
- JDBC笔记-Meta类
- jsp内置对象——request
- 呼叫流程
- 【数论】bzoj3601一个人的数论
- byte order form webrtc soure file
- POJ 1061 一次同余方程
- AFNetworking 3.0 使用详解和源码剖析
- docker ubuntu apt-get update 出错Hash Sum mismatch
- ubuntu16.04硬盘重装
- Mysql使用(1)
- java带界面数据库连接简易Atm项目
- 一个80后程序员的转行创业日记14--两个月一塌糊涂