lcp计数 (Trie)

lcp

思路:
看着挺唬人的,其实是道水题。
考虑每个子串作为lcp的贡献,就是包含这个子串的字符串的两两组合。

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#include <iostream>
#define N 1000005
#define LL long long
using namespace std;

int n, rot, cnt; 
int c[N][26], sum[N];
char s[N];
LL ans;

void insert(char *s){
    int p = rot, len = strlen(s);
    for(int i=0; i<len; i++){
        int x = c[p][s[i] - 'a'];
        if( !x ) c[p][s[i] - 'a'] = ++cnt;
        p = c[p][s[i] - 'a'];
        sum[p]++;
    }
}

void query(){
    int p = rot;
    queue<int> q;
    q.push( p );
    while( !q.empty() ){
        int u = q.front(); q.pop();
        for(int i=0; i<26; i++){
            int v = c[u][i];
            if( !v ) continue;
            ans += (1LL * sum[v] * (sum[v] - 1) / 2);
            //一共有sum[v]个字符串中,包含v这个子串
            //那么v作为lcp的贡献就有C(2,sum[v]) 
            q.push( v );
        }
    }
}

int main(){
    scanf("%d", &n);
    rot = cnt = 1;
    for(int i=1; i<=n; i++){
        scanf("%s", s);
        insert( s );
    }
    query();
    cout << ans << endl;
}
# P9218 「TAOI-1」Apollo ## 题目背景 Execution. 这题原来叫 `std::cout << std::fixed << std::setprecision(1) << 6.5 << '\n';`。 [被当事人删掉的图片.jpg] **【Upd 2023/04/15 21:42】添加了一组 Hack 数据位于 Subtask 2,#13。现在所有赛时的 $\bm{50}$ 分提交理论上均只能获得 $\bm{30}$ 分。** ## 题目描述 给出 $n$ 个十进制小数 $a_1 \dots a_n$。 对于一个 **数** $a$,定义其精度 $f(a)$ 表示最小的非负整数 $k$ 使得 $10^k\times a$ 为整数;对于整数 $a$,定义 $f(a) = 0$。对于两个十进制小数 $a, b$,定义 $g(a, b)$ 为对于所有 **数** $c \in [\min(a,b), \max(a,b)]$ 的 $f(c)$ 的最小值。 对于所有 $1 \leq i \leq n$,你需要求出 $\sum\limits_{j=1}^ng(a_i, a_j)$ 的值并输出。 定义十进制小数是一个含有整数部分和小数部分的数,其小数部分不为 $0$。之所以描述得这么愚蠢是因为保证输入的每个数都有小数点,但是好像无论什么写法都会有人提出不满,真是抱歉了。~~所以还是得看看被当事人删掉的图片。所以我在这里写闲话有人看得到吗。~~ ## 输入格式 第一行一个整数 $n$。 接下来 $n$ 行,每行一个十进制小数 $a_i$。 ## 输出格式 $n$ 行,每行一个整数,分别表示 $i = 1 \dots n$ 对应的答案。 ## 输入输出样例 #1 ### 输入 #1 ``` 5 11.4514 11.4778 11.1338 11.1236 11.4789 ``` ### 输出 #1 ``` 10 11 9 9 11 ``` ## 输入输出样例 #2 ### 输入 #2 ``` 8 1.1145 1.114 1.1145 1.514 1.19198 1.1154 1.114 1.1145 ``` ### 输出 #2 ``` 24 21 24 10 18 22 21 24 ``` ## 说明/提示 ### 数据范围 **本题采用捆绑测试。** 令 $\sum\limits_{i=1}^n f(a_i) = t$。 - Subtask 1(15 points):$a_i \leq 10$,$n \leq 5$,$t \leq 10$。 - Subtask 2(15 points):$a_i \leq 10$,$n \leq 100$,$t \leq 1000$。 - Subtask 3(20 points):$n \leq 1722$。 - Subtask 4(15 points):$a_i \leq 1$。 - Subtask 5(35 points):无特殊限制。 对于所有数据,$0 \lt a_i \lt 10^9$,$1 \leq n \leq 10^5$,$1 \leq t \leq 3 \times 10^6$,**保证 $\bm{a_i}$ 没有后导 [$\color{black}\bm 0$](https://cdn.luogu.com.cn/upload/image_hosting/pg66fm1z.png),不保证 $\bm{a_i}$ 互不相同**。 ### 样例解释 以 $i = 1$ 为例: + $j = 1$:取 $c = 11.4514$,$f(c) = 4$; + $j = 2$:取 $c = 11.46$,$f(c) = 2$; + $j = 3$:取 $c = 11.2$,$f(c) = 1$; + $j = 4$:取 $c = 11.3$,$f(c) = 1$; + $j = 5$:取 $c = 11.47$,$f(c) = 2$。 故 $\sum\limits_{j=1}^n g(a_1, a_j) = 4 + 2 + 1 + 1 + 2 = 10$。对于同一个 $j$,上文给出的所有 $c$,都可能有其它的不同的 $c$ 满足 $f(c)$ 同样最小。 #include <bits/stdc++.h> #define ll long long using namespace std; const ll N=1e5,M=3e6; struct node{ string x,y; }a[N+5]; ll n; ll tr[3*M+5][10],vis[2*M+5],root[N+5],ed[2*M+5],sz=1; map<string,ll> mp; vector<ll> v[N+5]; void insert(ll &rtt,string s){ if(!rtt) rtt=++sz; ll rt=rtt; for(char ss:s){ if(!tr[rt][ss-'0']) tr[rt][ss-'0']=++sz; rt=tr[rt][ss-'0']; vis[rt]++; } ed[rt]++; return; } ll searchh(ll rt,string s){ ll sum=0,cnt=1,tot=0; for(char ss:s){ // cout<<rt<<"rt\n"; for(ll j=0;j<=9;j++){ if(j==ss-'0') continue; sum+=vis[tr[rt][j]]*cnt; } cnt++; rt=tr[rt][ss-'0']; if(ed[rt]) sum+=(cnt-1)*ed[rt],tot+=ed[rt]; } bool flag=false; for(ll i=0;i<=9;i++){ if(tr[rt][i]){ flag=true; break; } } if(flag) sum+=(vis[rt]-tot)*(cnt-1); return sum; } void work(ll rt,string s){ for(char ss:s){ cout<<rt<<endl; for(ll j=0;j<=9;j++){ if(tr[rt][j]) cout<<j<<" "<<tr[rt][j]<<endl; } rt=tr[rt][ss-'0']; } return; } int main(){ scanf("%lld",&n); ll sx=0; for(ll i=1;i<=n;i++){ bool ff=0; string ss; cin>>ss; for(ll j=0;j<(ll)ss.size();j++){ if(ss[j]=='.') ff=1; else if(ff==0) a[i].x+=ss[j]; else a[i].y+=ss[j]; } if(!mp[a[i].x]) mp[a[i].x]=++sx; v[mp[a[i].x]].push_back(i); // insert(root[0],a[i].y); } for(ll i=1;i<=sx;i++){ for(ll j:v[i]){ insert(root[i],a[j].y); } } // work(root[1],"4514"); // cout<<searchh(root[1],"1145"); // return 0; for(ll i=1;i<=n;i++){ // ll xx=searchh(root[0],a[i].y); ll yy=searchh(root[mp[a[i].x]],a[i].y); printf("%lld\n",yy+n-v[mp[a[i].x]].size()); } return 0; } /* 8 1.1145 1.114 1.1145 1.514 1.19198 1.1154 1.114 1.1145 */ 找错误
最新发布
08-01
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值