Codeforces Round 948 (Div. 2) D. XORificator(哈希)

题目

n*m(n*m<=3e5)的矩阵,

实际为t(t<=1e4)组样例,但保证sum n*m不超过3e5

你可以选一行把所有01翻转,问最多可以让多少列只有一个1,然后把你翻转的行输出

思路来源

其实题还挺裸的,教了一下潘老师,看在这题分数还不低的情况下,写一下题解吧

题解

枚举aij,即枚举一下哪列的哪个数字在答案里,这样其他行翻转的状态是唯一确定的,

哈希记录一下,并记录取得答案的aij,每次转移是o(1)的,只要哈希不冲突即可

单哈希是会被卡的,所以可以双哈希或者mt_19937

代码1(mt_19937)

#include <bits/stdc++.h>
#define endl '\n'
#define ll long long
#define vint vector<int>
#define pii pair<int, int>
#define rep(i, a, b) for (int i = int(a); i < int(b); i++)
#define rng(i, a, b) for (int i = int(a); i <= int(b); i++)
#define per(i, a, b) for (int i = int(b) - 1; i >= int(a); i--)
#define gnr(i, a, b) for (int i = int(b); i >= int(a); i--)
using namespace std;
#define MT
#ifdef MT
void flip(int& c)
{
    c=!c;
}
 
const int M = 1e9+7;
 
void solve()
{
    mt19937_64 rnd((unsigned int) chrono::steady_clock::now().time_since_epoch().count());
    int n, m;
    cin >> n >> m;
    vector<string> mp(n);
    vector<uint64_t> ran(n);
    ll t = 1;
    rep(i, 0, n)
    {
        cin >> mp[i];
        ran[i]=rnd();
    }
   
    map<uint64_t, int> cnt;
    pair<int, pair<int,int>> ans = { 0, {0, 0} };
    rep(i, 0, m)
    {
        uint64_t hash=0;
        vint state(n);
        if (mp[0][i] == '0'){
            state[0] = 1;
            hash+=ran[0];
        }
        rep(j, 1, n)
        {
            if (mp[j][i] == '1') {
                state[j] = 1;
                hash+=ran[j];
            }
        }
        ans = max(ans, { ++cnt[hash], {i,0}});
        rep(j, 1, n)
        {
            hash+=(!state[j-1])*ran[j-1]-state[j-1]*ran[j-1];
            flip(state[j - 1]);
            hash+=(!state[j])*ran[j]-state[j]*ran[j];
            flip(state[j]);
            ans = max(ans, { ++cnt[hash], {i,j}});
        }
    }
    cout << ans.first << endl;
    auto [x,y] = ans.second;
    //cout<<x<<':'<<y<<endl;
    rep(j,0,n){
    	bool b = mp[j][x]-'0';
    	if((b&&(j!=y)) || ((!b)&&(j==y))){
    		cout<<1;
    	}else{
    		cout<<0;
    	}
    }
    cout<<endl;
}
#endif
 
int main(void)
{
    std::ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    
#ifdef MT
    int t;
    cin >> t;
    rep(i, 0, t) solve();
#endif
}

代码2(双哈希)

#include <bits/stdc++.h>
#include <iostream>
#include <map>
#define endl '\n'
#define ll long long
#define vint vector<int>
#define pii pair<int, int>
#define rep(i, a, b) for (int i = int(a); i < int(b); i++)
#define rng(i, a, b) for (int i = int(a); i <= int(b); i++)
#define per(i, a, b) for (int i = int(b) - 1; i >= int(a); i--)
#define gnr(i, a, b) for (int i = int(b); i >= int(a); i--)
using namespace std;
#define MT
#ifdef MT

void flip(int& c)
{
    c=!c;
}

const int M = 1e9+7,M2=1e9+9;

void solve()
{
    int n, m;
    cin >> n >> m;
    vector<string> mp(n);
    srand(time(NULL));
    vector<ll> ran(n),ran2(n);
    ll t = 1, t2=1;
    rep(i, 0, n)
    {
        cin >> mp[i];
        ran[i]=rand()*t%M;
        ran2[i]=rand()*t2%M2;
        t<<=1;t2<<=1;
        t%=M;t2%=M2;
    }
   
    map<pair<ll,ll>, int> cnt;
    pair<int, pair<int,int>> ans = { 0, {0, 0} };
    rep(i, 0, m)
    {
        ll hash=0,hs2=0,cur=0;
        vint state(n);
        if (mp[0][i] == '0'){
            state[0] = 1;
            hash+=ran[0];
            hs2+=ran2[0];
            hash%=M;
            hs2%=M2;
        }
        rep(j, 1, n)
        {
            if (mp[j][i] == '1') {
                state[j] = 1;
                hash+=ran[j];
                hs2+=ran2[j];
                hash%=M;
                hs2%=M2;
            }
        }
        ans = max(ans, { ++cnt[make_pair(hash,hs2)], {i,0}});
        rep(j, 1, n)
        {
            hash+=M+(!state[j-1])*ran[j-1]-state[j-1]*ran[j-1];
            hs2+=M2+(!state[j-1])*ran2[j-1]-state[j-1]*ran2[j-1];
            hash%=M;
            hs2%=M2;
            flip(state[j - 1]);
            hash+=M+(!state[j])*ran[j]-state[j]*ran[j];
            hash%=M;
            hs2+=M2+(!state[j])*ran2[j]-state[j]*ran2[j];
            hs2%=M2;
            flip(state[j]);
            ans = max(ans, { ++cnt[make_pair(hash,hs2)], {i,j}});
        }
    }
    cout << ans.first << endl;
    auto [x,y] = ans.second;
    //cout<<x<<':'<<y<<endl;
    rep(j,0,n){
    	bool b = mp[j][x]-'0';
    	if((b&&(j!=y)) || ((!b)&&(j==y))){
    		cout<<1;
    	}else{
    		cout<<0;
    	}
    }
    cout<<endl;
}
#endif

int main(void)
{
    std::ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    
#ifdef MT
    int t;
    cin >> t;
    rep(i, 0, t) solve();
#endif
}

潮汐研究作为海洋科学的关键分支,融合了物理海洋学、地理信息系统及水利工程等多领域知识。TMD2.05.zip是一套基于MATLAB环境开发的潮汐专用分析工具集,为科研人员与工程实践者提供系统化的潮汐建模与计算支持。该工具箱通过模块化设计实现了两大核心功能: 在交互界面设计方面,工具箱构建了图形化操作环境,有效降低了非专业用户的操作门槛。通过预设参数输入模块(涵盖地理坐标、时间序列、测站数据等),用户可自主配置模型运行条件。界面集成数据加载、参数调整、可视化呈现及流程控制等标准化组件,将复杂的数值运算过程转化为可交互的操作流程。 在潮汐预测模块中,工具箱整合了谐波分解法与潮流要素解析法等数学模型。这些算法能够解构潮汐观测数据,识别关键影响要素(包括K1、O1、M2等核心分潮),并生成不同时间尺度的潮汐预报。基于这些模型,研究者可精准推算特定海域的潮位变化周期与振幅特征,为海洋工程建设、港湾规划设计及海洋生态研究提供定量依据。 该工具集在实践中的应用方向包括: - **潮汐动力解析**:通过多站点观测数据比对,揭示区域主导潮汐成分的时空分布规律 - **数值模型构建**:基于历史观测序列建立潮汐动力学模型,实现潮汐现象的数字化重构与预测 - **工程影响量化**:在海岸开发项目中评估人工构筑物对自然潮汐节律的扰动效应 - **极端事件模拟**:建立风暴潮与天文潮耦合模型,提升海洋灾害预警的时空精度 工具箱以"TMD"为主程序包,内含完整的函数库与示例脚本。用户部署后可通过MATLAB平台调用相关模块,参照技术文档完成全流程操作。这套工具集将专业计算能力与人性化操作界面有机结合,形成了从数据输入到成果输出的完整研究链条,显著提升了潮汐研究的工程适用性与科研效率。 资源来源于网络分享,仅用于学习交流使用,请勿用于商业,如有侵权请联系我删除!
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

小衣同学

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值