StkOvflow

STACK OVERFLOW!

一言(ヒトコト)

即得易见平凡,由上自证显然,留作习题答案略,读者自证不难。
反之亦然同理,推论自然成立,略去过程Q.E.D,由上可知证毕。
——lca

AcWing. 1165.单词环

传送门

解题思路

我们一开始可以得出一个建图的思路,对于每个字符串我们把它们当成图中的一个点,然后能“接龙”的字符串之间连一条边,在这张图上跑01分数规划。这样貌似是可以的,但是我们看一下这张图有多么庞大:它的点数N的上限可以达到105,而边数应该是AN2近似于N2也就是差不多1010,这显然是一个光建图就会超时的时间复杂度,我们不可以采用。

稍加思考我们可以得到另一个思路:将能接龙的字符串的首尾2位看做是点,连有向边,有向边的边权是这个字符串的长度。

我们经过观察,可以得到这样一个结论:
在本题中,只要两个字符串的首尾和长度相同,就可以看做是同一个字符串

上述结论的推导:

假设有两个字符串,s1s2,我们的图中(第一个TLE礼包图)中存在这样的一个环:

c1>c2>c3>c4>s1>c5>c6

当我们将s2带进去的时候,由于s2的首尾和s1都是相同的,所以 s2必然可以和 c4以及 c5连接。然后由于点的数量不变,而边的权重又都相同,那自然平均值不变,所以这个环:

c1>c2>c3>c4>s2>c5>c6

和上面的环是等价的。

01分数规划的解题思路

这题要求出一个环,使得这个环(点的个数是pcnt)的wipcnt最大。

一样的用二分的思路,我们想要二分出这个最大值,假设当前二分的值是mid

mid<ans的时候,也就是mid<wipcnt的时候,check函数返回真值

将上式进行一点变形:
mid×pcnt<wi

mid×pcntwi<0

然后它就成了熟悉的负环题了,用SPFA判个环就好了

代码

#include <iostream>
#include <cstring>

using namespace std;

const int N = 1e5 + 10, M = 700;
int h[M], e[N], w[N], ne[N], idx;
double dist[M];
int st[M], ring, n;
string str;

void add(int a, int b, int c) 
{
    e[idx] = b, ne[idx] = h[a], w[idx] = c, h[a] = idx ++ ;
}

void dfs(int cur, double x) 
{
    st[cur] = 1;
    for (int i = h[cur]; ~i; i = ne[i]) 
    {
        int j = e[i]; double W = x - w[i];
        if (dist[j] > dist[cur] + W) 
        {
            dist[j] = dist[cur] + W;
            if (st[j] == 1) {
                ring = true ;
                break ;
            }
            dfs(j, x);
            if (ring) return ;
        }
    }
    st[cur] = 2;
}

bool check(double x) 
{
    ring = 0;
    memset(st, 0, sizeof st);
    memset(dist, 0, sizeof dist);
    
    for (int i = 1; i <= 676; i ++ ) 
    {
        if (st[i]) continue ;
        if (ring) break ;
        dfs(i, x);
    }
    
    return ring;
}

int main() 
{
    while (scanf("%d", &n), n)
    {
        idx = 0;
        memset(h, -1, sizeof h);
        
        for (int i = 1; i <= n; i ++ ) 
        {
            cin >> str;
            int len = str.size();
            int a = (str[0] - 'a') * 26 + str[1] - 'a';
            int b = (str[len - 2] - 'a') * 26 + str[len - 1] - 'a';
            add(a + 1, b + 1, len);
        }
        
        if (!check(0)) puts("No solution");
        else 
        {
            double l = 0, r = 1e5;
            while (r - l > 1e-4) 
            {
                double mid = (l + r) / 2;
                if (check(mid)) l = mid;
                else r = mid;
            }
            printf("%.2lf\n", r);
        }
        
    }
    
    return 0;
}

posted @   StkOvflow  阅读(41)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· 上周热点回顾(2.24-3.2)
点击右上角即可分享
微信分享提示