H题
题解:
挂一下原题的链接原题
网上题解很详细了。代码和网上题解思路一样。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 2e5 + 10;
//int nxt[maxn][10],fa[maxn],len[maxn],cnt[maxn],pos[maxn];
//int last = 1,tot = 1;
//void insert(int x){
// int p = last, np = ++tot;
// last = np; len[np] = len[p] + 1;
// for(;p && !nxt[p][x];p=fa[p]) nxt[p][x] = np;
// if(!p) fa[np] = 1;
// else{
// int q = nxt[p][x];
// if(len[p] + 1 == len[q]) fa[np] = q;
// else{
// int nq = ++tot; len[nq] = len[p] + 1;
// memcpy(nxt[nq],nxt[q],sizeof(nxt[q]));
// fa[nq] = fa[q]; fa[q] = fa[np] = nq;
// for(;nxt[p][x] == q;p=fa[p]) nxt[p][x] = nq;
// }
// }
//}
int vis[maxn],prime[maxn],cnt;
void init()
{
for(int i=2;i<maxn;i++){
if(!vis[i]) prime[++cnt]=i;
for(int j=1;j<=cnt&&i*prime[j]<maxn;j++){
vis[i*prime[j]]=1;
if(i%prime[j]==0) break;
}
}
}
//char str[maxn];
typedef pair<int,int> pp;
int used[maxn];
signed main(){
init();
int t; scanf("%d",&t);
while(t--){
vector<int> g;
vector<pp> res;
int n; scanf("%d",&n);
for(int i=2;i<=n;i++) used[i]=0;
for(int i=1;i<=cnt;i++){
if(prime[i]*2<=n) g.push_back(prime[i]);
else break;
}
reverse(g.begin(),g.end());
for(auto it:g){
queue<int> tmp;
tmp.push(it);
used[it]=1;
for(int i=it*3;i<=n;i+=it){
if(!used[i]) {
tmp.push(i);
used[i]=1;
}
}
if((int)tmp.size()%2==1&&it*2<=n){
if(!used[it*2]){
tmp.push(it*2);
used[it*2]=1;
}
}
while((int) tmp.size()){
int a=tmp.front(); tmp.pop();
if(tmp.size()==0) break;
int b=tmp.front(); tmp.pop();
res.push_back({a,b});
}
}
queue<int> tmp;
for(int i=2;i<=n;i+=2){
if(!used[i]){
tmp.push(i);
used[i]=1;
}
}
while((int)tmp.size()){
int a=tmp.front(); tmp.pop();
if(tmp.size()==0) break;
int b=tmp.front(); tmp.pop();
res.push_back({a,b});
}
printf("%d\n",res.size());
for(auto it:res) printf("%d %d\n",it.first,it.second);
}
}
C题:
题解:
我们从后往前插入,如果遇见前面位置
i
i
i的字母比后面位置
j
j
j的字母大,所以我们可以在当前
i
i
i位置插入自动机last-i个这个字母。最后统计子串数目。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 2e6+10;
char s[maxn];
struct SuffixAutoMation
{
int last=1,cnt=1;//cnt表示状态
int trans[maxn<<1][26],slink[maxn<<1],l[maxn<<1];
int t[maxn<<1],a[maxn<<1];
int num[maxn<<1];
inline void add(int x)
{
int p = last,np = ++cnt;last = np;l[np] = l[p]+1;
for(;p && !trans[p][x];p=slink[p]) trans[p][x] = np;
if(!p) slink[np] = 1;
else{
int q = trans[p][x];
if(l[p]+1 == l[q]) slink[np] = q;
else
{
int nq = ++cnt;l[nq] = l[p]+1;
memcpy(trans[nq],trans[q],sizeof(trans[q]));
slink[nq] = slink[q];
slink[q] = slink[np] = nq;
for(;trans[p][x] == q; p = slink[p]) trans[p][x] = nq;
}
}
num[np] = 1;
}
void build()
{
scanf("%s",s+1);
int len = strlen(s+1);
last = cnt = 1;
for(int i=1;i<=len;i++) add(s[i]-'a');
}
void topsort()
{
int ans = 0;
for(int i=1;i<=cnt;i++) t[l[i]]++;
for(int i=1;i<=cnt;i++) t[i]+=t[i-1];
for(int i=1;i<=cnt;i++) a[t[l[i]]--] = i;
for(int i=cnt;i>=1;i--){
int x = a[i];
num[slink[x]]+=num[x];
if(num[x]>1) ans = max(ans,l[x]*num[x]);
}
printf("%d\n",ans);
}
}sam;
int pos[maxn];
signed main()
{
scanf("%s",s+1);
int n=strlen(s+1);
pos[n+1]=1;
stack<int> st;
for(int i=n;i>=1;i--){
while (!st.empty()&&s[st.top()]<s[i]) st.pop();
int k=st.empty()?n+1:st.top();
sam.last=pos[k];
for(int j=i;j<k;j++) sam.add(s[i]-'a');
pos[i]=sam.last;
st.push(i);
}
ll ans=0;
for(int i=2;i<=sam.cnt;i++) ans+=sam.l[i]-sam.l[sam.slink[i]];
printf("%lld\n",ans);
}