提交记录 21475


用户 题目 状态 得分 用时 内存 语言 代码长度
ETO_leader 1002i. 【模板题】多项式乘法 Accepted 100 30.226 ms 8648 KB C++17 2.71 KB
提交时间 评测时间
2024-04-07 22:39:53 2024-04-07 22:39:57
#include<bits/stdc++.h>
#define cir(i,a,b) for(int i=a;i<b;++i)
using namespace std;
using complf=complex<double>;
/*namespace fcomplex{
    template<typename _Ty>
    class fcomplex{
    private:
        _Ty a,b;
    public:
        auto&operator+
    };
}*/
namespace poly{
    static constexpr auto pi=acos(-1);
    auto change(vector<complf>&a,const int n){
        vector<int> rev(n);
        cir(i,0,n){
            rev[i]=(rev[i>>1]>>1)|((n>>1)*(i&1));
        }
        cir(i,0,n) if(i<rev[i]) swap(a[i],a[rev[i]]);
    }
    template<const int type>
    auto fft(vector<complf>&a,const int n){
        change(a,n);
        for(auto h=2;h<n+1;h<<=1){
            const auto wx=complf(cos(pi*2/h),sin(pi*2*type/h));
            const auto midh=h/2;
            for(auto i=0;i<n;i+=h){
                auto u=complf(1,0);
                cir(k,i,i+midh){
                    const auto wy=u*a[k+midh];
                    a[k+midh]=a[k]-wy;a[k]+=wy;
                    u*=wx;
                }
            }
        }
        if(type==-1) for(auto&i:a) i/=n;
    }
    auto mul(vector<int> a,vector<int> b){
        const auto n=1<<((int)(log2(a.size()+b.size()))+1);
        vector<complf> fx(n);
        cir(i,0,(int)(a.size())) fx[i]+=complf(a[i],0);
        cir(i,0,(int)(b.size())) fx[i]+=complf(0,b[i]);
        fft<1>(fx,n);
        for(auto&i:fx) i*=i;
        fft<-1>(fx,n);
        vector<int> ans(n);
        cir(i,0,n) ans[i]=round(fx[i].imag()/2);
        return ans;
    }
}
namespace flash_man{
    static array<char,5000000> buf;
    auto p1 = buf.begin(), p2 = buf.begin();
    #define getchar() p1 == p2 && (p2 = (p1 = buf.begin()) + fread(buf.begin(), 1, 1000000, stdin), p1 == p2) ? EOF : *p1++
    template <typename type>
    inline auto read(type &x) {
        x = 0;
        char c = getchar();
        while (!isdigit(c)) {
            c = getchar();
        }
        while (isdigit(c)) {
            x = (x << 1) + (x << 3) + (c ^ 48);
            c = getchar();
        }
    }
    inline auto read(){
        int x;read(x);
        return x;
    }
    array<char,100> stc;
    inline auto write(int x,char c){
        if(x<0) putchar('-'),x=-x;
        auto tp=0;
        if(!x) stc[tp++]=0;
        while(x) stc[tp++]=x%10,x/=10;
        while((--tp)>-1) putchar(stc[tp]+'0');
        putchar(c);
    }
    #undef getchar
}
int main(){
    //ios::sync_with_stdio(false),cin.tie(0);
    int n,m;
    flash_man::read(n);flash_man::read(m);++n;++m;
    //cin>>n>>m;
    vector<int> a(n),b(m);
    for(auto&i:a) flash_man::read(i);//cin>>i;
    for(auto&i:b) flash_man::read(i);//cin>>i;
    const auto ans=poly::mul(a,b);
    cir(i,0,n+m-1) flash_man::write(ans[i],' ');//cout<<ans[i]<<' ';
    putchar('\n');
    //cout<<'\n';
    return 0;
}

CompilationN/AN/ACompile OKScore: N/A

Subtask #1 Testcase #144.93 us40 KBAcceptedScore: 100

Subtask #1 Testcase #230.059 ms8 MB + 376 KBAcceptedScore: 0

Subtask #1 Testcase #313.09 ms3 MB + 788 KBAcceptedScore: 0

Subtask #1 Testcase #413.107 ms3 MB + 776 KBAcceptedScore: 0

Subtask #1 Testcase #538.98 us40 KBAcceptedScore: 0

Subtask #1 Testcase #637.32 us40 KBAcceptedScore: 0

Subtask #1 Testcase #736.57 us40 KBAcceptedScore: 0

Subtask #1 Testcase #829.169 ms7 MB + 796 KBAcceptedScore: 0

Subtask #1 Testcase #929.147 ms7 MB + 796 KBAcceptedScore: 0

Subtask #1 Testcase #1028.274 ms7 MB + 188 KBAcceptedScore: 0

Subtask #1 Testcase #1130.226 ms8 MB + 456 KBAcceptedScore: 0

Subtask #1 Testcase #1227.147 ms7 MB + 336 KBAcceptedScore: 0

Subtask #1 Testcase #1337.24 us40 KBAcceptedScore: 0


Judge Duck Online | 评测鸭在线
Server Time: 2024-05-04 06:45:07 | Loaded in 1 ms | Server Status
个人娱乐项目,仅供学习交流使用