生活随笔
收集整理的這篇文章主要介紹了
筛法 V - 一道水题
小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
一天，szb 在上學(xué)的路上遇到了灰太狼。
灰太狼：幫我們做出這道題就放了你。
szb：什么題？
灰太狼：求一個(gè)能被 [1,n] 內(nèi)所有數(shù)整除的最小數(shù)字，并對(duì) 100000007 取模。
szb：這題太水了，就讓我小弟來(lái)做好了。
然后你就光榮的接受了這個(gè)任務(wù)。
Input
一行一個(gè)數(shù) n。
Output
一行一個(gè)數(shù) ans。
Example
樣例輸入
10
樣例輸出
2520
Hint
n≤108
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<algorithm>
#include<queue>
#include<map>
#include<cstdio>
#include<vector>
#include<set>
#include<cstring>
#include<cstdlib>
#include<time.h>
#include<stack>
#define INF 0x3f3f3f3f
using namespace std
;
const int maxn
=100;
typedef long long ll
;
ll mod
=1e8+7;
ll num
,f
[maxn
];
struct mat
{ll m
[maxn
][maxn
];
}ans
,a
;
ll
mode(ll a
,ll b
){ll sum
=1;a
=a
%mod
;while(b
){if(b
&1)sum
=(sum
*a
)%mod
;b
/=2;a
=(a
*a
)%mod
;}return sum
;
}
ll
inv(ll a
,ll b
){return(a
*mode(b
,mod
-2))%mod
;
}
ll
gcd(ll a
, ll b
) {return b
?gcd(b
, a
%b
):a
;
}
ll
lcm(ll x
, ll y
){return x
/ gcd(x
, y
)*y
;
}
ll
exgcd(ll a
,ll b
,ll
&x
,ll
&y
){if(b
==0){x
=1;y
=0;return a
;}int r
=exgcd(b
,a
%b
,x
,y
);int temp
=y
;y
=x
-(a
/b
)*y
;x
=temp
;return r
;
}
mat
mul(mat a
,mat b
,int n
){mat c
;memset(c
.m
,0,sizeof(c
.m
));for(int i
=1;i
<=n
;i
++)for(int j
=1;j
<=n
;j
++)for(int k
=1;k
<=n
;k
++)c
.m
[i
][j
]=(c
.m
[i
][j
]+(a
.m
[i
][k
]*b
.m
[k
][j
])%mod
)%mod
;return c
;
}
mat
_power(mat a
,int b
,int n
){memset(ans
.m
,0,sizeof(ans
.m
));for(int i
=1;i
<=n
;i
++)ans
.m
[i
][i
]=1;while(b
){if(b
&1)ans
=mul(a
,ans
,n
);a
=mul(a
,a
,n
);b
>>=1;}return ans
;
}void init(int n
){int vis
[maxn
];for(int i
=2;i
<=n
;i
++){if(vis
[i
])continue;for(int j
=i
*2;j
<=n
;j
+=i
)vis
[j
]=1;}
}
int oula(int x
)
{int i
,j
;int num
= x
;for(i
= 2; i
*i
<= x
; i
++){if(x
% i
== 0){num
= (num
/i
)*(i
-1);while(x
% i
== 0){x
= x
/ i
;}}}if(x
!= 1) num
= (num
/x
)*(x
-1);return num
;
}
const int mx
=1e8+5;
bool vis
[mx
];
int z
[1000000];
int main()
{int n
,cnt
=0;scanf("%d",&n
);ll ans
=1;memset(vis
,0,sizeof(vis
));for(int i
=2;i
<=n
;i
++){if(!vis
[i
]){z
[cnt
++]=i
;for(ll j
=i
;j
<=n
;j
*=i
)ans
=ans
*i
%mod
;}for(int j
=0;j
<cnt
;j
++){ll v
=i
*z
[j
];if(v
>n
)break;vis
[v
]=1;if(i
%z
[j
]==0)break;}}printf("%lld\n",ans
);return 0;
}
總結(jié)
以上是生活随笔為你收集整理的筛法 V - 一道水题的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。
如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。
