E: Number
题意:
给一个整数n(5<=n<=59),求1/n mod x,其中x为将n分解若干个正整数然后使得乘积最大的那个积。
题解:
求n逆元!
首先要将x给算出来,
有两种方法:
1.直接套公式:
当n=3k,k>1的时候,最大的积是3^k
当n=3k+1,k>1的时候,最大的积是3^(k-1)4
当n=3k+2,k>1的时候,最大的积是3^k2
2.递推
首先可以得出1、2、3的x分别为1、2、3,细想一下之后的数一定是之前已经有过的数乘出来,并且为两个数相乘:
x(4)=x(2)x(2)=22=4
x(5)=x(2)x(3)=23=6
x(6)=x(3)x(3)=33=9
x(7)=x(3)x(4)=34=12
x(8)=x(6)x(2)=92=18
…
就可推出所有的x(n)了。
然就行是求逆元了,用一个扩展欧欧几里得求出n的逆元即可,然后ans=*inv[n] mod x即可。
#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll n;
ll x,y;
ll poww(ll a,ll b)
{
ll ans=1;
while(b)
{
if(b&1) ans=ans*a;
a=a*a;
b>>=1;
}
return ans;
}
void ex_gcd(ll a,ll b)
{
if(b==0)
{
x=1,y=0;
return ;
}
ex_gcd(b,a%b);
ll X=x,Y=y;
x=Y;
y=X-a/b*Y;
}
int main()
{
scanf("%lld",&n);
ll xx;
/*
当n=3k,k>1的时候,最大的积是3^k
当n=3k+1,k>1的时候,最大的积是3^(k-1)*4
当n=3k+2,k>1的时候,最大的积是3^k*2
*/
if(n%3==0) xx=poww(3,n/3);
else if((n-1)%3==0) xx=poww(3,(n-1)/3-1)*4;
else if((n-2)%3==0) xx=poww(3,(n-2)/3)*2;
ex_gcd(n,xx);
x=(x%xx+xx)%xx;
printf("%lld",x);
return 0;
}
#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll dp[65];
ll n;
ll x,y;
void ex_gcd(ll a,ll b)
{
if(b==0)
{
x=1,y=0;
return ;
}
ex_gcd(b,a%b);
ll X=x,Y=y;
x=Y;
y=X-a/b*Y;
}
int main()
{
scanf("%lld",&n);
dp[1]=1,dp[2]=2,dp[3]=3;
for(int i=4;i<=n;i++)
{
dp[i]=0;
for(int j=1;j<=i/2+1;j++)
dp[i]=max(dp[i],dp[j]*dp[i-j]);
}
ex_gcd(n,dp[n]);
if(x<0)
x+=(-x/dp[n]+1)*dp[n];
else if(x>0)
x%=dp[n];
//x=(x%dp[n]+dp[n])%dp[n];
printf("%lld\n",x);
return 0;
}