两道数论题推得我快吐了
1001
这题说实话第一眼看到,马上想到去网上找直角四面体的公式,然后找到了:
1/h^ 2 =1/a ^ 2+1/b ^ 2+1/c ^ 2;
E(1/h^ 2)=cnt*(1/1 ^ 2+ 1/2 ^ 2+……1/n ^ 2)/总可能性(n^3);
求(1,n)中每个数字出现了多少次
一共a,b,c三个位置每个位置有n种选择,总可能性为n^3.
总可能性有n^3,出现的数的总数是3n ^ 3,每个数出现的的个数就是3n ^ 2;
最后得到公式E(1/h^2)=3(1/1 ^ 2+ 1/2 ^ 2+……1/n ^2)/n;
记得求逆元时先求逆元再平方,不要求平方之后在逆元
#include<bits/stdc++.h>
#define ll unsigned long long
using namespace std;
const int N=6e6+10;
const int inf=0x3f3f3f3f;
const ll mod=998244353;
inline int rd()
{
register int s=0,w=1;
register char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-') w=-1;
ch=getchar() ;
}
while(ch<='0'&&ch<='9')
{
s=s*10+ch-'0';
ch=getchar();
}
return s*w;
}
ll f[N],re[N];
ll fpow(ll a,ll b,ll m)
{
ll ans=1,base=a%m;
while(b)
{
if(b&1) ans=(ans*base)%m;
base=(base*base)%m;
b>>=1;
}
return ans;
}
void ff()
{
for(int i=1;i<=N;i++)
{
//ll u=(ll)(i*i)%mod;
f[i]=fpow(i,mod-2,mod)%mod;
f[i]=(f[i]*f[i])%mod;
re[i]=(re[i-1]+f[i])%mod;
}
}
int main()
{
int t;
scanf("%d",&t);
ff();
while(t--)
{
ll n;
scanf("%llu",&n);
ll x1=fpow(n,mod-2,mod);
ll ans=((re[n]*3)%mod*x1)%mod;
printf("%llu\n",ans);
}
return 0;
}
1009
说是4种折法,其实就2种,上下和左右,上下折k次会出现2^k条横向剪痕,左右折n-k次会出现
2^ (n-k)条剪痕,于是会剪出(2 ^ k+1)(2 ^ (n-k)+1)个纸片。之后求期望,我做了半天都没做出
该死的数学
大佬的数学推导过程:https://blog.csdn.net/oampamp1/article/details/107796100
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+10;
const int mod=998244353;
const int inf=1e9+7;
ll quick(ll a,ll b){
ll ans=1;
while(b){
if(b&1)ans=ans*a%mod;
a=a*a%mod;
b>>=1;
}
return ans;
}
int main () {
int t;
scanf("%d",&t);
while(t--){
ll n;
scanf("%lld",&n);
ll ans=(1ll+quick(2,n)+(2*quick(3,n)%mod)%mod*quick(quick(2,mod-2),n)%mod)%mod;
printf("%lld\n",ans);
}
}