题目大意:求n!在k进制下末尾0的个数。
做法:我们对k进行质因数分解,然后用k的质因数去除 n!,然后取num2/num1的最小值就可以了。
看数据范围,对k的质因数分解是要用rho大整数分解来做的,然而数据好像比较弱,筛到1e7就可以过的QAQ。
对于rho的应用见我前边一篇博客 欧拉函数专项
#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#define LL long long
using namespace std;
LL n,k,ans,cnt,m,MOD;
bool np[10000009];
LL prime[10000009],num[10000009],q[10000009],r;
void su()
{
np[1]=1,np[0]=1;
for(int i=2;i<=10000001;i++)
{
if(!np[i]) prime[++cnt]=i;
for(int j=1;j<=cnt&&1ll*(prime[j]*i)<=10000001;j++)
{
np[i*prime[j]]=1;
if(i%prime[j]==0) break;
}
}
}
void pre()
{
LL p=k;
for(int i=1;i<=cnt;i++)
{
LL x=prime[i],f=0;
while(p%x==0&&p)
{
p/=x;
num[x]++;
f=1;
}
if(f) q[++r]=prime[i];
if(!p) break;
}
}
LL count(LL x)
{
LL s=0,p=x;
while(p<=n)
{
s+=n/p;
p*=x;
}
return s;
}
int main()
{
scanf("%lld%lld",&n,&k);
su();
pre();
ans=1e18;
for(int i=1;i<=r;i++)
{
ans=min(ans,count(q[i])/num[q[i]]);
}
if(ans==1e18) ans=0;
printf("%lld",ans);
return 0;
}
清晰版题目描述
小强拿到一个3×n的数组,要在每一列选一个数(或者不选),满足以下条件:
1.如果在第一行选,那它必须大于等于上一个数
2.如果在第二行选,那么必须小于等于上一个数
3.如果在第三行选,对于连续的一段在第三行选的数,必须满足方向相同(都小于等于上一个数或者都大于等于上一个数)
20分:暴力枚举每一个位置 不选,选1,选2,选3
60分:用dp[i][j]表示第 i 列选第 j 行达到的最长长度。
设 dp i,j 表示:第 i 个位置选择第 j 个序列的数,前 i 个位置
所能构成的最长子序列的长度。特别的,dp i,2 表示第 i 个位
置选择第三个序列的数,且大于等于前一个数;dp i,3 表示
第 i 个位置选择第三个序列的数,且小于等于前一个数。
dp i,1 = max(dp j,k ),(1≤ j < i,1≤ k ≤ 4,a j,k ≤ a i,1 ) + 1
dp i,2 = max(dp j,k ),(1 ≤ j < i,1 ≤ k ≤ 4,a j,k ≥ a i,2 ) + 1
dp i,3 = max(dp j,k ),(1≤ j < i,1 ≤ k ≤ 4,k≠4,a j,k ≤ a i,3 ) + 1
dp i,4 = max(dp j,k ),(1≤ j < i,1≤ k ≤ 4,k≠3,a j,k ≥ a i,4) + 1
边界条件:dp 0,j = 0。
状态数 O(n),转移复杂度 O(n),总复杂度 O(n^2 )。
100分:在60分的基础上,在选取最大值的时候,用树状数组来优化 (需要离散化)。
在这里要用8个树状数组,为了节省代码量,我们只用树状数组来计算前面的最大值 ( 所以需要加上一些特殊处理 ) 例如当需要a j,k>=a i,2时我们可以按照m-a i,2来查询(a j,k>=a i,2 即 m-a j,k <= m-a i,2)。
60分裸的dp
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define M 100009
#define LL long long
using namespace std;
int n,a[M][5],dp[M][5],ans;
int main()
{
scanf("%d",&n);
for(int j=1;j<=3;j++)for(int i=1;i<=n;i++)
scanf("%d",&a[i][j]);
for(int i=1;i<=n;i++) a[i][4]=a[i][3];
for(int i=1;i<=n;i++)
{
for(int j=1;j<i;j++)
{
for(int k=1;k<=4;k++)
{
if(a[j][k]<=a[i][1]) dp[i][1]=max(dp[i][1],dp[j][k]+1);
if(a[j][k]>=a[i][2]) dp[i][2]=max(dp[i][2],dp[j][k]+1);
if(a[j][k]>=a[i][3]&&k!=4) dp[i][3]=max(dp[i][3],dp[j][k]+1);
if(a[j][k]<=a[i][4]&&k!=3) dp[i][4]=max(dp[i][4],dp[j][k]+1);
}
}
}
for(int i=1;i<=4;i++) ans=max(ans,dp[n][i]);
printf("%d",ans);
return 0;
}
100分
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define M 100009
#define LL long long
using namespace std;
int n,a[M][5],dp[M][5],ans;
int te[10][M*3],t[M*3],tot,m;
inline int read()
{
register int x=0,f=1; char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}
return x*f;
}
int query(int id,int x)
{
int s=0;
for(int i=x;i;i-=(-i)&i)
s=max(te[id][i],s);
return s;
}
void update(int id,int x,int p)
{
for(int i=x;i<=m;i+=i&(-i))
te[id][i]=max(te[id][i],p);
}
int main()
{
scanf("%d",&n);
for(int j=1;j<=3;j++)
for(int i=1;i<=n;i++)
a[i][j]=read(),t[++tot]=a[i][j];
sort(t+1,t+tot+1);
m=unique(t+1,t+tot+1)-t-1;
for(int j=1;j<=3;j++)
for(int i=1;i<=n;i++)
a[i][j]=lower_bound(t+1,t+m+1,a[i][j])-t;//离散化
for(int i=1;i<=n;i++)
{
dp[i][1]=max(dp[i][1],query(1,a[i][1])+1);
dp[i][1]=max(dp[i][1],query(3,a[i][1])+1);
dp[i][1]=max(dp[i][1],query(5,a[i][1])+1);
dp[i][1]=max(dp[i][1],query(7,a[i][1])+1);
dp[i][2]=max(dp[i][2],query(2,m-a[i][2]+1)+1);
dp[i][2]=max(dp[i][2],query(4,m-a[i][2]+1)+1);
dp[i][2]=max(dp[i][2],query(6,m-a[i][2]+1)+1);
dp[i][2]=max(dp[i][2],query(8,m-a[i][2]+1)+1);
dp[i][3]=max(dp[i][3],query(1,a[i][3])+1);
dp[i][3]=max(dp[i][3],query(3,a[i][3])+1);
dp[i][3]=max(dp[i][3],query(5,a[i][3])+1);
dp[i][4]=max(dp[i][4],query(2,m-a[i][3]+1)+1);
dp[i][4]=max(dp[i][4],query(4,m-a[i][3]+1)+1);
dp[i][4]=max(dp[i][4],query(8,m-a[i][3]+1)+1);
update(1,a[i][1],dp[i][1]);
update(3,a[i][2],dp[i][2]);
update(5,a[i][3],dp[i][3]);
update(7,a[i][3],dp[i][4]);
update(2,m-a[i][1]+1,dp[i][1]);
update(4,m-a[i][2]+1,dp[i][2]);
update(6,m-a[i][3]+1,dp[i][3]);
update(8,m-a[i][3]+1,dp[i][4]);
}
for(int i=1;i<=4;i++) ans=max(ans,dp[n][i]);
printf("%d",ans);
return 0;
}